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(1) From a set of assumptions which seem reasonable 
a priori, the authors calculate the charge distributions in 
the molecules of mono-substituted benzenes, and the 
dipole moments of the latter from reaction rate data alone. 
These coincide very well with the observed moments. 
(2) From the dipole moments of aromatic and aliphatic 
compounds, the ortho-, meta- and para-percentages during 
nitration were calculated. These also accord satisfactorily 
with the experimental values. (3) It is concluded that 
these two facts give a quantitative basis to the usual 
explanation about the orienting power of substituents of 
benzene. (4) The nitration of benzene is an ionic reaction: 
The nitrating agent X- NO, dissociates into X~ and NO.+ 
in the activated complex, and the latter ion interacts with 


the reacting carbon. The activation free energy is reduced 
by the amount of the electrostatic interaction, i.e., 
4.80-10-%,/(rD). Here ¢, is the charge on the reacting 
carbon, 7, the distance between NO,* and the reacting 
carbon in the activated complex, D, the dielectric constant 
of the medium between NO,* and the reacting carbon. (5) 
The authors find that r=1.6A and D=1 approximately, 
independent of the nature of the reacting media. The con- 
sequences of this theory are outlined. (6) Further this 
interpretation of reaction rates provides, for the first time, 
a method of calculating the atomic charges in the molecules 
as distinguished from the usual dielectric measurements 
which lead only to dipole moments. 





1. INTRODUCTION 


T is a well-known fact that some substituents 

of benzene e.g., CH3, OH, NHe, halogens, etc. 
act as ortho- and para-orienting while NO», 
COOH, SO3H, etc. cause reaction in the meta- 
position when mono-substituted benzenes are 
subjected to substitution reactions, e.g., nitra- 
tion, halogenation, sulphonation, etc. This is 
explained in modern organic chemistry by the 
following two assumptions: 

(a) An ortho- and para-orienting substituent 
makes the electron density on the corresponding 
carbon atoms greater than that on the meta- 
atom by inductive and resonance effect, while a 
meta-orienting one acts in the opposite way. 


ete 


* Assistant Professor, Kyoto Imperial University, Japan. 
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(b) Substituting agents are electrophilic,! i.e., 
electron-seeking, because of their positive electric 
nature. Thus it reacts most easily with that 
carbon atom which is surrounded by the greatest 
number of electrons. 

Assuming this explanation to be correct then 
the ortho-, meta- and para-percentages, which 
are proportional to the reaction rates at the 
corresponding carbon atoms, should provide a 
means of determining the charge distribution and 
consequently also permit the calculation of the 
dipole moments of such molecules from rate 
data alone. If dipole moments calculated in this 
way from rates coincide with the observed 
values, then the above qualitative explanation 
for the orienting power of substituents will be 


1 Ingold, J. Chem. Soc. 1120 (1933). 
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shown to have a considerably more substantial 
basis. Also assumption (b) which organic chemists 
make from qualitative considerations only will 
be much more precisely formulated. 

In the present work we have calculated the 
moments of mono-substituted benzenes from the 
ortho-, meta- and para-percentages of the sub- 
stituents after nitration since they have been 
studied more completely than other substitution 
reactions. Afterwards the converse calculation is 
made, i.e., the percentages of the nitrated 
products are calculated from the dipole moments 
alone. As far as we know, this is the first attempt 
at such a calculation, although there is a qualita- 
tive statement regarding the relation between 
ortho, meta- and para-substitution and dipole 
moments, which has been known as ‘‘Sutton’s 
rule.’”? 


2. MopEL ASSUMED FOR THE 
ACTIVATED COMPLEX 


In nitration of mono-substituted benzene, the 
following assumptions are made: 

(1) The nitrating agent is ionized into a 
negative ion and an NO;* which in the activated 
complex lies adjacent to the carbon atom with 
which it is to react. 

(2) A’ substituent of benzene changes the 
electric charge on the six-carbon atoms, as is 
usually assumed, by both an inductive and a 
resonance effect. 

(3) The inductive effect induces equal charges, 


2 Sutton, Proc. Roy. Soc. 133A, 668 (1931), Groves and 
Sugden, J. Chem. Soc. 973 (1935). 


€;, on the six-carbon atoms of benzene and a 
charge —6e; on the inductive group. This charge 
may be either positive or negative. In first 
approximation, we assume further that the 
charge, ¢,, which also may be positive or nega- 
tive, accumulates on each of the two ortho- and 
one para-atoms by resonance with the three 
bond structures in Fig. 1, there being no corre- 
sponding effect at the meta-position. (In Fig. 1 
the curved arrows between X and the substi- 
tuted carbon atoms indicate the directions of 
electron movements.) 

(4) The charge 3e, on the ortho- and para- 
carbons is balanced by a charge of —#e, on the 
substituent group and —3e, on the carbon atom 
to which the latter is bonded. 

Therefore, the charge distribution in a mole- 
cule can be represented by Fig. 2 and the follow- 
ing formulas, combining assumptions (3) and (4), 


(a) €m=€i, 
(b) Ep = €0=Emt€r, 
(c) €1=€m—3(€p—€m), 


(d) €2= —6€,—3(€,—€m), 


where ¢, represents the charge on the carbon 
attached to the substituent group X, and e, the 
charge on X. The location of the other e’s is 
indicated by the subscripts. 

(5) We assume further that the electrostatic 
interaction between the reacting carbon atom 
and NO,* changes the activation free energy 
AF for the nitration of benzene by the amount, 
eyén/(rD). The velocity constant at a carbon 
atom for mono-substituted benzene is then given 


X -6&-#6, 
6-46 


Cte, Fig. 2. 
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by the equation 


kT 
By = K—e7 AF P+ eyeq/ (PD) ET (2) 


where e, and e, are the charges on the reacting 
carbon atom and NOs", respectively, the latter 
being taken as 4.80-10-" e.s.u. 7 is the distance 
between the carbon and NO,* which is always 
taken as 1.6A in this calculation. D is the di- 
electric constant of the medium between NO.* 
and the carbon atom which we take equal to 
unity. The other notation has the usual meaning.’ 


3. CALCULATION OF CHARGE DISTRIBUTIONS IN 
MOLECULES AND OF DIPOLE MOMENTS 


If ky is the rate constant at a carbon atom of 
benzene, the following equation follows immedi- 
ately from Eq. (2), 

ky 
— =e ltyen/ (rD))/kT (3) 
ku 


Here, k,/kw can be obtained from the fractions 
Xo, Xm and x, of ortho-, meta- and para-sub- 
stituents after nitration together with the rela- 
tive nitrated amounts ky’ and k,’ of benzene and 
substituted benzene respectively in a competi- 
tive reaction, since we can write down immedi- 
ately the following relations: 


(a) ko/Ru=3x0(kz'/kun’), 
(b) Rm/Rki=3xm(R2'/kn’), (4)* 
(c) ky/Ru =6xp(R2’/kun’). 


Thus, the charge €m and ¢€, are calculated from 
(3) and (4), and consequently the entire charge 
distribution in the molecule of mono-substituted 
benzene is obtained using Eq. (1). 

The dipole moment of the molecule is given by 
the equation 


m=ex(rz+ro) +(€1—€m)P0, (S) 
where r, is the length of a C—C bond in the 


* Eyring, J. Chem. Phys. 3, 107 (1935), Wynne-Jones 
and Eyring, J. Chem. Phys. 3, 492 (1935). 

* There are six positions for nitration of benzene while 
there are two for ortho- and meta- and one for paranitration 
in the case of substituted benzene. Consequently, ko/ky 
= (xole’/2)/(key’/6), Rm/ k= (xmkz’/2)/(kn'/6) and kp/ky 
=(xpkz’)/(ku’/6). 


benzene nucleus and 7, is the distance between 
the substituent X and the carbon to which it is 
attached. Thus, if we assume e)=€e, as first 
approximation, the Eq. (4a) is not necessary to 
compute dipole moments. The calculations of 
dipole moments to be given next, are based on 
this assumption. However this assumption is not 
strictly correct, as will be discussed later in the 
refinement of the theory, which is necessary to 
take account of the difference in rates in the 
ortho- and para-substitution, i.e., of the differ- 
ence between the charges in ortho- and para- 
position. We shall hereafter use the symbol (eo) 
for the charge obtained from the rate of the 
ortho-position, to denote that it is different from 
the charge obtained from the simple assumption 
€»=€9. When we use (eo) also for the calculation 
of dipole moments it changes our results only by 
about 5 percent, as will be shown later. This 
does not affect significantly the agreement with 
experiment. 

The values 7, and r, used in the calculation are 
tabulated in Table I.4 The r, values in this 
table represent the bond lengths between substi- 
tuted carbon atoms and the atoms in the sub- 
stituents which combined directly with the 
former. 

At first, results will be given, which are 
obtained by making use of the studies of Ingold 
and collaborators on competitive nitration. Gen- 


TABLE I. 7,=1.39A. 








CH:2Cl | CHCl: CCls 


1.50 | 1.50 | 1.50 





























COOCHs | COOC:Hs 
1.54 1.54 


























erally the percentages of the various substitution 
products are independent of the reacting media. 
This is a justification for taking D constant 
independent of the solvent. This result also 
indicates that the electrostatic effect of other 
carbons than the one being attacked do not 
contribute appreciably to the activation energy, 


4These values were taken from Brockway and Beach 


(Rev. Mod. Phys. 8, 231 (1936), Am. Ass. Adv. Sci. No. 7, 
p. 88 (1939)). 
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TABLE II.® Nitration of toluene. pcate= 


0.48; tors = 0.44. 








303°K 
kz'/kyy’ =22 


MEAN FREE ENERGIES AND 
CHARGE DISTRIBUTION 





€yéen 


rD 
(CAL./MOLE) 


€yén 


rD 
(CAL./MOLE) 


eyen 
rD 
(CAL./MOLE) 


€, 71019 
(E.S.U.) 





— 2071 
— 497 
— 2247 


— 2178 
— 611 
— 2354 


—2125 
— 554 
— 2301 


(—0.0491)* 
—0.0128 
—0.0532 




















0.0478 
0.137 

















* These numbers within the parentheses in this and the following tables are the charges (eo) of the ortho-positions calculated from the rates. 


TABLE III.° Nitration of ethylbenzoate. 
Keale = —1. oi Mobs = —1.91. 








° 


291°K 
kz'/K jy’ =0.00367. 





€yéen 
€y °1010 
(E.S.U.) 


(0.0772) 
0.0649 
0.0950 
0.0198 

—0.430 


rD 
(CAL./MOLE) 


3338 
2809 
4109 


k,y/Ky 


0.00303 
0.00759 
0.000793 























as such an interaction would be modified by the 
dielectric constant of the media. Because of this 
constancy we followed Bird and Ingold® in using 
mean percentages obtained from experiments in 
different media and at different temperatures. 
The method of averaging of the ortho-, meta- 
and para-percentages obtained under different 
conditions was as follows: By the equations 
te, way (6) 
rD Xm 


AF.t= ((€0) —em)en _ 


Xo 
—kT |In — (7) 
1D Xm 


5 Bird and Ingold, J. Chem. Soc. 918 (1938). 


obtained from Eqs. (2) and (4), the mean free 
energies AF ,,.t and AF,, are averaged regardless 
of the conditions under which the nitration was 
conducted. If the mean free energies are then 
put back into Eqs. (6) and (7), one can now 
calculate ‘‘mean percentages’ at any desired 
temperature. This method is the same in prin- 
ciple as that used by Scheffer,® by Bradfield and 
Jones,’ and by Bird and Ingold.* Ingold and col- 
laborators showed also that the rate ratios for 
substituted benzenes to benzene for nitration, 
i.e., kz’/Ry’ are independent of the reacting 
media. Thus mean values of k,’/ky’ in different 
media were used in our calculations. 

In the following Tables II, III, IV, V, VI and 
VII, the percentages were obtained by the above 
procedure making use of Table III of the paper 
of Bird and Ingold.* The rate ratios, k.’/ku’, 
were taken from the data of the literature as 
referred to in the following tables. 

To simplify our calculation of dipole moments, 
we have taken the total charge, ez, on the sub- 
stituent X as being located on that atom attached 
to the carbon atom 1 of the benzene ring. 

6 Scheffer, Proc. K. Akad. Wetensch. Amsterdam, 15, 


118 (1913). 
7 Bradfield and Jones, J. Chem. Soc. 1006, 3073 (1928). 


TABLE IV.® Nitration of chlorobenzene. pcate= — 1.54; bors = — 1.55. 








° 


273°K 2 
kz'/kyy’ =0.0269 ke’ /kyy’ =0.0322 


08°K MEAN FREE ENERGIES AND 
kz'/k yy’ =0.0364 CHARGE DISTRIBUTION 





€yén €yén 


rD 
ky/kH | (cAL./MoLE) 


rD 
ky/ky | (cat./MoLe) 


€yen €yéen 
rD 
(CAL./MOLE) 


rD 
ky/ky (CAL./MOLE) 





0.0228 
0.000164 
0.116 


2041 
4705 
1163 


0.0300 
0.000193 
0.133 























0.0351 
0.000218 
0.148 





























[Q10 
.U.) 


0491)* 
0128 
0532 
0478 
137 
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TABLE V.° Nitration of bromobenzene. peate= — 1.68; poss = — 1.53. 








° ° 


2 
kz! /k jy’ =0.0246 kz! /kyy’ =0.0305 


308°K MEAN FREE ENERGIES AND 
kz’ /kyy’ =0.0321 CHARGE DISTRIBUTION 



































€yén €yén €yéen €yéen 
: 7D 7D 7D YD €, *1010 
% ky/ky (CAL./MOLE) % ky/ky (CAL./MOLE) % ky/ky (CAL./MOLE) (CAL./MOLE) (E.S.U.) 
o | 39.0 |0.0288 1915 4.14 | 0.0379 1928 42.4 0.0405 1954 1932 (0.0446) 
m | (0.2) |0.000148 4761 (0.2) | 0.000183 5072 (0.2) 0.000193 5211 5015 0.116 
Pp | 61.0 | 0.0900 1300 58.6 | 0.107 1317 57.6 0.110 1345 1321 0.0305 
1 0.244 
X —0.568 




















The results for toluene and ethylbenzoate are 
given in Tables II® and III,® respectively, to- 
gether with k,’/ky’ used in the calculations. The 
agreement here is satisfactory. 

The dipole moment —1.31 calculated in this 
case is that due to resonance and induction in 
the benzene ring. The extra 0.60 Debye unit 
observed is to be attributed to the rest of the 
molecule. 

Next we shall give the results for halobenzenes. 
Usually halobenzenes are assumed not to pro- 
duce meta-substituted products. But this is not 
strictly true. Thus Holleman (see Bird and 
Ingold'®) found 0.3 percent of the meta-product 
in the nitration of bromobenzene. The authors 
have taken 0.2 percent for the amount of meta- 
products which lies within the usual experimental 
error. This value is used in Tables IV, V, VI 
and VII for the calculation. 

The agreement between the calculations and 
experiments in these cases also is very good. 
Here, however, it hinges on the approximate 
correctness of our estimation of the meta-per- 
centages. We also found, as will be shown in a 
later paper, the same e» from a study of acid 
and base-strengths and from side chain reactivi- 
ties of meta-halogenated compounds. Therefore, 
our estimation of the meta-products is quite 


- reasonable. 


For the compounds which have not yet been 
competitively nitrated in the presence of benzene 
the dipole moments were calculated by the 
following procedure. 





§Ingold, Lapworth, Rothstein and Ward, J. Chem. 
Soc. 1959 (1931); Bird and Ingold, ibid. 918 (1938). 
The dipole moments in Debye units in this and the follow- 
ing tables, are the mean values of the data in Trans. 
Faraday Soc. 30, Appendix (1934). The signs for the 
direction of the dipole moments are the customary ones. 

* Ingold and Smith, J. Chem. Soc. 905 (1938). 

” Reference 5, p. 923. 


TABLE VI.5* Nitration of iodobenzene. 
Keale = —1.27; Mobs >= ~ 1.30. 








° 


2 
kz'/ky’ =0.18 








€yéen 
7D €y *1010 
% ky/ky (CAL./MOLE) (E.S.U.) 
o | 43.1 0.232 841 (0.0194) 
m | (0.2) 0.00108 3931 0.0909 
Pp | 56.9 0.613 282 0.00651 
1 0.218 
aA —0.419 




















* In the competitive nitration of iodo- and fluorobenzene, there are 
considerable errors in kz’/k 7’. Especially in the case of iodobenzene, the 


values of kz’/kz7’ change within the range 0.1-0.3 with no distinction 
for the temperatures between 0° and 35°, under which the experiments 
were carried out. Therefore, in the calculation of these cases, we used 
the mean values of kz’/ky’ at 18°, which were given by these authors. 
Since they have not determined the ortho- meta- and para-percentages 
of the substituents of these compounds, we used the mean percentages 
at 18° obtained from Holleman’s values (Chem. Rey. 1, 187 (1924)). 


TABLE VII. Nitration of fluorobenzene. 
Keale = — 1.20; hobs = — 1.45, 








ov 


2 
kz'/kyy’ =0.15 








€yéen 
7D €y °10'0 
% ky/ky (CAL./MOLE) (E.S.U.) 

o | 14.6 0.0653 1570 (0.0363) 
m | (0.2) 0.00090 4036 0.0933 
p | 85.4 0.768 152 0.00351 
1 0.181 
X —0.473 




















First, Eq. (6) is used, where 2x»/xm is known 
experimentally. Therefore, if ¢» of a molecule is 
known by another method, e, can be calculated. 
We assume that the inductive part of the dipole 
moment of mono-substituted benzene is equal 
to the moment wai of the corresponding substi- 
tuted tertiary aliphatic compound, i.e., 


— 6€m(ro+7z) = HMali- (8) 


Equation (8) permits us to calculate ¢,, from 
Mali, and consequently we can now calculate the 
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TABLE VIII." WNitration of nitrobenzene. pari= —3.29 
(nitro-n-butane). weate= — 3.63; ors = — 3.81. 

















MEAN FREE ENERGIES 
AND CHARGE 
243°K 303°K DISTRIBUTION 
(ey —em)en (€y—em)en | (ey —em)en 
rD rD rD €y °1010 
% (CAL./MOLE)| % (CAL./MOLE)| (CAL./MOLE)]| (E.S.U.) 
o 6.4 1287 8.2 1443 1365 (0.223) 
m | 93.2 91.1 0.191 
7) 0.25 2514 0.7 2438 2476 0.248 
1 0.105 
xX —1.23 


























whole charge distribution in a molecule. The 
following results were obtained by this procedure. 
(When there is no datum for the tertiary com- 
pound, the value for a compound which has a 
structure as similar as possible to the latter is 
used.) In Tables VIII—XIII,"-"* the charges in 
the ortho-positions were calculated by Eq. (7). 

In these cases also, the calculated dipole 
moments accord very well with the observed 
values. 

We applied this method to the calculation of 
the dipole moments of ethylbenzoate using 
Maii= —1.78 for ethylbutyrate and the ortho-, 
meta- and para-percentages in Table III. The 
charges (€0), €m, €p, €1 and e, are (0.114), 0.101, 
0.131, 0.0562 and —0.653, respectively, and 
Mcal= —1.98. Considering this excellent agree- 
ment between the calculated and observed 
moment (obs=—1.91), the explanation con- 
cerning the discrepancy between the latter and 
the calculated moment from the rate data alone 
is apparently valid. 


11 Holleman and Buryn, Rec. trav. chim. 19, 79 (1900). 

12 Holleman, Rec. trav. chim. 18, 265 (1899). 

13 Holleman, Direkte einfuhrung von Substituenten in der 
Benzolkernen (1910), p. 125. 

144 Holleman, Vermeulen and de Mooy, Rec. trav. chim. 
33, 1 (1914). 


TABLE IX.” Nitration of benzoic acid. paii= —0.73 (butyric acid). peatc= —1.02; os = — 1.0. 
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4. REFINEMENT OF THE THEORY 


In the above calculations, we assumed €o= €>. 
But actually because the charge in the ortho- 
position is in closer proximity to the charge e€; on 
the substituted carbon, the corresponding struc- 
ture will occur to a somewhat different extent 
than does the para-structure. The quantum- 
mechanical consideration of this will be given 
later, when the theory of the charge distribution 
in the molecule of substituted benzene, will be 
discussed."5 Here we shall use the empirical 
formula 

€1—€p 


~ (9) 
15 





(€0) =ep+ 


or 


(€o) =€pt ’ 
8 





(9’) 


to correct the charges on the orthocarbon atoms 
of the molecule, when it is necessary to obtain 
the actual charges (€) on the orthocarbons from 
the charges e,. (Both Eqs. (9) and (9’) work 
almost equally well. We will consider which is 


18 Wheland and Pauling, J. Am. Chem. Soc. 57, 2086 
(1935); E. Hiickel, Zeits. f. Physik 72, 310 (1931); E. 
Hiickel and W. Hiickel, Nature 129, 937 (1932); E. Hiickel, 
Zeits f. physik. Chemie B35, 163 (1937). 

















TABLE X.%3 WNitration of methylbenzoate. pari= — 1.05 
(methylvalerate). tcate= — 1.83; tops = — 1.86. 
MEAN FREE ENERGIES 
AND CHARGE 
243°K 303°K DISTRIBUTION 
(ey —em)en (ey—em)en | (ey —em)en 
rD rD ., rD €y “101° 
% (CAL./MOLE)| % _ | (CAL./MOLE)] (CAL./MOLE)| (E.S.U.) 
°o 23.6 552 25.7 599 576 (0.1062) 
m | 74.4 69.8 0.0929 
Dp 2.0 1405 4.5 1227 1316 0.123 
1 0.0473 
XxX 0.603 




































273°K 


MEAN FREE ENERGIES AND 


303°K CHARGE DISTRIBUTION 



































(ey —€m)en (ey —em)en (ey —em)en (ey —em)en 
rD rD rD rD €, °10° 
% (CAL./MOLE) % (CAL./MOLE) % (CAL./MOLE) (CAL./MOLE) (E.S.U.) 
0 14.4 855 18.8 784 22.3 738 792 (0.0594) 
m 85.0 80.2 76.5 0.0411 
p 0.6 2049 1.3 1852 1.2 2074 1992 0.0871 
1 ‘ —0.0342 
A —0.316 
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TABLE XI." WNitration of benzylchloride. pai=—1.98 
(1-chloro-2-methylpropane). tcaic= — 1.69; ors = — 1.82. 


TABLE XII." WNitration of benzalchloride. wau=—2.06 
(1,1-dichloropropane). pcate= — 1.95; ors = — 2.04. 














298°K 298°K 
(ey —em)en (ey —em)en 
rD €y 101° rD €y °101° 
% (CAL./MOLE) (E.S.U.) % (CAL./MOLE) (E.S.U.) 
0 40.9 1342 (0.083) 0 23.3 220 (0.122) 
m 4.2 0.114 m | 33.8 0.117 
p | 54.9 — 1924 0.070 p | 42.9 — 549 0.104 
1 0.180 1 0.136 
xX —0.617 X — 0.683 


























best from a theoretical point of view at a later 
date.) 

We assume the difference in charge between e, 
and (€9) accumulates on the substituted carbon 1 
and substituting group X equally, according to 
our assumption (4) made in the model for the 
activated complex. Then the corrected charges 
on these two atoms will be given by the following 
equations: 

(€1) =e. +€p— (€0), 
(€2) = er tep— (€0). 


Here, parentheses around the symbol for charge 
indicate the charge has been corrected in the 
way which follows from Eqs. (9) or (9’) and (10). 
The corrected moments together with (¢€),* (1) 
and (e,) are given in Table XIV. 

From Table XIV, it will be concluded that this 
correction often makes the calculated dipole 
moments slightly worse. The difference is usually 
of the order of 5 percent, and consequently it 
does not affect significantly our previous con- 
clusions. A slight change in the value 1.6A for 
the distance between NO.* and the carbon being 
substituted would reduce the discrepancy be- 
tween the calculated and observed dipole mo- 
ments. 

Because the use of the refined theory to 
calculate the ‘‘refined” dipole moments does not 
improve them but, if anything, makes them 
slightly worse, it is interesting to note that one 
can assume the “approximate” calculation of 
dipole moments is correct. In this case one 
interprets the difference between the ortho- and 
Para-percentages not as due to different charges 
at these positions but as due to the difference in 
charge in the two neighboring positions. Thus a 


(10) 


* The charges (€) in Table XIV are the values obtained 
from rates, and not the values obtained by Eq. (9) or (9’). 


TABLE XIII." Nitration of phenylchloroform. pati = — 1.54 
(methychloroform). pcatc= — 1.48; bors = —2.11. (This value 
ts lower than the observed value, probably because of the lack 
of a satisfactory value for pati.) 














298°K 
(ey —em)en 

rD €y +1010 

% (CAL./MOLE) (E.S.U.) 
0 6.8 1328 (0.118) 
m | 64.5 0.0876 
p 28.7 — 68.7 0.0892 
1 0.0900 

xX —0.524 

















para-position is flanked by two meta-positions 
while an ortho-position is flanked by a meta- 
position on one side and ¢; on the other. If the 
distance, r’, from the charge e€, to the substi- 
tuting group and D’, the effective dielectric 
constant, satisfies the equation r’D’=8X1.6, 
then the free energy of activation in the ortho- 
position will differ from that in the para-position 
by the amount, (€1—€m)en/(8X1.6), or we can 
say that the apparent charge in the ortho-position 
is given by Eq. (9’). Actually this interpre- 
tation leaves us with our original “‘approximate’”’ 
dipole moments. One can, however, only justify 
neglecting the effect of the larger charge on the 
original substituent in the benzene on the basis 
of a larger effective dielectric constant and a 
somewhat greater distance. This assumption 
that the flanking charges contribute in an im- 
portant way to the free energy of activation 
would necessitate other obvious changes in the 
calculation of the dipole moments. The results 
are in less satisfactory agreement with experi- 
ment, so that we prefer the assumption that 
there is a different charge in the ortho- and para- 
position. 








































































T. RI AND H. EYRING 
TABLE XIV. 

CeHsCHs CeHsCOOC2Hs CsHsCl CeHsBr CeHsI CcHsF 
(€o) -10%° —0.0491 0.0772 0.0474 0.0446 0.0194 0.0363 
(e1) - 101° 0.0437 0.0376 0.227 0.230 0.205 0.148 
(ez) + 101° 0.133 —0.412 —0.578 — 0.582 — 0.432 — 0.506 
Brefined 0.46 —1.25 — 1.62 —1.74 — 1.33 — 1.34 
Hobs 0.44 —1.91 — 1.56 —1.53 — 1.30 —1.45 
Mapprox 0.48 —_ 1.31 _= 1.54 = 1.68 _= 1.27 = 1.20 

CeHsNOz2 CsHsCOOH CsHsCOOCH3 CeHsCH:Cl CeHsCHCle CsHsCCls 
(eo) - 101° 0.223 0.0594 0.106 0.083 0.122 0.118 
(€1)- 101° 0.130 —0.0065 0.0643 0.163 0.118 0.0610 
(ez) - 101° — 1.20 — 0.288 — 0.586 — 0.630 —0.701 —0.553 
Prefined —3.50 —0.910 — 1.76 — 1.76 — 2.03 — 1.64 
Mobs —3.81 — 1.00 — 1.86 — 1.82 — 2.04 —2.11 
“andl —3.63 —1.02 — 1.83 —1.69 —1.95 —1.48 


























5. CALCULATION OF ORTHO-, META- AND PARA- 
PERCENTAGES OF NITRATED PRODUCTS 
FROM DIPOLE MOMENTS ALONE 


From the above results, it is clear that there 
is the possibility of calculation of ortho-, meta- 
and para-percentages of nitrated products from 
the dipole moments of aromatic and aliphatic 
compounds. To do this, the authors assumed 
that paro—Mali is equal to the resonance part of 
Maroy i.e., 


(11) 


Eq. (11) is easily obtained from the charge 
distribution in Fig. 2. Combining (11) and (8), 
we can calculate «, and €»,, consequently, the 
whole charge distribution in a molecule follows 
from a knowledge of waii and paro. However, it 
must be kept in mind that in this case the charge 
(€9) on the ortho-atom must be used, which is 
obtained from (9) or (9’). Introducing the 
charges (€o), €m and e, in Eqs. (6) and (7), the 
fractions xo, xm and x, are calculated. In the 
following, the results are given together with 
Mali ANd paro— Mali USed in the calculation.* 

From the above results, it can be understood 
that the agreement between the calculated and 
observed percentages is quite good considering 
that there is a considerable percentage of error 
iN Maro— Mali, due to its small value. Especially, 
in the case of nitrobenzene, the ortho-percentage 
is smaller than the para-percentage in the 
calculation, while the opposite is true experi- 


Maro — Mali = —3(e,— Em) (%2+2rp). 


* The charges (eo) in the following tables were obtained 
using Eq. (9). 





mentally. This is, of course, due to the error in 
Maro—Mali in this case. We did not try this 
calculation in the other cases, partly because of 
the lack of the values of aii as in the cases of 
fluorobenzene and toluene, and partly because 
of the small values of paro—Mali, Which are 
between 0.1—0.2, and consequently contain large 
percentage errors. 

It is noteworthy here that the e,’s of the 
chloro-, bromo- and iodobenzene in Tables XV, 
XVI and XVII accord very well with those 
given in Tables IV, V and VI, respectively, 
which were calculated making the assumption 
that the percentages of meta-nitrated products 
of these compounds is in all cases 0.2 percent. 
This agreement indicates this assumption is not 
unreasonable. 


6. DISCUSSION 


We have shown above that the correct dipole 
moments of substituted benzenes can be calcu- 
lated from the ortho-, meta- and para-percent- 
ages of the substituents after nitration, and that 
conversely the latter can be calculated from the 
dipole moments alone. This agreement shows the 
substantial correctness of our assumptions. How- 
ever, we will now discuss them in more detail. 

Benford and Ingold'* studied kinetically the 
nitration of benzene and substituted benzenes 
in CH;NOsz solution, and found that this reaction 
is zero order. They explained this fact by the 
assumption that the rate determining step is (b) 


16 Benford and Ingold, J. Chem. Soc. 929 (1938). 
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DIPOLE MOMENTS AND NITRATION RATES 





TABLE XV. Nitration of chlorobenzene. pari = —2.15 (tertiary chlorobutane). paro— ati = 0.59. 








(ey —em)en 


273°K 


298°K 
















































































D J or oO 
€, °101° (cat/mote) (%)carc (%) ops (% care (%) ons 
0 (0.0429) — 3160 37.4 28.2 39.8 31.0 
m 0.116 0.1 0.2 
p 0.0282 — 3810 62.5 71.8 60.0 69.0 
1 0.248 
X —0.564 
TABLE XVI. Nitration of bromobenzene. pari = —2.15 (tertiary bromobutane). paro— Mati = 0.62. 
273°K 298°K 
(ey —em)en 
D - - , 
€y °101° ecas. Amand (%)carc (%) ons (%o)carc (%)ons 
0 (0.0389) — 3203 37.4 39.0 39.8 41.4 
m 0.113 0.1 0.2 
p 0.0240 — 3853 62.5 61.0 60.0 58.6 
1 0.247 
X —0.544 . 
TABLE XVII. Nitration of todobenzene. pati = —1.95.* pbaro— Mati =0.65. 
273°K 298°K 
(ey —em)en 
rD Oo oO oy, or 
€y *1010 (CAL./MOLE) 7o)CALC (%)oss (%o)carc (%) ons 
0 (0.0199) — 3229 Ky Pe 41.7 40.0 43.7 
m 0.0945 0.1 0.2 
p 0.0050 — 3874 62.2 58.3 59.8 56.3 
1 0.229 
X —0.433 


























* The dipole moment of tertiary iodobutane is —2.13, and the dipole moment of phenyliodide is —1.30, so warg —Hajj =0-83. 


charges em and ep of CeHslI, using this ugjj, and warg —Hali, then ep becomes negative. 
we know that the nitration of iodobenzene is slower than that of benzene. 


If we calculate the 


This is in disagreement with the experimental fact, because 
Because of the uncertainty of the dipole moment of tertiary buty] iodide, 
we used —1.95 for ua3;, which is the mean value of the dipole moments of the latter and that of 1-iodobutane (u = —1.76). 


TABLE XVIII. Nitration of nitrobenzene. wari = —3.29 (nitro-n-butane). paro— Mali = —0.52. 






































273°K 298°K 
(ey —em)en 
“7D 
€y °1010 (cat./oLe) (Jocare (%)ors (%)carc (%)ors 

0 (0.260) 2944 0.4 6.4 0.7 7.9 

m 0.192 96.0 93.4 94.7 91.6 

p 0.274 3550 3.6 0.2 4.6 0.5 

1 0.069 

X — 1.28 
in the following series of reactions: They concluded also qualitatively, having 
CH;NO.—CH, : NO-OH (a) studied solvent effects upon the nitration, that 


CH, : NO-OH+HO-NO, 
—CH, ° NO-O-NO.+H.0, (b) 
CH: : NO-O-NO.+C,Hs 
—C,H;-NO.+CH:2 : NO-OH. (c) 


the efficiency of the nitrating agents X-NO, 
(X=CH2 : NO-O-— in the above case) depends — 
upon the electron affinities of the X groups for 
the electrons of the bonds X — NOs, as indicated, 
for instance, by acid strength of X-H. Thus 


T. RI AND H. EYRING 


Kékulé structures 


SO.H - NOzis more reactive than CH3;CO-O-NOsz, 
and the latter is more reactive than HO- NOs. 

Basing our results on these experimental facts, 
we assumed that X-NOz dissociates into X~ 
and NO," in the activated complex (assumption 
(1)), and that the latter acts on the charge e, of 
the reacting carbon atoms, reducing the activa- 
tion free energy by the amount of the electro- 
static interaction, e,en/(rD) (assumption (5)). 
These assumptions, therefore, seem reasonable 
a priort. However, some remarks must be made 
on the values of 7 and D, which were taken as 
1.6A and unity, respectively, in our calculations. 

Calculations of potential surfaces for activated 
complexes indicate that those atomic distances 
which are changing are often around 10 percent 
larger than the normal distances. Therefore we 
have chosen 1.6A for the distance C— NO,* in the 
activated complex, since the normal distance 
C—NO, is 1.46A. Our justification for taking 
D constant and independent of solvents has 
already been given, i.e., the constancy of ortho-, 
meta- and para-percentages regardless of the 
reacting media. In general, the macroscopic 
dielectric constants measured by the usual 
experimental method have no meaning for 
phenomena occurring at atomic distances, since 
most of the interatomic lines of force do not 
pass through the solvent. Actually, Debye" takes 
D equal to 2 or 3 independent of solvent when 
he discusses phenomena at molecular distances. 
In our case, moreover, there is no solvent 
molecule between the reacting carbon and NO2t 
in the activated complex. Therefore, it is 
reasonable to take D equal to unity. The 
agreement between the moments calculated and 
observed shows these values are highly satis- 
factory. Interactions between non-neighboring 
atoms will be smaller both because of the 
distance and because of the dielectric constant 
and in this approximation are neglected. 


17 Debye, Dipolmoment und chemische Struktur (Leipzig, 
1929); Polar Molecules (1929). 


Thus we can conclude from the agreement 
between the experiments and our results that 
nitration is an ionic process, since these calcula- 
tions rest on assumption (1) and (5). We have also 
found that organic reactions such as esterification, 
saponification (hydrolysis), benzoylation, and 
quaternary salt formation can be treated in the 
same way from the study of side chain reac- 
tivities. This material will be presented for 
publication soon. Thus all these reactions are 
ionic as organic chemists usually assume. That 
there is a close correlation between the reactions 
considered here and many other processes is 
shown for example by the interesting results of 
Hammett." Interesting but quite different treat- 
ments of these questions are also to be found in 
the recent reports of Jenkins" and of Westheimer 
and Shookhoff.”° 

We assume the same inductive effect on the six 
carbon atoms of the benzene nucleus (assumption 
(3)). This is based on findings that can be 
discussed in sufficient detail only in a subsequent 
paper. However, a double bond in a side chain 
of benzene ring acts as a conductor as shown by 
the induction of charge from the ring on atoms 
in the side chains. By resonance with the two 
Kékulé bond structures in Fig. 3, all the C—C 
bonds partake of the nature of double bonds as 
is evidenced by the fact that the length of a 
C—C bond in benzene is longer than that of an 
aliphatic C=C bond, but shorter than that of 
an aliphatic C—C bond. For this reason we 
assume the benzene ring acts as a conductor in 
the induction, i.e., the charges ¢; produced by 
induction on each of the six carbon atoms are 
equal to each other. We assumed further as a 
first approximation that the same charges «, 
accumulate on each of the two ortho- and one 
para-atoms by resonance. The approximate 
correctness of this assumption has already been 
discussed, as well as a refinement of this simple 
assumption. Although this refinement does not 
effect significantly the calculation of dipole 
moments from rates, it is important in the 
converse calculation, i.e., in the calculation of 


18 Hammett, J. Am. Chem. Soc. 59, 96 (1937); Trans. 
Faraday Soc. 34, 156 (1938). 

19 Jenkins, J. Chem. Soc. 1780 (1939). 

20 Westheimer and Shookhoff, J. Am. Chem. Soc. 62, 
269 (1940). 





DIPOLE MOMENTS AND NITRATION RATES 


TABLE XIX. 








SUBSTITUENTS 
CHARGE 


F 


cl 


Br 


I 





€m* 10!° 
(€p—€m) “102 


0.0933 (n)* 
—0.090 (n) 


0.115 (n) 
—0.088 (n) 


0.116 (n) 
—0.085 (n) 


0.0909 (x) 
—0.084 (n) 





SUBSTITUENTS 
CHARGE 


COOH 


COOCHs 


COOC?2Hs 


NO2z 





€m* 101° 
(€p—€m) + 102 


0.0411 (m)* 
0.0460 (n) 


0.0929 (m) 
0.0300 (n) 


0.0649 (7) 
0.0301 (x) 


0.191 (m) 
0.057 (n) 





SUBSTITUENTS 
CHARGE 


CHs 


CH:Cl 


CHCl. 


CCls 





€m* 101° 
(€p—€m) ” 101° 





—0.0128 (n) 
—0.0404 (n) 





0.114 (m) 
—0.044 (n) 





0.117 (m) 
—0.013 (n) 





0.0876 (m) 
0.0016 (n) 








* Here (m) and (m) represent the values obtained from the nitration and the dipole moments, respectively. 


ortho-, meta- and para-percentages from dipole 
moments alone. 

It is not likely that our simple rules for charge 
distribution in terms of €» and e, will remain un- 
changed as we accumulate more experimental evi- 
dence and improve the theory. However, in their 
present unfinished form these rules are surpris- 
ingly successful. 

We collect the values of €,n and €,—é€m in 
Table XIX. Here, as stated above, e, represents 
the charge induced on each of the carbon atoms 
of the benzene nucleus, and ¢e,—€m represents 
the charge accumulated on one of the ortho- 
or para-carbons by the resonance effect. 

€m and €p—é€m in Table XIX, therefore, give 
the inductive and resonance effects in the 
substituted benzenes both as regards direction 
and magnitude. Thus toluene is the only com- 
pound in Table XIX which accumulates negative 
charge in the benzene ring both by the inductive 
and resonance effects, and it has the greatest 
reactivity for nitration. 

(Here the negative signs in this table show 
the charges are negative; they must not be 


confused with the usual positive and negative 
notations indicating the directions of the induc- 
tive and resonance effects, since positive induc- 
tion, for example, means that negative charge is 
induced in the benzene nucleus as in toluene, 
and vice versa. Thus in halobenzenes, the 
inductions are negative, and resonance positive, 
while in COOH, COOCH;, COOC:H; and NO» 
groups they are both negative.) 

It is very interesting to notice that €,—é€n<0 
in the case of ortho- and para-orienting groups, 
and €p,—€m>0O in the case of meta-orienting 
groups. These relations show clearly that the 
former groups make the electron densities on 
ortho- and para-atoms greater than those on 
meta-atoms, while the latter act in the opposite 
way. Thus, we can see that the usual explanation 
on the orienting power of substituents is reason- 
able, being established by the quantitative 
evidence given above. Our theory, however, 
goes beyond such qualitative pictures to a 
definite model which permits us to calculate 
quantitative results. 












































JUNE, 1940 


JOURNAL OF CHEMICAL PHYSICS 





VOLUME 8 


Far Ultraviolet Absorption Spectra of Some Aliphatic Ketones 
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Absorption spectra of methyl n-propyl, methyl isopropyl and diethyl ketones between 2000 
and 1500A were found. The electronic transitions and their vibrational structure are discussed 
in comparison with acetone and methyl ethyl ketone. It is concluded that most transitions are 


probably of a Rydberg type, leading to ionization potentials in the neighborhood of 10 ev. 





HE absorption spectra of acetone and 

methyl ethyl ketone in the Schumann 
region have been studied in some detail,! but 
spectra of the higher ketones have not been 
reported. It should be expected that the general 
features of the spectra in this region would be 
the same for all ketones, but it seemed of interest 
to study some of the details of the spectra of 
diethyl, methyl isopropyl and methyl n-propyl 
ketones. The vibration frequencies? in the normal 
electronic state, as well as other physical con- 
stants remain remarkably constant as we go 
from acetone to these ketones. 


APPARATUS AND EXPERIMENTAL PROCEDURE 


The first plates were taken in a one-meter 
spectrograph at Brown University.'* The source 
was the usual Lyman continuum, and the 
ketone vapors in the spectrograph were separated 
from the light source only by the narrow slit. 
In this way the wave-length region from 2200 to 
800 was covered. The pressure of the vapors 
was very uncertain and varied during an 
exposure, due probably to absorption by the 
wax gaskets of the spectrograph. This resulted 
in loss of definition in the vibrational structure. 

The final plates were made with a normal 
incidence vacuum spectrograph built here re- 
cently by the writer. It contains a 2-meter 
radius 120,000-line grating, ruled on glass by 
R. W. Wood. The dispersion in the first-order 
ultraviolet is about 4.145A/mm. The light source 
was an uncondensed discharge in hydrogen 
operated usually with 1 amp. at 4000 volts. 
~ 1 (a) Noyes, Duncan and Manning, J. Chem. Phys. 2, 
717 (1934). (b) A. B. F. Duncan, J. Chem. Phys. 3, 131 
(1935). (c) Duncan, Ells and Noyes, J. Am. Chem. Soc. 
58, 1454 (1936). 


2 Kohlrausch and Koppl, Zeits. f. physik. Chemie B24, 
370 (1934). 
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Exposure times were 20-30 minutes with East- 
man III-—O ultraviolet sensitive plates. In most 
cases a fused quartz absorption cell 20 cm long 
was used between source and slit. The plane 
polished windows which were fused on trans- 
mitted to 1560A. In other cases (especially with 
diethyl ketone) the 4-meter column of the 
spectrograph tube was used. To avoid possible 
decomposition products during exposure the 
vapor could be passed through the cell in a 
current of pure nitrogen at different fixed rates. 
Pressures were fixed by immersion of the liquid 
ketone in a CO,-ether hand-regulated thermostat 
or in ice-salt mixtures in which the nitrogen was 
also precooled. The actual vapor pressures were 
measured, and the useful ranges were 0.3—2.5 mm 
for the cell and 0.0025-—0.5 for the longer columns. 
These ranges were covered systematically with 
a large number of exposures. 

The ketones, obtained from Eastman Kodak, 
all gave tests for acids, oxidizing agents (with 
Fe+++CNS_), and biacetyl. A negative chemical 
test showed that the latter could not have a 
partial pressure greater than 210-5 mm, but 
the four strongest bands* were found on the 
first plates. The ketones were washed with 
dilute NaHSO; and Na2CQ; solutions, then with 
many small portions of water. After drying they 
were fractionated in a current of nitrogen. All 
samples used boiled within 0.2° of the reported 
boiling point and showed none of the above 
impurity. While biacetyl is probably exceptional 
in its large absorption coefficient in this region, 
this shows that experimental data on the spectra 
of organic compounds must be examined rather 
critically. 

Lines of NI at 1745.246 and 1742.734A 


3V. R. Ells, J. Am. Chem. Soc. 60, 1864 (1938). 
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appeared on most plates and were used as 
primary standards. In addition Sil 1989.01 and 
CI 1930.900 * appeared on all plates. From these 
a few additional lines of silicon were identified. 


RESULTS 


Although photographs were made down to 
800A for the three ketones over considerable 
range of pressure, nothing but general absorption 
could be found below 1500A. Not even maxima 
which might be correlated with high electronic 
terms were apparent. Between 1500 and 2000A 
three well-marked transitions were found except 
for diethyl ketone, for which the second (expected 
at about 1720A) is unaccountably missing. The 
longer wave-length transition of this ketone is 
also unaccompanied by vibrational structure, in 
contrast to all other ketones studied. 

Before discussing the spectra of these ketones, 
the spectrum of acetone must be mentioned, 
since the more numerous electronic energy levels 
are more clearly marked and serve as a basis for 
comparison with the higher ketones. Reconsider- 
ation of the results published in 1934” still does 
not explain all the details; to do this, probably 
additional experiments with the higher dis- 
persion available here are necessary and are 
planned. But some features not stressed in this 
earlier paper should be mentioned. (1) The most 
prominent feature of the spectrum as a whole is 
the Rydberg series reported leading to an 
ionization potential of 10.20 ev (using more 
recent values of the fundamental constants). If 
the Rydberg denominator is written (7+0.505) 
instead of (n—0.495) retaining v.=82,767 we get 
for n=1 the value 34,318 cm. The maximum 
in the near ultraviolet absorption is at 35,400 
(+ about 700) cm-'. The intensity for n=2 is 
about 10* times the intensity of »=1, but there 
may be an explanation for this not included in 
recent theory. (2) The second most prominent 
features in the far ultraviolet, appearing at 
pressures only slightly greater than this Rydberg 
series, are the band systems beginning at 51,171 
and 69,300. These have similar vibrational 
structure though that of the former is much more 
elaborate. Also the vibrational structure in the 
region 72,000—78,000 is complicated and probably 





* Boyce and Robinson, J. Opt. Soc. Am. 26, 133 (1936). 





TABLE I. Electronic term values (v..—v in cm), 











OXYGEN (ATOM) ACETONE FORMALDEHYDE 
36,069 3s(5S) (46,367) ~— 
33,043 35(8S) 31,601 ~ 
14,358 4s(5S) 17,554 23,440 
13,612 4s(%S) 13,456 16,230 
7721 + ~=5s(8S) 8933 10,040 
7425 5s(3S) 7396 8410 
4818 6s(5S) 5407 5960 
4671 6s(S) — 5110 








associated with more than one electronic transi- 
tion. The result is that a second series can be 
proposed (which was found in 1934 but rejected 
because of intensity anomalies): »=82,772—R/ 
(1+0.865)?, n=1, 2, 3, ---. Calculated values 
for n=1, 2, 3, 4 are 51,171, 69,353, 75,376 and 
78,085. Observed values are 51,171, 69,316, and 
75,376. Objections to the series are (a) failure to 
observe 78,085 while bands at 77,646 and 77,483 
are observed and (b) 69,316 is less intense than 
75,376 or 51,171. In favor of this second series 
is the obvious mathematical fit (perhaps an 
objection) and the fact that Price® observed two 
series in formaldehyde converging to the same 
ionization potential. 

If the two series are real (including »=1 of 
the first series), it is possible to explain by 
Rydberg transitions all the electronic transitions 
of acetone except the weak absorption between 
60,000 and 69,000 (two electronic transitions 
with complicated vibrational structure about 
10-* as strong as the other absorption below 
51,000) and a strong band appearing at low 
pressures at 72,529. If all the Rydberg transitions 
are of the type 2p —ms, involving practically 
atomic orbitals, we may expect a very rough 
agreement with term values in oxygen. Table I 
shows these values® and the results of Price in 
the same range of m values. These values are 
obtained by subtracting the origins of the 


TABLE II. Electronic transitions of the ketones (cm). 








I II III IV 





Acetone 36,400 | 51,171 | 60,086 | 65,213 
Methy] ethyl 36,100 | 50,634 | 58,110 | 63,525 
Methyl n-propyl 35,840 | 50,699 | 58,200 | 63,500 
Methyl isopropyl 35,460 | 50,460 | 57,566 | 63,030 
Diethyl 35,970 | 50,700 — 63,300 




















5 W. C. Price, J. Chem. Phys. 3, 256 (1935). 
6 Bacher and Goudsmit, Atomic Energy States (1932). 






























TABLE III. Vibrational transitions of the ketones. 








METHYL ETHYL METHYL ISOPROPYL METHYL 1-PROPYL 





50,634 50,460 50,699 
50,857 vot we 
51,209 vot ws 50,836 vot w4 51,313 vot ws 
51,915 w+ 51,791 w+; 51,889 w+ 
53,119 = v9+21 51,279 vo+we 
52,426 vot wits 52,099 vot2we 
wi = 1281 w,=1231 w= 1190 
wa= 575 w= 388 w= 614 
we= 220 we= 819 








electronic transitions from the spectroscopic 
ionization potential J for formaldehyde and 
acetone. The agreement in order of magnitude 
lends some support to the idea that of these 
transitions (except possibly 36,400 cm in 
acetone) are really Rydberg transitions of the 
above mentioned type. 

The observed transitions of the other ketones 
are shown in Table II, where corresponding 
values for acetone and methyl ethyl ketone are 
included. Transitions number I, in the near 
ultraviolet, were measured in the liquid state in 
most cases, where only maxima occur. These 
values are quite uncertain while measurements 
for transitions II, III and IV have at most an 
uncertainty of about 20 cm~. The shift toward 
longer wave-lengths of the absorption maximum 
which occurs when the hydrogens of acetone are 
substituted by simple radicals has been discussed 
for transitions I.” The shifts in the other transi- 
tions are to be associated with changes in the 
pure electronic terms as a consequence of change 
in structure. 

In all of the latter transitions, substitution of 
the first hydrogen by methyl (acetone to methyl 
ethyl ketone) produces the greatest shift. A 
further shift is produced only when an a hydro- 
gen is substituted to give an unsymmetrical 
structure (methyl ethyl to methyl isopropyl). 
Substitution to produce a symmetrical structure 
(diethyl) or 6 substitution indeed seems to show 
shifts mostly in the opposite direction, but these 
are relatively small. It is to be noted also that 
shifts are larger for III and IV than for II. 

A qualitative explanation of these shifts can 
be given® if it is admitted that all these transi- 
7K. L. Wolf, Zeits. f. physik. Chemie B2, 39 (1929). 

8 Compare the discussion of the alkyl halides by W. C. 


Price, J. Chem. Phys. 4, 547 (1936), where references to 
discussion by R. S. Mulliken will be found. 
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tions except III are 2p ,—ms, since because of 
their spherical distribution the s orbitals pene- 
trate the alkyl groups and are affected by changes 
in them. Shifts in IV would be greater than in 
II, since II are 3s while IV are 4s terms. While 
it is evident that shifts in III are even larger 
than in IV, any attempt to revise the present 
analysis in order to include transitions III in 
2po—ns would cause greater intensity difficulties 
and necessarily exclude a transition which was 
10* more intense. Possibilities of excitation of 
other electrons are obviously not excluded when 
the electronic structures of ketones suggested by 
Mulliken® are considered. While the ionizations 
potentials corresponding to the removal of 2), 
cannot be predicted quantitatively for these 
ketones, they should be only slightly less than 
10.2 ev, and it might be guessed that methyl 
isopropyl would have the lowest. 

Transition II is associated with vibrational 
transitions in the propyl ketones. No vibrational 
structure could be found for diethyl ketone, and 
no structure for the propyl ketones was found 
in III and IV. The long wave edges, which 
appear sharper, are given in Table III with 
probable assignments. The data for methyl 
ethyl ketone is not new, but is included for 
comparison. w; is found in all ketones studied 
and is probably the carbonyl frequency. We 
hope to decide if this is true in acetone’ from 
proposed measurements on acetone-dg. Little can 
be said about ws, w2 and ws except to point out 
that because of their magnitude they probably 
represent deformation of the carbon-oxygen 
framework. It is perhaps stretching a point to 
assign C2, symmetry to acetone though for 
normal vibrations involving little hydrogen 
motion it is a permissible approximation. How- 
ever it is doubtful if the higher ketones possess 
even C, (plane) symmetry but probably only Ci. 
Hence all vibrations (and overtones) should be 
permitted in the excited electronic states. There 
is little hope at present of calculating relative 
intensities, so we cannot say why so few 
vibrations actually appear or which are to be 
expected. But it is a striking fact that very few 


appear. 
9R.S. Mulliken, J. Chem. Phys. 3, 564 (1935). 
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Note on the Infra-Red Absorption Spectra of the Vapors of DCOOH and DCOOD 


R. C. HERMAN AND V. WILLIAMS 
Palmer Physical Laboratory, Princeton University, Princeton, New Jersey 


(Received February 28, 1940) 


N order to complete the set of infra-red spectra 

of the deuterated formic acids begun by 
Bonner and Hofstadter! (HCOOH) and Hof- 
stadter? (HCOOD), the dimer and monomer 
spectra of the vapors of DCOOH and DCOOD 
have been obtained in the region 1 to 15u. We 
hoped that sufficiently regular shifts in band 
positions would appear from one spectrum to 


1L. G. Bonner and R. Hofstadter, J. Chem. Phys. 6, 
531 (1938). 
2 R. Hofstadter, J. Chem. Phys. 6, 540 (1938). 
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= i) bad 
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the next so that approximate calculations would 
enable us to make assignments of bands whose 
frequencies lie below 1700 cm—. 

Heavy formic acid (DCOOD) was prepared 
by the thermal decomposition of anhydrous 
deuterated oxalic acid [(COOD)> ] in vacuum at 
about 180°C. The deuterated oxalic acid was 
obtained by exchanging anhydrous oxalic acid 
(Kahlbaum) with 99.6 percent DO and was 
dehydrated by heating under vacuum. The 


Fic. 1. Absorp- 
tion spectrum of 
DCOOH. Cell 
length 22 cm. 
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DCOOD was purified by several vacuum distil- 
lations after pumping off the COs. The DCOOH 
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Fic. 2. Absorp- 
tion spectrum of 
DCOOD. Cell 
length 22 cm. 
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several times with HO. Anhydrous copper 
sulfate was used to dehydrate the acid between 














was prepared by exchanging the DCOOD exchanges. 
TABLE I. 
DCOOD DCOOH 
DIMER MONOMER DIMER MONOMER a 
e v(ODO) 2680 cm7 v(OD) 3175 cm v(OHO) 3633 cm7 v(OH) 
2325 cm™ { »(CD) 2265 (CD) 2248 (CD) 2250 »(CD) 
1719 v(C=O) 1730 v(C =O) 1730 v(C=O) 1730 v(C=0) 
1243 1178 1358 1206 
1178 1143 1224 1146 
1150 1052 1150 1114 
1052 1015 1121 982 
974 957 990 946 
927 926 944* 678 
887 




















* Broad. 








INFRA-RED SPECTRA OF DCOOH AND DCOOD 
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Fic. 3. Line graph for the monomers of formic acid. 
Heights of lines indicate intensities approximately. 
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A preliminary experiment was performed on 
the thermal decomposition of light oxalic acid 
in vacuum at about 180°C. The resultant 
HCOOH was tested for water content by com- 
parison of vapor pressure measurements with 
those made by Coolidge* for formic acid-water 
mixtures. Further, the infra-red spectrum of the 
light acid showed no absorption band in the 
region of 1600 cm— where water absorbs strongly. 
These tests showed that there was less than 
3 percent of HO present. The high yield of 
formic acid obtained substantiates the work of 
Wobbe and Noyes? on the thermal decomposition 
of oxalic acid. Unfortunately, it was not possible 
to determine the extent of the deuteration of the 
formic acids by means of the mass spectrograph. 
However, the spectra themselves indicate that the 
CD, OD, and OH groups are at least 95 percent 
pure. The experimental technique was the same 
as that employed in previous investigations.° 

The monomer and dimer spectra of DCOOH 
and DCOOD are shown in Figs. 1 and 2. Parts 
of the monomer spectra of HCOOH and HCOOD 
were restudied and a few corrections and addi- 
tions are made. In the spectrum of HCOOH 
v(C=O) is found at 1760 cm-, the double band 


3A. S. Coolidge, J. Am. Chem. Soc. 50, 2166 (1928). 

4D. E. Wobbe and W. A. Noyes, Jr., J. Am. Chem. 
Soc. 48, 2856 (1926). 

5R. C. Herman and R. Hofstadter, J. Chem. Phys. 6, 
534 (1938). 
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listed at 1093 has minima at 1111 and 1075 cm~, 
and there is also a true monomer doublet with 
minima at 1221 and 1190 cm~. In the monomer 
spectrum of HCOOD the following minima were 
redetermined : 1178, 1142, and 1018 cm. Fig. 3 
shows the positions of all the bands below 
1800 cm in the spectra of the four isotopic 
monomer forms of formic acid. Table I gives 
the positions and some assignments for the bands 
in the monomer and dimer spectra of DCOOH 
and DCOOD. 

We are unable to give a consistent interpreta- 
tion of the bands whose frequencies lie below 
1700 cm. The apparently irregular shifts that 
occur in the spectra upon deuterium substitution 
seem to indicate either accidental degeneracies 
or strong interactions which are different for 
the different molecules. The simple triatomic 
calculation suggested by Bonner and Hofstadter 
is certainly inadequate to explain the spectra. 

It is interesting that practically all the bands 
in the monomer spectra below 1700 cm- show 
a doublet structure with a separation of about 
34 cm which is roughly the separation to be 
inferred from the spectrum of the 3v(OH) band 
in formic acid at 10,203A studied by Bauer and 
Badger. It seems difficult to understand why 
these bands seem to show a persistent doublet 
structure since any vibrational mode which has 
an appreciable change in dipole moment parallel 
to the least moment of inertia axis should 
exhibit a strong gQ branch. It would seem worth 
while to study the fundamental bands of the 
various formic acids under higher resolution. 

We wish to express our appreciation to the 
Research Department of the Calco Chemical 
Division of American Cyanamid Company for 
the interest shown in this work. 


6S. H. Bauer and R. M. Badger, J. Chem. Phys. 5, 852 
(1937). 
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Effect of Electric Fields on the Decomposition of Nitrous Oxide by Alpha-Rays 


A. D. KoLuMBAN AND H. Essex 
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Ion yields expressed in molecules of gaseous products uncondensed by liquid air, per ion 


pair are 
30° no field 


half-saturation 
full saturation 


100° no field 


half-saturation 
full saturation 


4.35 
4.15 
4.61 
4.75 
4.65 
5.76. 


If these results are interpreted on the assumption that ions are neutralized at the electrodes 
without decomposition, then it follows that in the absence of an electric field, 9 percent of the 
reaction, is initiated by combination of ions at 30°C, and 4 percent of the reaction is so initiated 
at 100°C. The increase in ion yield at high field strength is probably due to electron collisions 
with molecules. As this occurs at field strengths at which the current-field strength curve is 
horizontal, the new reaction is not initiated by ionization. Therefore, as the major part of the 
energy of alpha-rays remains in the primarily ejected electrons, electron collisions without 
ionization should initiate part of the reaction occurring in the absence of electric fields. 





N chemical reactions in the presence of 
ionizing radiations, it is possible to separate 
that portion of the reaction consequent on ion 
recombination from that portion of the reaction 
independent of ion recombination by comparing 
the rates of the reaction in the absence and 
presence of electric fields. In the presence of 
an electric field between metal plates, ions can 
be removed from the gaseous mixture before 
recombination takes place and neutralized upon 
the electrodes without decomposition. If the 
field strength is sufficient to produce a saturation 
current all of the ions are removed from the gas 
phase before recombination. Any reaction still 
taking place under such conditions must be 
initiated by other mechanisms than ion 
recombination. 

In this laboratory a method has been worked 
out for determining the ion yield in a reaction 
vessel in which the number of ions produced in 
unit time can be measured directly by the 
saturation current between metal electrodes and 
compared with the measured rate of the reaction. 
This method of determining ion yield is also 
directly applicable to the determination of the 
ion yield in electric fields and has been so used 
in the decomposition of ammonia by alpha-rays 
by Smith and Essex.! In that reaction 40 percent 


1 Smith and Essex, J. Chem. Phys. 6, 188 (1938). 


of the reaction was found to be initiated by ion 
recombination. The same technique has now been 
applied to the determination of the ion yields in 
the presence and absence of electric fields in the 
decomposition of nitrous oxide. 


EXPERIMENTAL 


The experimental technique of these measure- 
ments was essentially the same as that used in 
the decomposition of ammonia and described in 
considerable detail in the article cited. The 
nitrous oxide used was a _ cylinder-contained 
product of the Ohio Chemical Company, and 
stated by them to have a minimum nitrous oxide 
content of 99.5 percent. It was passed into the 
reaction system over soda lime and barium oxide 
to remove water and carbon dioxide and in the 
reaction system was separated from noncon- 
densible gases by evacuation after freezing out 
in a trap with liquid air.2 To insure complete 
removal of noncondensible gases, the cycle— 
condensation, evacuation, vaporization, conden- 
sation was repeated six times with each sample. 

The reaction system is shown in Fig. 1. The 
cylindrical reaction vessel was provided with 


2 The liquid air used in this investigation was provided 
by the General Electric Company through the courtesy 
of Dr. L. M. Willey, to whom the authors express their 
gratitude. 
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platinum plate electrodes placed very close (2 
mm) to the ends of the cylindrical vessel so that 
a minimum amount of ionization was produced 
back of the electrodes where some of the ions 
might escape inclusion in the saturation current. 
The radioactive material, mesothorium, was 
placed in the depression in the bottom of the 
tube. In one series of runs the reaction system 
was sealed after removal of the noncondensible 
gases and the reaction rate measured by the 
change in the residual pressure after freezing out 
the nitrous oxide by liquid air in trap I. The 
number of ions produced in unit time was 
determined by measuring the saturation current 
between the electrodes. In making the runs at no 
field the experiments were continued for a week 
or more. Saturation current measurements, and 
pressure measurements on the McLeod gauge 
after freezing out the nitrous oxide with liquid 
air were made initially after 24 hr., and after 
that, at intervals of several days each during the 
run. Similar technique was employed in making 
the measurements in a field: but for the measure- 
ments in a field except when determining the 
saturation current, or measuring the residual 
pressure, a constant potential was applied be- 
tween the electrodes during the entire run. In 
order to prevent arcing between the electrodes 
and the glass walls of the reaction vessel thin 
copper shields were sealed to the outside of the 
glass ends of the reaction vessel with nitro- 
cellulose cement in which fine silver filings were 
mixed. The silver filings effected electrical con- 
tact between glass and shield. The shield at the 
high potential end was electrically connected to 
the high potential lead. The shield at the low 
potential end was grounded and provided with 
a hole through which passed the lead from the 
electrode to the ion-current measuring galva- 
nometer. In this way, the current through or 
over the glass was not included in the ion- 
current. The reaction vessel was enclosed in a 
constant temperature bath and runs were made 
at 30°C and 100°C. Pressures were in the 
neighborhood of atmospheric. Preliminary meas- 
urements showed that the error introduced due 
to desorption of noncondensible gas from the 
walls was negligible. 

In a second series of runs the reaction system 
was connected through a mercury sealed glass 
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Cc 
6 Oe 
, + 
B 


Fic. 1. B—Constant temperature air bath, R—reaction 
vessel, G—McLeod gauge, E—electrodes, P—Q—thin 
copper shields to prevent arcing, 7—trap for condensing 
out nitrous oxide, C—capillary for sealing up system. 





























stopcock to a mercury vapor vacuum pump 
backed by an oil pump. In these runs after the 
residual pressure had reached a value too great 
to be measured by the McLeod gauge, the NO 
was frozen out in the trap by liquid air, the 
system evacuated by opening the stopcock, and 
after closing the stopcock a second run was 
made. No appreciable difference was observed in 
the experimental data obtained with the stopcock 
enclosed, and with the completely sealed appa- 
ratus: Therefore the data are not separated in 
the table of experimental results. With the 
stopcock enclosed system, pressures slightly 


‘greater than atmospheric, were used, to prevent 


leak in of gases, not condensible by liquid air. 
Slight outward leak of N.O is unimportant. 

The production of noncondensible gases under 
the conditions of these experiments, exhibited a 
considerable induction period lasting at least 24 
hours. The current potential curves obtained 
across the reaction vessel became horizontal 
above 16,000 volts at 30°C and above 9000 volts 
at 100°C indicating complete removal of ions, 
i.e., no combination of ions in the gaseous phase. 

The experimental results are presented in 
Table I. Ion yields are expressed in terms of 
moles of noncondensible gases produced per ion 
pair. Partly because of the small amount of 
decomposition occurring in such a comparatively 
brief time, the ion yields obtained during the 
first 24 hours show much greater divergence from 
each other than the ion yields obtained after 
24 hours. In Fig. 2 the average ion yields are 
plotted against the ratio of the plate current to 
the saturation current. In Fig. 3 the ion yields 
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and current are plotted against the field strength. 
Between the experimentally determined points, 
the curves are given their probable form. 


DISCUSSION OF RESULTS 


The decomposition of N.O by alpha-rays has 
been studied by Wourtzel.* The ion yield 
obtained as recalculated by Mund? is 2.7 at 
18°C and 4.0 at 315°C. 

G. R. Gedye® has investigated the N,O 
decomposition by cathode rays and obtains an 
ion yield of 3.9 at 18°C. 

There is general agreement among the investi- 
gators who have studied the decomposition of 
N.O in the presence of alpha-rays and also in 
the presence of cathode rays that at a particular 
stage of the reaction the products are nitrogen, 
nitric oxide and oxygen. The balanced equation 
for their production is 


4N,0-3N2+2NO+0,2 (1) 
(Oz produced partly as O). 


In the experiments of other investigators in 
which the partial pressures of NO and Os» were 
considerable, reaction may be expected to occur 
according to the equation 2NO+O2=2NOs. The 
kinetics of this reaction have been studied by 
Bodenstein.® The reaction is termolecular and 


3 Wourtzel, Le Radium 11, 289, 332 (1919). 

4 Mund, J. Phys. Chem. 30, 890 (1926). 

5 Gedye, J. Chem. Soc. 2, 3016 (1931). 

6 Bodenstein, Zeits. f. physik. Chemie 100, 68 (1922). 


TABLE I. Jon yields at 30°C. 








Ion YIELD 
No FIELD 


Dur- 
LENGTH ING FRoM 24 LENGTH 
EXPpERI- | OF RUN First | Hours | Experi-| or RUN 
MENT IN 24 TO END MENT IN 
NuMBER | Hours Hovurs| oF RuN | NUMBER| Hours 


Ion YIELD FULL 
SATURATION 


Ion YIELD HALF- 
SATURATION 


Dur- . Dur- 
ING FROM 24 LENGTH ING From 24 
First | Hours | Experi-| oF RuN | PREs- | First | Hours 
24 TO END MENT IN SURE 24 TO END 
Hours| oF RUN | NUMBER| Hours | IN MM | Hours| OF RUN 





11 313 3.82 | 4.02 | 13A 127.5 
11B 236 3.38 | 4.42 | 13B-1 | 240.5 
11C 214 3.06 | 4.17 | 13B-2 | 182 
11D 217.5 4.53 | 4.45 
11#-1 | 148.5 4.17 | 4.66 
11H-2 | 244 3.85 | 4.26 
11Z-3 | 241 4.45 


3.10 | 4.14 12A 109.5 | 487 | 4.31 | 4.32 
3.63 | 4.24 12B | 240 764 | 3.91 | 4.90 
3.69 | 4.08 





Average 3.75 | 4.35 

















14 
14A 
14B 
14BB 














Average 












































—- (FP — A) 


s has 
yield 
7 at 


N:O 


IS an 


vesti- 
on of 
so in 
cular 
ogen, 
ation 


(1) 


rs in 
were 
oecur 
. The 
d by 
and 


DECOMPOSITION OF NO BY"”ALPHA-RAYS 


the velocity constant is so small at the tempera- 
tures of these measurements, that appreciable 
combination is not possible at the very low 
partial pressures of NO and Oz, reached in these 
experiments. On the other hand, part of the 
oxygen must originally be produced as atoms 
and some of these atoms of oxygen may react 
with NO to give NO». Other oxygen atoms will 
react with each other to give Os. From the NO» 
produced in this way the following subsequent 
reactions may be expected to occur 


2NO2+2Hg = 2HgNO,(s), (2) 
2HgNO2(s) = 2HgO+2NO (3) 


as in the photochemical decomposition of N,O.? 

HgNO2sz is a solid and its decomposition is very 
slow accounting for the induction period ob- 
served in our experiments. From the initial and 
final ion yields (in moles of noncondensible gases 
formed per ion pair) it is possible to estimate 
the ion yield in moles of NxO decomposed per 
ion pair for comparison of these results with 
those of Mund.‘ It is apparent from the Eqs. 
(1), (2) and (3) that if X is the number of 
moles of NO which react with O to give NOs, 
the number of moles of noncondensible gas 
formed from 4 moles of N,O in the initial stages 
of the reaction is 


x 3x 
34+2-—x+1--=6-—— 
2 2 


and the number of moles of noncondensible gas 
formed from 4 moles of N,O after the stationary 
state or “‘radioactive’’ equilibrium is established 
is 


zs sf 
3+2+1--=6--. 
2 2 


* may therefore be obtained from the experi- 
mental data from the relation 





6—3x/2 ion yield in the first 24 hours 


6— 1x/2- ion yield after 24 hours 


At 30° and 100°C these values are 
30°C 100°C 


lon yield first 24 hr. 3.75 4.14 
Ion yield after 24 hr. 4.35 4.75 


"W. A. Noyes, Jr., J. Chem. Phys. 5, 807 (1937). 
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At 30°C x=0.776 and at 100°C x=0.723 from 


which the ion yield in terms of moles of N2O 
decomposed per ion pair is, 


at 30°C 


4 
—— X4.35 = 3.10 
6—0.388 


+ 
at 100°C ———— X4.75 =3.37. 
362 


The value 3.10 obtained at 30° may be compared 
with 2.75, the value at 30° from Mund’s calcu- 
lations from Wourtzel’s experimental data. As 
the intensity of ionization in our experiments was 
a great deal less than that of Wourtzel, the 
difference between our ion yield and that of 
Wourtzel is in agreement with the frequently 
observed fact that the ion yield in alpha-ray 
decompositions increases with a decrease in 
intensity of ionization. In our experiments on 
the decomposition of ammonia by alpha-rays at 
low intensities of ionization, the ion yield agreed 
excellently with the extrapolated ion yield for 
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this reaction at zero intensity of ionization as 
calculated from the measurements of Luyckx 
on the change in the ion yield with intensity of 
ionization. No data is available in the literature 
on the decomposition of N.O at low intensities 
of ionization, with which the results of this 
work might be directly compared. 

The temperature coefficient of the ion yield is, 
according to Mund 0.0044 per degree and from 
these measurements 0.0039 per degree. 

The ion yield consequent on recombination of 
ions is calculated from the ion yield in no field 
and the ion yield in a field of sufficient strength 
to produce a plate current equal to half the 
saturation current by doubling the difference 
between the ion yield at zero field strength and 
at half saturation. The results of such a calcula- 
tion show that at 30° 9.2 percent of the reaction 
is initiated by combination of ions and 90.8 
percent due to other mechanisms and at 100°C 
4.2 percent is due to recombination of ions and 
95.8 percent to other mechanisms. The ion yields 
at voltages sufficient to produce saturation can 
not be used for such calculations because at high 
voltages decomposition is produced by a new 
mechanism, probably by collision of electrons 
with molecules. The amount of decomposition 
due to this new reaction which is apparent at 
high field strength can be safely neglected at 
voltages sufficient to produce half-saturation as 
these voltages are very low compared with the 
voltages necessary for saturation, i.e., 900 volts 
compared with 17,000 volts at 30°C and 500 
volts compared with 9000 volts at 100°C. 

In Fig. 3 are plotted the curves of ion yield 
versus field strength at 30° and 100°C together 
with the current-field strength curves at the 
same temperatures. Comparison of the curves 
shows that decomposition due to high field 
strength is apparent far below the field strength 
necessary for ionization by electron collision. 
Similar relations are apparent in the work on 
the alpha-ray decomposition of ammonia.! 

These facts point definitely to the decomposi- 
tion of N,O by electrons having far less than 
sufficient energy for ionization and therefore to 
this mechanism in the decomposition of N,O by 


alpha-rays in the absence of electric fields, due 
to ejected electrons. 

Ion yields of the observed magnitudes are 
energetically possible even though no reaction 
accompanies the recombination of ions. In air 
the energy of an alpha-particle consumed per 
ion pair produced is 35 electron volts.* The 
relative ionization in N2O compared with air is 
1.05 ;® from which figures the alpha-ray energy 
consumed per ion pair in N2O is 35X1/1.05 
= 33.3 ev. Most of the energy lost in the passage 
of an alpha-particle must be lost in inelastic 
collisions by the primarily produced electrons. 
Quantitative data for the fraction of energy lost 
by inelastic collisions in N2O are not available 
but for the passage of low energy electrons 
through hydrogen, according to Lunt and Meek” 
92 percent of the energy is so lost. Assuming 
the same efficiency in N,O, the energy so lost 
per ion pair is 33.3X0.92=30.6 ev. The most 
abundant positive ions produced in N;O by elec- 
tron collisions are N.O+ for which the ionization 
potential is 12.9.1. The energy lost per ion pair 
by other inelastic collisions is therefore 30.6 
—12.9=17.7 ev, sufficient to split 10 molecules 
of N2O by the mechanism 


N2O=N:2+0O; AE=1.7 ev. 


Part of this reaction probably occurs with elec- 
tron attachment to oxygen 


N,O+e= N2+0-. 


Electrons having energies in excess of 1.7 volts 
are capable of attachment to N,O in this manner 
according to Bradbury and Tatel.!? Formation 
of negative ions by this mechanism also probably 
accounts for the increased ion yields at high field 
strengths. 

Studies of the alpha-ray decomposition of 
nitrous oxide and ammonia at higher field 
strengths are now being made in this laboratory. 


8 Rutherford, Chadwick and Elliss, Radiations from 
Radioactive Substances, p. 80-81. 

9R. D. Kleeman, Proc. Roy. Soc. A79, 220 (1907). 

10 Lunt and Meek, Proc. Roy. Soc. A157, 146 (1936). 

1 E. C. G. Stueckelberg and H. D. Smyth, Phys. Rev. 
36, 478 (1930). 

12 Bradbury and Tatel, J. Chem. Phys. 2, 835 (1934). 
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Note on the Ultraviolet Absorption Systems of Benzene Vapor 


G. NORDHEIM AND H, Sponer, Department of Physics, Duke University, Durham, North Carolina 


AND 


E. TELLER, Department of Physics, The George Washington University, Washington, D. C. 
(Received March 21, 1940) 


An attempt has been made to give a coherent interpretation of the absorption systems of 
benzene vapor in the far ultraviolet including the Rydberg series. The following assignments are 
proposed: A forbidden transition of symmetry 1A 1,—>'B,, is suggested for the bands at 2050 
to 1850A and an allowed !A,,—'E~, transition for the much more intense bands at 1850 to 
1650A. The continuous background in the latter region may possibly be due to a transition 
involving a C—H dissociation. The observed two Rydberg series can both be assigned to 


allowed transitions of symmetry 'A1,—'E~y. 





1. INTERPRETATION OF THE SYSTEMS AT 2000A 
AND 1750A 


HE first absorption region of benzene lies at 
2700—2200A. It represents an electronic 
transition 1A ,,—'B2, which is forbidden! by sym- 
metry and is made allowed when vibrations of 
type” e+, are excited. The analysis* of the vibra- 
tional structure of this absorption system has 
confirmed these assumptions. The interacting et, 
vibration is the Raman active vibration of 
frequency 606 cm. 

A stronger absorption spectrum has been ob- 
served‘ in the region from 2050A to 1850A. It 
consists of a number of diffuse broad groups 
showing some structure. The spectrum retained 
its diffuseness even when the absorption tube 
was surrounded by dry ice.’ Below 1850A the 
bands are overlapped by a still stronger absorp- 
tion system which has the appearance of a 
continuous background upon which bands are 
superimposed, which are much sharper than the 
bands of the system just mentioned. 

The strong intensity of the absorption below 
1850A suggests an allowed transition. According 
to calculations® of electronic states an allowed 
~ TAL, Sklar, J. Chem. Phys. 5, 669 (1937). 

*Symmetry symbols for vibrations are denoted with 
small Greek letters. 

*H. Sponer, G. Nordheim, A. L. Sklar and E. Teller, 
J. Chem. Phys. 7, 207 (1939). 

(1938) P. Carr and H. Stiicklen, J. Chem. Phys. 6, 55 

’ Personal communication from Dr. H. Stiicklen for 
which we are much indebted. 

°E. Hiickel, Zeits. f. Physik 70, 204 (1931); Zeits. f. 
Elektrochem. 43, 752 (1937); A. L. Sklar, reference 1; 


M. Goeppert-Mayer and A. L. Sklar, J. Chem. Phys. 6, 
645 (1938). " J < . 


, 


transition is to be expected in this region. The 
transition can be described by the excitation of a 
double bond electron of character e~, into an et, 
level.” The resulting electronic configuration can 
be represented by a hole (e-,)— left in the A, 
ground level multiplied by the wave function of 
an e+, electron. The product (e-,)—e+, yields 
molecular wave functions of symmetry B,,, Bou, 
E-,. Transitions to the first two levels are for- 
bidden while A1,—(e~,)—“'e+.E~, is allowed. (All 
discussed transitions concern singlet terms.) We 
tentatively interpret the relatively sharp bands 
below 1850A as this transition. It has been 
shown’ that in a degenerate electronic state the 
symmetrical nuclear configuration does not corre- 
spond to a stable equilibrium. This holds for the 
upper E~-, state whose degeneracy can be re- 
moved by a nuclear displacement of e+, sym- 
metry. Hence the equilibrium configuration 
differs from the symmetrical one by an et, dis- 
placement. The deviation from the symmetrical 
configuration will, however, be small because the 
single bond electrons which are predominant in 
determining the shape of the molecule, are left 
undisturbed in this transition and only one 
double bond electron is excited. The slight 
change in equilibrium configuration may account 
for the somewhat irregular vibrational structure 
of the A,,—E~, transition. We are inclined to 
believe that the apparently continuous absorp- 
tion which forms the background corresponds to 


7If the symmetry symbols refer to single electrons, 
capital letters are eee by small ones. 

8 A. H. Jahn and E. Teller, Proc. Roy. Soc. Lond, A161, 
220 (1937). 
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another electronic transition; namely, to a re- 
pulsive state leading to the rupture of one C—H 
bond.® 

The absorption region between 2050A and 
1850A represents the forbidden transition A,,— 
(e~,)—e*+,Biu. The transition A\,—(e~,)—~e+,Bou 
has been ascribed to the longest wave-length ab- 
sorption as mentioned in the beginning of this 
note. Both transitions are made allowed by 
vibrations of type e+, and the A,,—B,, transition 
also by a Bx, vibration.’ The assignments of the 
observed absorption systems to the two transi- 
tions are based upon calculations" of the relative 
positions of the B;, and Bz, levels. The observed 
occurrence of 965 and 160 progressions in the 
system A,,—B,, may be explained in an analo- 
gous way as for the near ultraviolet bands. It is 
interesting to note that the 965 progression 
which, as Carr and Stiicklen (cf. 4) have done, is 
best interpreted by an excitation of the totally 
symmetrical C—C vibration in the upper state, 
shows unusually strong convergence so that large 
anharmonicity of this vibration must be assumed. 
As the transition is forbidden, both a 1-0 and 
0—1 transition should be observed for the vibra- 
tion which makes the transition allowed. Ap- 
parently only the latter is present and will be one 
band of the first group. The former will be weaker 
on account of the Boltzmann factor and the 
corresponding band, being diffuse, may escape 
observation. In addition, it may be obscured by 
the tail of the predissociation observed in the 
A,,—B2, system. On the other hand, the neigh- 
borhood of the allowed E-, level, and the C—H 
repulsion curve may very well explain the ap- 
parent irregularity of the intensity of the 160 
progressions; bands in which excitation of a Be, 


®It cannot be excluded, however, that the strong 
apparently continuous region in reality represents the 
Aig(e~,) “e+. E~. transition and that the sharp bands 
are somewhat accidentally superimposed on it belonging 
to another electronic state which could represent the first 
member of one of the Rydberg series. This interpretation 
is favored by Dr. R. S. Mulliken (private communication, 
see also J. Chem. Phys. 7, 20 (1939)). 
_ 7A Bs, vibration causes mixing of the B;, state with 
an Ag state, and an Aj,—>Ao2, transition is allowed. 
However, there is no information about the energies of 
Ax, states in benzene, and therefore, it is not possible to 
estimate the intensity with which 82, vibrations could 
appear in the Ai,—B}, system. 

11M. Goeppert-Mayer and A. L. Sklar, J. Chem. Phys. 
6, 645 (1938). These authors mention already the possibility 
of assigning the B,, and E~, levels to the observed ab- 
sorptions in the way proposed here. 
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vibration rather than an e+, vibration makes the 
transition allowed may contribute to the irregular 
appearance of the system. A more detailed 
analysis of the spectrum seems impossible at the 
present moment owing to the diffuseness of the 
whole system. 

That the spectrum at 2050A is much stronger 
than the 2600A bands can be explained by the 
fact that the former lies closer to the allowed 
absorption system at 1770A. The intensity of a 
forbidden transition is proportional” to (1/AE)? 
where AE is the energy distance to the next 
allowed level of appropriate symmetry. Calcu- 
lating the AE values from the long wave ends of 
the systems A,,—E~u, Aty—Biu, A1g—Bou, one 
finds that the A,,—B,, transition should be 
about ten times stronger than the A,,—B2, 
system. The actual ratio can easily be larger 
because we have used in fact AE values for the 
equilibrium configuration of the nuclei in each 
state while the hypersurfaces of potential energy 
may approach closer for vibrating nuclei. This 
will be of greater relative importance for the 
higher lying B,, surface. A further contribution 
may arise if the upper state of the continuum at 
1770A is used instead of the E~, in calculating 
the AE values. If this continuum extends farther 
towards the red, then (1/A£)? will again show a 
greater increase for B,, than for Bo. The 
observed intensity ratio is indeed larger" 
than 1 : 10. 

The 2050A system may alternatively be as- 
cribed to a forbidden A1,—£*, transition, since 
an Et, level has been predicted by Sklar (cf. 1) 
in this energy region on the basis of the Heitler- 
London approximation method. The above con- 
siderations may as well be applied to an A1,—E*, 
transition, but the vibration which makes the 
transition allowed will be «~, rather than €*>. 
Moreover, the degeneracy of the E+, state may 
cause a slight asymmetry of the equilibrium 
configuration and hence a more involved vibra- 
tional structure. 


2G. Herzberg and E. Teller, Zeits. f. physik. Chemie 
B21, 410 (1933). See also considerations on intensities by 
R. S. Mulliken, J. Chem. Phys. 7, 20, 353 (1939). 

13 According to private communication from Dr. Carol 
A. Rieke, for which we are greatly indebted, the intensity 
ratio of about 1 : 1000 calculated by N. Q. Chako, J. 
Chem. Phys. 2,'644 (1934) is much too high. 
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—?T C-H DISTANCE 


Fic. 1. Schematic potential curves illustrating predissoci- 
ation and diffuseness in benzene spectra. 


2. PREDISSOCIATION AND DIFFUSENESS OF 
THE BANDS 


Predissociation and diffuseness of some of the 
bands may be explained by the rupture of a 
C—H bond. In Fig. 1 potential energies of the 
different levels are drawn as functions of one of 
the C—H distances denoted by r. The symmetry 
symbols attached are significant only in the 
neighborhood of the minima where the C—H 
distance in question is equal to the other C—H 
distances. The repulsive curve which we assume 
to represent the upper state of the continuous 
absorption at 1770A must have the symmetry 
E-, or Ao in the vicinity of the equilibrium 
distance in order to give rise to a strong ab- 
sorption. If we assume that the excitation con- 
sists in transferring a C—H bonding electron into 
an antibonding state, then the transition moment 
will lie in the molecular plane and the repulsive 
state must be E-,. Such a state can be con- 
structed by superimposing excitation in the 
different C—H bonds with appropriate phases. If 
now one C—H distance is increased, the super- 
imposed amplitudes change until, for great r 
values, the excitation is localized in that C—H 
bond which is to be broken. 

We shall expect predissociation in the transi- 
tion Ay,—Bz,, when the vibrational energy of 
B»,, reaches the intersection with the repulsive 
curve. Though we know from the spectrum that 
the vibrational energy is mostly accumulated in 
the totally symmetrical C-vibration, one may 
expect that on account of anharmonic forces the 
energy of oscillation may go over into the 


stretching of a C—H bond. The diffuseness of 
the A1,—B,, system is readily explained from 
Fig. 1 where the intersection between the re- 
pulsive curve and the B,, curve takes place at 
much lower vibrational energies. On the other 
hand, we explain the sharpness of the A,,—-E-, 
system by placing the E~-, curve just above the 
repulsive curve, hereby avoiding intersection. 
Fig. 1 gives the impression that the maximum of 
continuous absorption lies at slightly longer 
waves than the strongest bands of the A:1,—E~, 
system. Actually the maxima of the continuous 
and discrete systems seem to coincide. But our 
figure, taking into account only one degree of 
freedom, is an extremely poor representation of 
the polydimensional energy surfaces, and it is 
difficult to predict the exact location of such 
maxima within a few hundred wave numbers. 


3. INTERPRETATION OF THE RYDBERG SERIES 


In the far ultraviolet, very strong sharp ab- 
sorption bands have been observed and arranged 
into two Rydberg series.'* They begin at about 


1550A and extend till 1360A. The Rydberg 
formulae which have been proposed are as 
follows: 


y=74,495 —R/(n+0.97)? (1) 
y=74,590—R/(n+0.55)? n=3,4,5+++. (ID) 


In these formulae the two ionization limits 
almost coincide and correspond to a mean value 
of 9.19 volts. The most recent determination of 
the ionization potential by the electronic impact 
method gives 9.8+0.1 volts, a value which 
represents an upper limit according to the nature 
of these experiments. 

The intensity of both the Rydberg series 
definitely points to the fact that they are due to 
allowed transitions (i.e., the upper state must 
have the symmetry E~-, or Ao). The members of 
one series will belong to electronic states of the 
same symmetry. If we extrapolate series I toward 
the red, we come into the region of 1600A (n=2) 
and 2100A (m=1). In both these regions we 
observe only diffuse and structureless absorption 


14 7. C. Price and R. W. Wood, J. Chem. Phys. 3, 439 
1935). 

16 A, Hustrulid, P. Kusch and John T. Tate, Phys. Rev. 
54, 1037 (1938). 
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which corresponds to the tail ends of the above- 
mentioned 1850A and 2600A systems, respect- 
ively. Therefore, we are led to the conclusion that 
the system at about 1480A is the first observed 
member of series (I). 

On the other hand, the extrapolation of series 
(II) leads into the region of the 1770A system 
observed by Carr and Stiicklen and discussed 
above. Since we have earlier identified this 
strongest of all absorption systems with an 
allowed transition A,,—E~,, it seems suggestive 
to us to assign tentatively series (II) to transi- 
tions A 1,—(e~,)~!(ne+,,) E~, with the 1850-1650A 
system as its first member. That the intensity of 
the first member’® is much larger than that of the 
succeeding ones may be explained by the fact 
that it takes place much more within the 
molecular dimensions than the higher members.” 

Since the other Rydberg series (I) leads to the 
same ionization state, an (e-,) hole has to be 
expected in the upper states of series (I) as well 
as of series (II). The intensity of series (1) points 
also to allowed transitions. The only transi- 
tions satisfying both these requirements are 
A 1y—(€~g)“"(Mdou)E~y and A 1g—>(€~,)"(ne~4) A ou. 
The first of these two assignments seems to us 
more likely. In fact, in A1,—E~, the degenerate 
upper states would have a somewhat unsym- 
metrical nuclear equilibrium configuration, and 
thus a relatively complicated vibrational struc- 
ture should be expected. In the A1,—A 2, transi- 
tion, on the other hand, the occurrence of only 
totally symmetrical vibrations would be insuff- 


16See photograph p. 209 in W. C. Price and D. W, 
Walsh, Proc. Roy. Soc. Lond. A174, 207 (1940). 

17 The same intensity relations hold for atomic Rydberg 
series where also the first member is more influenced by 
the extended structure of the atom. 
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cient to account for the observations of Price, 
who finds a 960 and 670 energy difference associ- 
ated with series (I) in addition to smaller energy 
differences. If these energy differences correspond 
to a vibrational structure, they would be ex- 
plained more readily by an A,,—E~, transition. 

But the A1,—(e~,)~!(me~) Ao assignment can- 
not be excluded definitely. If the (e-,) hole and 
the (ne~,) electron would be coupled very weakly, 
the (e-,) hole would determine the equilibrium 
configuration of the nuclei. The degenerate e~, 
orbit would then lead to an unsymmetrical 
configuration and a complicated vibrational struc- 
ture, so that the results would be essentially the 
same as for A;,—E~,. The coupling described in 
this section obtains, however, only if the separa- 
tion of the various levels of the (e-,)—'(ne-,) 
configuration is less than the splitting of the 
(e-,)—! wave function by the zero point vibrations 
of the molecule. Analogy with atomic spectra 
leads us to an estimate of about 0.1 volt for the 
former ‘quantity in the first member of series (I). 
On the other hand, the splitting by the zero 
point vibration will probably be rather smaller 
than the total width of the vibrational structure 
which is about 1000 cm—!. Therefore, it is unlikely 
that A1,—A», could give rise to a sufficiently 
complicated vibrational pattern. 

There is, of course, also the possibility that 
some of the bands forming the pattern and 
especially some of the violet companions could be 
arranged into independent Rydberg series and in 
this case both transitions A1,—(e~,)—!(d2u) Eu 
and Ai,—(e~,)~'(ne~u)A2u. could be present. 
Theoretically this would be the most satisfactory 
alternative, since it would account for all possible 
allowed transitions in which the most loosely 
bound e-, electron is excited. 
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Spectroscopic Measurements of Gaseous CN 


II. Thermal Dissociation of Cyanogen* 


Joun U. Wuitet 
Department of Physics, University of California, Berkeley, California 


(Received January 22, 1940) 


Free CN radicals are formed from cyanogen at high temperature. Their absorption spectrum 
was observed in a furnace at 1500°K. From the lower limit of the absolute absorption coefficient 
found in Part I of this paper, their partial pressure and hence the equilibrium constant of the 
dissociation of cyanogen to free radicals are calculated. A lower limit for the heat of the re- 
action 

C:N:>2CN Kis00°K = 1.1 10- Ho=138.1(+1) kcal. 


follows directly. f=0.1 is estimated to be the most probable f value for the *2->*2 transition 
of CN. The most probable value of the heat of dissociation of cyanogen is calculated from this 
to be Ho=146(+4) kcal. By calculating the partial pressure of CN at about 400°C, the thermal 
polymerization of cyanogen is explained as the addition of the free radicals formed to cyanogen. 
molecules. From thermodynamic cycles, it is shown that the most probable value of the heat of 
dissociation of CO is 11.054 volts, and of the heat of vaporization of graphite 169.2 kcal., 
although 9.797 volts and 140.2 kcal. respectively, for these values cannot be completely 
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eliminated. 





INTRODUCTION 


HE absolute absorption coefficient of CN, 

which was found in the preceding paper,! 
gives a quantitative measure for the amount of 
CN radical present wherever its absorption is 
observed. In this paper the measure is con- 
sidered as a lower limit and applied to the CN 
formed by the thermal dissociation of cyanogen 
at high temperatures. When the partial pressures 
of CN and cyanogen are known, the equilibrium 
constant and free energy of the dissociation can 
be calculated directly. Lower limits and best 
values of the thermodynamic quantities are 
calculated from the absorption coefficient found 
previously and from a somewhat larger one 
estimated to be closer to’ the correct value, 
respectively. 


APPARATUS 


The same apparatus was used in this experi- 
ment as in the preceding one, with the substitu- 
tion of a high temperature furnace for the 
apparatus producing electrical dissociation. The 





*This and the preceding paper are condensed from a 
thesis submitted in partial fulfillment of the requirements 
for the degree of Doctor of Philosophy in the Department 
of Physics. 

t Now with the Standard Oil Development Company, 
Elizabeth, New Jersey. 

J. U. White, J. Chem. Phys. 8, 79 (1940). 


21-foot grating with the cylindrical lens in the 
second order was unchanged. Some of the same 
cyanogen was used. As there was no question 
of timing, sparks were run in the source as fast 
as power could be supplied to the condenser. 
Half a minute, at five sparks a second, was 
needed to give a good exposure through the 
absorption cell. The photographic plates and 
methods of calibration and measurement were 
identical. The furnace was made from an 
Alundum tube 70 cm long. The middle 60 cm 
were wound with 88 turns of carefully spaced 
Nichrome ribbon, 0.04 ohm per foot, set in 
Alundum cement. The 5 cm at each end were 
wound with number 18 wire and connected in a 
separate circuit. Quartz windows cemented on 


TABLE I. Temperature distribution in furnace. 











DISTANCE 
FROM CENTER ABSOLUTE 
OF FURNACE TEMPERATURE 
3.5 cm 1516°K Position of thermocouple 
during exposure 
5 1515 
8 1514 
10 1514 
13.5 1518 
17.5 1520 
19.5 1518 
21.5 1512 
23.5 1503 
25 1493 End of absorption cell 
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Fic. 1. Microphotometer trace of CN absorption at 1500°K. 


the ends of the furnace, in addition to those of 
the absorption tube, prevented convection, and 
asbestos cloth served as thermal insulation. The 
temperature distribution was adjusted by pulling 
asbestos off the middle of the furnace and 
draping it over the ends until the set of readings 
shown in Table I was obtained. Quartz windows 
were fused on a quartz tube 50 cm long and 
2.5 cm in diameter to make the absorption cell. 
When it was placed in the middle of the furnace, 
a correction of —1° was needed to reduce the 
temperature at the position of the thermocouple 
to the mean for the absorption cell. 

The thermocouple was made of platinum and 
platinum-rhodium wires, and readings were 
taken with a sensitive galvanometer. The melting 
and freezing points of samples of pure silver and 
gold were measured under the same conditions 
that existed later in the other temperature 
determinations. The observed differences be- 
tween the melting and freezing points were 3° 
and 2° for silver and gold, respectively. The 
averages were taken as the best values in drawing 
the calibration curve for the galvanometer. 


MEASUREMENTS 


After the furnace had been heated to 1100°C 
and allowed to come to equilibrium, the absorp- 
tion cell was pumped out and cyanogen admitted 
to a pressure of about 60 cm. An exposure was 
taken ten seconds later. The gas was then 
pumped out, some fresh admitted to } the 
pressure, and another exposure taken. This was 
repeated until the absorption disappeared. The 
furnace was then heated to a higher temperature 
and another set of exposures was taken. This 


was repeated at intervals of about 30° to the 
highest temperature, 1515°K. At this point, four 
exposures at pressures ranging from 64 cm down 
to 1 cm showed measurable absorption. Three 
exposures were also taken at about the same 
pressure at intervals ranging from 10 seconds to 
5 minutes after admitting the cyanogen to the 
cell. The final results show that there was no 
delay in establishing equilibrium. Only the 
exposures taken at the highest temperature 
were worked up. However, a visual comparison 
with the exposures at other temperatures showed 
that they would give about the same results if 
they were used. A microphotometer trace of one 
of the exposures, showing the high temperature 
distribution of intensity in the band and the 
formation of a head, is shown in Fig. 1. 


CALCULATIONS 


The method of measuring the plates and 
reducing the measurements to absorption coefh- 
cients was exactly the same as that described 
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Fic. 2. Relation between area and maximum 
intensity of the lines. 
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SPECTROSCOPIC MEASUREMENTS OF GASEOUS CN 


TABLE II. Equilibrium constants and partial pressures at 1515°K. 








TIME AFTER 


ADMITTING PRESSURE 





CYANOGEN CYANOGEN PP.CN K 

10 seconds 669 mm 680( +200) x 10-* mm 9.40(+7) X 107% 
1 minute 616 670 9.8 
5 minutes 655 710 10.1 

10 seconds 176 410 12.9 

10 seconds 44.5 200 12.3 

10 seconds 12.0 100 11.4 








before.t The areas of 20 different lines were 
measured on large scale microphotometer traces 
to obtain the new relation between area and 
maximum absorption. It is shown in Fig. 2. 
A new set of curves relating the amount of light 
absorbed by each line to the absorption coeffi- 
cient of its stronger component was calculated 
for 1500°K. It is shown in Fig. 3 with scales 
both of A/Jos and of the observed minimum of 
the fraction transmitted. Values of a were 
determined from it for each line and then }> xB 
was calculated for each exposure. If we recall 
that it was shown in Part I of this paper that 
an upper limit for the partial pressure of CN at 
300°K is 
0.0035 


Pox = —— B, 
x K 


and multiply by the ratio of the two tempera- 
tures, then 


0.0147 
Peonisisex ahaa 
x K 
RESULTS 


The partial pressures of CN calculated for the 
different observed intensities of absorption are 
tabulated in Table II. The pressure of cyanogen 
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Fic. 3. Relation between absorption coefficient and 
amount of light absorbed. 





and the corresponding equilibrium constants are 
also shown. The probable errors indicate the 
absolute uncertainty including that in the f value 
used. The consistency of the readings shows that 
their relative errors compared to each other are 
much smaller. The extent of the agreement of 
the first three values in the table shows the 
fluctuations to be expected between the different 
readings and indicates that there is no trouble 
from delay in reaching equilibrium. However, 
the variation with pressure is larger than can be 
accounted for by photometric errors. I have no 
explanation for this, but point out that in any 
case its effect on the results is small. The value 
of K=1.1(+0.7)X10-" in atmospheres was 
taken as the best value. At 1515°K a lower 
limit for the free energy of the dissociation 
C.N.2-2CN, is 


AF=-—RT |In K=82.9 kcal. 


The entropies of C2.N2 and CN were next 
calculated in the usual way.”* The vibration 
frequencies and moment of inertia of CN are 
well known, those of C2N»2 were taken from the 
paper of Stevenson‘ with the exception of 4e. 
His value of 230 cm~! was raised to 240 cm—, 
on the suggestion of Professor W. F. Giauque, 
to make the calculated entropy agree better 
with the thermal value that he and Ruhrwein 
obtained.» When these frequencies are substi- 
tuted in the table given by Gordon and Barnes 
and the proper calculations are made, the molal 
entropy of CN is found to be 60.6 E.U. at 
1500°K. That of C.Ne2 is 84.9 E.U. As the 


nuclear spin does not change in a chemical 


2 W. F. Giauque, J. Am. Chem. Soc. 52, 4808 (1930). 

9a) R. Gordon and C. Barnes, J. Chem. Phys. 1, 297 
(1933). 

4D. P. Stevenson, J. Chem. Phys. 7, 171 (1939). 

5R. A. Ruhrwein and W. F. Giauque, J. Am. Chem. 
Soc. 61, 2940 (1939). 
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TABLE III. Appearance potentials, excitation potentials, and 
relative abundances of particles present in the discharge. 








APPEARANCE RELATIVE EXCITATION 





PARTICLE POTENTIAL ABUNDANCE POTENTIAL 
C2oNe2* 14.1 volts 3.4 
C2N* 21.3 0.7 
Cot 18.6 0.2 
CNt 19.8 0.1 
CN 6.3 3.2, 1.3 volts 
At 15.5 
A 3.5 
C+N from CN 6.4 








reaction, it has been neglected in these calcula- 
tions. The entropy change in the dissociation is 
AS=36.3 E.U., from which we can reduce the 
free energy to 1500°. 


AF 4500° = 82.3(+ 1.5) kcal. 


The heat of dissociation of C2Nz2 into CN is, 
therefore, 
Ai500° = AF+ TAS 
= 136.8(+1.5) kcal. 


This is reduced to AH at absolute zero by adding 
the difference in the heat capacities of the 
reactants and products. If we assume rotation 
to be classical and the heats of the vibration 
terms to be given by the Einstein formula, the 
correction is +1.3 kcal. Then AH) =138.1 kcal. 
In other units AH) is 6.00 electron volts. Since this 
calculation is based on an upper limit to the 
concentration of CN found previously, it only 
yields a lower limit to the heat of dissociation. 
To decide how close the actual value is to this 
limit, we must consider the method of deter- 
mining it. 


DEGREE OF DISSOCIATION OF CYANOGEN IN 
THE ELECTRIC DISCHARGE 


We must now discuss the assumption of com- 
plete dissociation of the cyanogen in the electric 
discharge at great dilution, which was made in 
Part I of this paper.’ The first evidence comes 
from a study of the mechanism of the discharge. 
From Table III, which is reproduced from 
Part I, it is evident that at the start of a dis- 
charge in a mixture of cyanogen and argon the 
process requiring least energy is the formation 
of CN radicals, thereafter it is their excitation. 
The low excitation potentials of CN explain 
why only its spectrum was seen in emission. 
In the case where there is very little cyanogen 
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present one would expect it to be almost com- 
pletely dissociated, which would be indicated 
by a linear relation between the amount of 
cyanogen present and the intensity of CN 
absorption. This has been found. The amount of 
cyanogen needed to supply the ions postulated 
in the section on reaction kinetics is less than 
one percent of the total amount present at the 
greatest dilution when any rates were measured. 

Since the publication of Part I,! it has been 
pointed out to the author that as soon as some 
CN radicals have been formed, their dissociation 
to carbon and nitrogen is energetically about as 
easy as the dissociation of the remaining cyano- 
gen to form more radicals. Therefore, even when 
the cyanogen present has been very largely 
dissociated, the number of radicals actually 
present is less than it would be if there were no 
secondary dissociation by an amount depending 
on the relative collision cross sections of C2N» 
and CN. Lacking experimental knowledge of 
these cross sections, we can only assume from 
their energy requirements that they are about 
equal. In that case, a calculation of the reaction 
rates shows that the maximum concentration of 
CN attained is 1/e of that predicted for no 
secondary dissociation. As the measurements 
were probably not made exactly at the time of 
maximum concentration, there is a further, but 
unknown correction to be applied in the same 
direction. 

We are guided by three other considerations 
in deciding how much allowance to make for 
this. The first is rigorous and sets a lower limit 
to the concentration of CN that could have 
been present. Dispersion theory states that for 
any molecule the sum of the f values for all the 
transitions out of a given state is one. Therefore, 
in the case of the greatest possible error the 
concentration of CN cannot be less than 0.026 
of the calculated value, because in that case the 
f value for that one transition would be greater 
than unity. Since CN has other transitions, the 
lower limit must be somewhat higher than this. 
The second check lies in the comparison of the 
f value with that calculated by Mulliken*® from 
theoretical considerations. His value, which is of 
the order of 0.1, indicates that not all the CN 


6 R.S. Mulliken, private communication with the author. 
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that could have been formed was observed in 
the present experiment. The third is that the 
observed dissociation energy lies slightly on the 
short wave side of the beginning of the absorption 
spectrum of cyanogen. From a study of the 
photopolymerization of cyanogen Hogness and 
Liu-Sheng T’sai? have found that it is polymer- 
ized appreciably by the light from a hydrogen 
tube transmitted through air. To get an appreci- 
able effect, light of wave-lengths considerably 
longer than 1850A must be effective. Their 
conclusion from their measurements, that the 
first. step in polymerization is dissociation to 
free radicals, is confirmed by the calculations 
made in the section on thermal polymerization. 

In view of these facts it appears that the 
f value of 0.026(+0.006) calculated previously 
is too low. It remains as an experimentally 
observed lower limit, but the best value is 
probably three or four times as large. Corre- 
spondingly, the heat of dissociation given above 
should be taken as an accurate lower limit. 
f=0.1(+0.05) is used in calculating the most 
probable value. Since this method of calculation 
is insensitive to variations in the f value, this 
change only increases Hy by 8 kcal. to make the 
most probable value AH)=146(+4) kcal., or 
in other units AHy)=6.33(+0.17) volts. 


COMPARISON WITH PREVIOUS WORK 


Kistiakowsky and Gershinowitz® found AH 
=77 kcal. by measuring the variation of K with 
temperature on a quartz spectrograph. It is the 
purpose of this paragraph to show that their 
incomplete resolution and neglect of pressure 
broadening effects may introduce large enough 
errors to account for the difference. Their 
fundamental measurement was the comparison 
of the intensities of absorption of the band head 
at 43883 for different temperatures and partial 
pressures of cyanogen. Effectively they measured 
K in arbitrary units and substituted in the 
van’t Hoff equation. For example, if K and K’ 
correspond to the temperatures of JT and 7” 
respectively, 





"T. R. Hogness and Liu-Sheng T’sai, J. Am. Chem. 
Soc. 54, 123 (1932). 

*G. B. Kistiakowsky and H. Gershinowitz, J. Chem. 
Phys. 1, 432 (1933). 
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If the value of K’ defines the arbitrary units in 
which it is expressed, all the error in the deter- 
mination must appear in K. The error in K 
needed to produce the observed difference in AH 
is found by differentiating and dividing by the 
original equation. If we use the lower limit 
found for AH, 


6InK 6AH 137-77 
en =0.44. 
InK/K’ H 137 





Then 
6K /K =0.44 In K/K’. 


For T—T’=46°, which is larger than most of 
their temperature differences, K/K’=4. Since 
K=[(CN }/[C2N2], and [C2Ne] is accurately 
measurable, 


46K /K =8[CN ]/[CN]=0.30. 


Thus, an error of only 30 percent in estimating 
the equality of two absorption lines would 
account for the discrepancy. Their measurements 
were all made at atmospheric pressure, in a 
furnace about a meter long, using a spectrograph 
of medium resolving power. From the present 
work with very high resolving power I have 
calculated that the absorption had to become 
very strong before they could observe it at all. 
All their measurements lie in the region where 
the centers of the lines are almost completely 
absorbed out. There the amount of light absorbed 
changes slowly with the concentration of the 
absorbing molecules but is very sensitive to 
broadening of the lines. If the CN band lines are 
broader in an atmosphere of cyanogen than in 
one of nitrogen, it is understandable that halving 
the partial pressure of cyanogen and substituting 
nitrogen for it might decrease the width of the 
lines by 30 percent. This would explain the 
difference between their results and mine. 


CHEMICAL APPLICATIONS 


Thermal polymerization 


The change of free energy with temperature of 
any reaction at atmospheric pressure is given by 
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TABLE IV. Possible heats of dissociation of CO and CN and heats of vaporization of graphite on the basis of predissociations 
observed in the CO spectrum. 








DISSOCIATION PRODUCTS 
oF CO at 11.054 voLts 
CARBON OXYGEN 


D(CO) 
VOLTS 


D(CN) 


VOLTS 


KCAL. 





‘p sp 
1D ad 
a 1D 
1S sp 
1p 1p 
sp 1S 
5S sp 


11.054 
9.797 
9.097 
8.384 
7.840 
6.886 
6.738 


6.38(+0.13) 
5.12 


147(+3) 
118 


4.42 
3.71 
3.17 
2.21 
2.06 








the partial differential equation 
(6AF/5T),=—AS. 


The entropy change in the dissociation of 
cyanogen has, as explained above, been calcu- 
lated as a function of temperature. When an 
approximation, valid in the neighborhood of 
1200°K is made for the vibrational entropy 
change, the following expression is obtained: 


AS=4.56+7/2R In T—0.013T. 


Substituting this in the equation above and 
integrating at constant pressure, an expression 
is obtained for AF. The constant is evaluated 
by using the known free energy change at 
1500°F. From the values of AF around 400°C, 
the order of magnitude of the partial pressure 
and density of CN radicals is calculated. At 
375°C, there is 1/30 of a radical per cc, at 
425°C, 3, and at 475°C, 300. If each one formed 
adds immediately to a cyanogen molecule, the 
polymerization process should begin to become 
appreciable somewhere in this temperature range. 
Actually, it is observed to start slightly above 
400°C, which is in satisfactory agreement with 
the calculations made above. This is further 
evidence that the first step in polymerization is 
the formation of free radicals and that at 
temperatures where polymerization occurs, free 
radicals of CN may add to cyanogen directly 
instead of recombining with each other. 


HEAT OF VAPORIZATION OF GRAPHITE AND HEAT 
OF DISSOCIATION OF CO 


To show the application of the present work 
to the determination of these two quantities, 
it is necessary to give an outline of the other 


relevant data. Herzberg® has given a very clear 
and detailed treatment of the problem in his 
review of this field, from which the important 
points in the following discussion are taken. 
By two thermodynamic cycles involving data 
that was incontestible, except for the heat of 
dissociation of cyanogen to two CN radicals, he 
obtained two relations between L, the heat of 
vaporization of graphite, D(CO), and D(CN). 
After the one involving D(CN) has been changed 
to correspond to the new value for cyanogen, 
they are: 
L=D(CO) —85.6(+0.5) kcal. 
D(CN) = D(CO) —108(+3) kcal. 


Thus, if any one of these three quantities can be 
determined, the other two follow directly. The 
most accurate information is that based on the 
predissociations observed in the CO spectrum. 
As Herzberg has pointed out, the predissociation 
limit at 11.054 volts is an extremely accurate 
value for a possible heat of dissociation. Un- 
fortunately it is not necessarily the correct value. 
The data only tells us that the heat of dissocia- 
tion has this value if the dissociation products 
are not excited, otherwise it is this value minus 
the excitation energy of the products. As the 
position of the low lying levels of carbon and 
oxygen are well known, it is possible to calculate 
the dissociation energies corresponding to all the 
permissible combinations of products. Herzberg’s 
table of all these quantities compatible with the 
known vibration levels of CO is reproduced in 
Table IV with the changes necessitated by the 
new cyanogen value. The values of L and 
D(CN) corresponding to each possibility are also 
given. Unfortunately there is not enough evi- 
dence known from the CO spectrum to determine 


9G. Herzberg, Chem. Rev. 20, 145 (1937). 
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which possibility is the correct one. The selection 
must be made from supplementary data about 
CN or L. 

The heat of dissociation of CN has not been 
determined accurately. However, the vibration 
levels of the lowest *2 and °II states have been 
observed directly up to 3.4 volts and those of 
the *II have been clearly followed by using the 
perturbations up to 4.2 volts. As all the evidence 
indicates that these two states dissociate to the 
same products, the heat of dissociation must be 
larger than this. A Birge-Sponer extrapolation 
of the levels gives around 7 volts for the heat of 
dissociation of the normal state, which is almost 
certain to be too high. It would, however, be 
rather difficult to reconcile it and the known 
levels with a value less than 5.5 volts. 

If CN dissociates to normal atoms, C(*P) 
+N(4S), the highest possibility of Table IV is 
the only acceptable one, although the next is so 
close to the arbitrary limit we have set that it 
cannot be completely excluded. If CN goes to 
the first permitted combination of excited atoms, 
1.35 volts above the ground state, the second, 
third, and fourth possibilities are acceptable. 
The high excitation energy of all the other 
combinations of products makes them incom- 
patible with the known vibration levels of CO. 
The possible energies for the two pairs of 
products just mentioned can be further restricted 
by considering the dissociation energies allowed 
the B® state of CN on the basis of the possible 
dissociation products remaining for it. In the 
first case the highest value gives a very reason- 
able heat of dissociation for the excited state, 
but the next highest is too low. This corroborates 
the conclusions drawn from the ground state 
and further decreases the likelihood of dissocia- 
tion to normal atoms at 5.12 volts. In the 
second case, the third and fourth highest values 
are only compatible with the dissociation of the 
excited state to combinations of products two 
and three pairs above the lowest possible. The 
usual rule that a state dissociates to the lowest 
permissible pair of products indicates that these 
two possibilities are unlikely compared to the 


second. By the same rule, dissociation to excited 
atoms at 5.12 volts is less likely than to normal 
ones at 6.38 volts. 

The effect of errors in the f value on these 
arguments should be considered. The greatest 
errors possible, corresponding to f=1 and 
f=0.026 increase and decrease AHoc,n, by 0.60 
and 0.33 volt, respectively, which would decrease 
or increase the values of D(CN) by half as much. 
The estimated probable error is 0.13 volt. All 
the possibilities mentioned above are permissible 
throughout this range except the third, which 
is impossible for the smaller f values. Dissociation 
to the first combination of normal atoms is 
always the most acceptable. In summarizing, 
we may say that if CN dissociates to normal 
atoms, the most reasonable value is 6.38 volts, 
and if it dissociates to excited atoms, 5.12 volts 
is the best one. As Sponer!® and Mulliken" both 
give CN normal dissociation products, the most 
likely possibility is D(CN) =6.38 volts, D(CO) 
= 11.054 volts, L=169.2 kcal. 

This value of L is in excellent agreement with 
the thermal measurement of 177 kcal. made by 
Marshall and Norton” and with almost all the 
other evaporation method values. It does not 
agree with 150 kcal., the only equilibrium rather 
than dynamic value.” The discrepancy between 
these has been ascribed to experimental error 
and to a real difference between the equilibrium 
and dynamic vapor pressures. While the value 
chosen as most likely on the basis of the present 
work indicates that it is experimental error, the 
fact that the second choice gives a value of 
L=140.2 kcal. prevents us ffom making a posi- 
tive selection between these two. 

The author is indebted for helpful discussions 
and criticisms to Professor F. A. Jenkins, at 
whose suggestion the problem was undertaken, 
and to Professor G. K. Rollefson. 


10H. Sponer, Molekulspecktren, p. 39. 

1 R. S. Mulliken, Rev. Mod. Phys. 4, 48 (1932). 

22 A. L. Marshall and F. I. Norton, J. Am. Chem. Soc. 
55, 431 (1933). 

13H. Kohn and M. Guckel, Zeits. f. Physik 27, 305 
(1925). 





JOURNAL OF CHEMICAL PHYSICS 


VOLUME 8 


The Photolysis of Methylethyl Ketone 
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The nature of the products from the photochemical decomposition of methylethyi ketone at 
temperatures from 90 to 200° has been re-examined. In contrast with earlier results, the forma- 
tion of ethane, propane, and butane parallels the frequency of the corresponding free radical 
collisions. At the higher temperatures, however, production of ethane and especially of methane 
becomes predominant. The influences on the product distribution of secondary reactions, the 
lifetimes of the free radical collision complexes, and the activation energies and entropies, are 


considered. 





HE photochemical decomposition of methyl- 

ethyl ketone has previously been investi- 
gated by Norrish and Appleyard! using radiation 
in the near ultraviolet absorption region around 
3000A, and more recently by Ells and Noyes? 
who worked in the region around 1850A. Both 
these investigations suggested that equivalent 
amounts of ethane, propane, and butane occurred 
in the reaction products, together with an 
amount of methane which increased markedly 
with increasing temperature. This result was 
indeed surprising, since if the hydrocarbons were 
formed by interaction of free radicals produced 
in the primary process, one would expect the 
ratio ethane : propane : butane to be approxi- 
mately 1 : 2:1. The present experiments were 
undertaken to check this anomalous result, using 
a method of analysis for hydrocarbons which 
had proved reliable in previous studies of free 
radical reactions.’ 


EXPERIMENTAL DETAILS 


The photolyses were carried out in a static 
system identical with that used for our studies 
of metal alkyl decompositions.* Ketone at forty 
mm pressure was irradiated with the unfiltered 
output of a high pressure arc of the mercury- 
in-quartz type, operating hot at 55 volts and 
3.5 amp. d.c., until about fifty percent was 
decomposed. The methylethyl ketone used 
was a product of the Eastman Kodak Company, 
purified by treatment with anhydrous copper 
sulfate, and careful fractional distillation. A 


1 Norrish and Appleyard, J. Chem. Soc. 874 (1934). 


2 Ells and Noyes, J. Am. Chem. Soc. 60, 2031 (1938). 
3 Moore and Taylor, J. Chem. Phys. 8, 396 (1940). 


fraction boiling within a tenth of a degree at 
around 79.5° was retained; its refractive index 
was mp'=1.38118, compared with a reported 
value for the pure compound! of mp = 1.38140. 

The products, amounting to about five cc of 
gas, were analyzed as follows: 

1. A fraction pumped off from liquid air was 
burned over cupric oxide at 310°, converting 
carbon monoxide to dioxide, from which methane 
could be pumped off at liquid-air temperature; 
the carbon dioxide was pumped off at — 78°, and 
any residual vapor was assumed to be water. 

2. A fraction pumped off at —138° (m.p. of 
ethyl chloride) was hydrogenated on nickel at 
100° to determine ethylene, then cracked over 
nickel at 310° to methane, thus giving the aver- 
age carbon content of the fraction. This fraction 
is made up of C2 and C3; hydrocarbons. 

3. A fraction at —78°, consisting of C3 and Cy 
hydrocarbons, was similarly analyzed. 


EXPERIMENTAL RESULTS 


The observations were extended over a 
temperature range of from 90 to 200°. At the 
lowest temperature there was observed a small 
amount of diketone formation, evidenced by the 
appearance of a yellow oil and ‘the typical odor 
of diketones in the reaction vessel. At the higher 
temperatures, however, no diketone formation 
was detected. Tests for acetaldehyde with 
Schiff’s reagent were negative. There was evi- 
dently some polymer formation at all tempera- 
tures, and after pumping off products volatile 
at —78°, an odor resembling that of mesityl 


4]. Timmermans and F. Martin, J. chim. phys. 25, 411 
(1928). 
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oxide was distinctly noticeable ; there was insuff- 
cient residue, however, to make possible the 
identification of the responsible compound or 
compounds. 

The results of the analyses of the volatile 
products are given in Table I. 

We note that at the lower temperatures 
considerably more propane is formed than either 
ethane or butane, but that at 200° the formation 
of methane and ethane is of primary importance. 


GENERAL DISCUSSION 


If, as a first approximation, we disregard the 
reactions of acetyl and propionyl radicals 
possibly arising in the primary process, the 
reactions with which we are concerned are as 
follows: 


CH3+CH3=CoH,g, (1a) 
CH3;+C.H;=C3Hs, (1b) 
C3:H5+C2H;=CaH no, (1c) 


CH3;+CH;COC.H; =CH,+CH3;COC2H,, (2a) 
CsH;+ CH 3COC>2H; = C.Hs+CH;COC:;H 4: (2b) 


At 110° the hydrocarbon products probably 
arise principally from the interactions of free 
radicals. The ratio of the numbers of collisions 
of types la : 1b : 1c is roughly 1.00 : 1.74 : 0.83, 
while the observed formation of ethane, propane, 
and butane are in the ratio 1.00 : 1.54 : 1.00. 
Assuming equal activation energies for the 
reactions (1), there are two other factors which 
may influence the relative recombination rates 
of the free radicals. 

The lifetime of the quasi-molecule formed in 
the free radical collisions will increase from (1a) 
to (1c). The entropy of activation should 
decrease in going from (1a) to (1c), since it is 
less for reactions between more complex radicals.® 
A proper evaluation of these factors would 
perhaps suffice to explain the observed deviations 
from the purely kinetic theory ratios. 

In addition, and especially at higher tempera- 
tures, the secondary reactions (2) will be of 





* Bawn, Trans. Faraday Soc. 31, 1536 (1935). 


OF METHYLETHYL KETONE 





TABLE I. Percentage composition of photolysis 
products of MeEtCO. 











Temp. °C 90 110 155 200 
He. 0.4 0.1 0.3 0.0 
CO 47.5 49.0 47.0 48.5 
CH, 3.4 4.4 9.0 18.0 
CH, ye 0.5 1.7 3.0 
C2H¢ 11.5 13.0 14.5 12.5 
C;Hs 19.1 20.0 18.6 11.0 
CsHio 13.2 13.0 9.0 7.0 
C,H, 1.5 — 0.0 — 
























importance, leading to the formation of methane 
and ethane at the expense of the other hydro- 
carbons. A detailed kinetic analysis of these 
reactions is difficult and not very informative, 
since we can already follow qualitatively the 
effect of the decreased radical concentrations in 
decreasing the amount of recombination of 
radicals as compared with secondary reaction 
with ketone molecules. 

It is interesting to note that methane is 
formed from methylethyl ketone considerably 
more easily than from acetone.® This difference 
may most reasonably be ascribed to the lowered 
energy of activation for the reaction of a methyl 
radical with —C.H; to form C:H,4, as compared 
with the reaction with —CH; to form CHg. 

It appears also that methyl radicals are more 
effective in polymerizing ethylene than are 
ethyl radicals, since much greater amounts of 
ethylene remain in the photochemical decompo- 
sition of diethyl ketone than are found with 
methylethyl ketone.2” This is in accord with 
the results of experiments by Yeddanapalli 
and Jungers* on the relative efficiencies in 
ethylene polymerization of methyl and ethyl 
radicals prepared by reaction of the correspond- 
ing iodides with mercury vapor. 

The formation of an oxygen-containing poly- 
mer from such fragments as CH3COC2Hy, would 
also account for the small amount of ethylene 
formation, this possibility being also suggested 
by the characteristic odor of the nonvolatile 
photolysis products. 

6 Taylor and Rosenblum, J. Chem. Phys. 6, 119 (1938). 

7 Noyes and Ells, J. Am. Chem. Soc. 61, 2492 (1939). 


8 J. C. Jungers and L. M. Yeddanapalli, Trans. Faraday 
Soc. 36, 483 (1940). 
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The Thermal Reaction Between Hydrogen and Oxygen 


I. Comparison of Experiments and Theory 
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The present theory of the slow thermal reaction between hydrogen and oxygen, while 
satisfactorily describing the existence of the various explosion limits, fails to describe the 
following features: the sudden transition from the slowest reaction to the explosion (‘‘second 
explosion limit’’), and some effects of the vessel size and of a KCI film on the rate and the 
second explosion limit. Conclusions regarding the theory are discussed. 





I. PROBLEM 


OR certain spectrographic investigations 
which we shall report on in a subsequent 
paper, it is desirable to produce the maximum 
steady reaction rate of the thermal reaction of 
hydrogen and oxygen. As our work requires a 
large reaction vessel and pressures near atmos- 
pheric, the danger of a violent explosion prompted 
us first to study the explosion limits in smaller 
containers. 

Hinshelwood and his co-workers! have de- 
scribed this reaction by means of the pressure- 
temperature diagram of Fig. 1. Explosion occurs 
for any combination of pressure and temperature 
located on the peninsula of the diagram; above 
and below this peninsula a steady reaction 
takes place. (On the scale of Fig. 1, the ‘“‘Con- 
tinuous Reaction” below the ‘First Limit” of 
the explosion peninsula does not show up.) 

(Hinshelwood, investigating the explosion pen- 
insula, called its borderlines ‘‘upper”’ and ‘‘lower”’ 
explosion limits. In our discussion the name 
“upper” limit would be misleading, since our 
attention is focused on the range of the steady 
reaction between this ‘“‘upper” and the still 
higher ‘‘third”’ limit (Fig. 1). Therefore we shall 
use the term “second” instead of the conven- 
tional ‘‘upper”’ explosion limit.) 

While older theories restricted themselves to 
rather vague outlines of the part played by the 
gas phase and the surface reactions, recently a 
detailed reaction scheme has been worked out by 


1See the book by C. N. Hinshelwood and A. T. Wil- 
liamson, The Reaction between Hydrogen and Oxygen 
(Oxford, Clarendon Press, 1934). 
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Kassel and Storch, and Lewis and von Elbe.’ 
This theory accounts for a great many striking 
features of the reaction, in particular the various 
explosion limits and the steady reaction between 
the second and the third limit. In our experi- 
mental work we found many results explained by 
the theory, but also other results which we failed 
to reconcile with it. It is the purpose of the 
present paper to discuss whether this theory may 
need some modification. 


II. THE PRESENT STATE OF THE THEORY® 


There is general agreement that the reaction 
above the second explosion limit, at sufficiently 
high temperatures, is a chain reaction and that 
the walls play a part in breaking the chains. 
This hypothesis is necessary in order to explain 
the fact that the rate increases with the vessel 
size or with the addition of a foreign gas, which 
is expected to slow down the diffusion of the 
chain carriers to the walls. 


2L. S. Kassel and H. H. Storch, J. Am. Chem. Soc. 57, 
672 (1935); L. S. Kassel, Chem. Rev. 21, 331 (1937); 
B. Lewis and G. von Elbe, J. Am. Chem. Soc. 59, 656 
(1937); B. Lewis and G. von Elbe, Combustion, Flames 
and Explosions of Gases (Cambridge University Press, 
1938). See the recent textbook of H. J. Schumacher, 
Chemische Gasreactionen (Dresden and Leipzig, Theodor 
Steinkopff, 1938), p. 462. In some respect the most recent 
form of the theory, carried through in full detail by Lewis 
and von Elbe, differs from the earlier form, summarized 
by Kassel in Chem. Rev.; Kassel assumes that the chains 
ordinarily start on the walls with but a few originating in 
the gas phase, while Lewis and von Elbe find that the 
starting of most chains at the walls is incompatible with 
their scheme (book, p. 68). 

3 A detailed discussion of the theory is given in the book 
of Lewis and von Elbe, reference 2. In order to avol 
confusion we apply the same numbers to the single 
reactions. 
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THERMAL REACTION 


The detail of the theory is based on Semenoff’s* 
conception of branching chains whose breaking 
exceeds the branching in the steady reaction and 
just equals the branching at the explosion limits. 
This is accomplished by the following chain: 


OH+H.—-H,0+H, I 
H+0:—-0H+0, Il 
O+H.-OH+H. III 


Since the branching reactions II and III are of the 
first order in the chain carriers, the breaking 
reaction is, to a certain extent, restricted by the 
theorem of Kassel and Storch.” This states that 
branched chain explosions can occur only if 
chain branching and chain breaking processes 
are of the same order with respect to the chain 
carriers. This requirement is fulfilled by the 
breaking reaction: 


H+0.+M—HO.+M. VI 


However, VI actually breaks the chain only if 
the radical HO which is formed by VI is sure 
to be destroyed at the wall: 


HO.+ wall—destruction. XII 


The chain breaking reaction VI requires triple 
collisions. Therefore it becomes increasingly rare 
with decreasing pressure. This failure of chain 
breaking at low pressures explains the second 
explosion limit on the low pressure side of the 
steady reaction. 

On the other hand, it is observed that at 
higher pressures the rate of the steady reaction 
increases with the pressure so strongly that for a 
certain high pressure limit another type of 
explosion results (third explosion limit). In order 
to explain this it is assumed that at higher 
pressure VI is not absolutely chain breaking, 
but that instead HO: has a chance to continue 
the chain by 


X must be supposed to have a high energy of 
activation so that it comes into play only after 
numerous collisions but is ineffective at low 
pressures because of too rapid a diffusion of HO: 
to the wall. (The alternative continuation of the 
chain by HO.+H2—-H.,0.+H is disregarded 


*N. Semenoff, Chemical Kinetics and Chain Reactions 
(Oxford, Clarendon Press, 1935). 
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Fic. 1. Explosion limits of reaction of 2H»+Os. First 
and second limits after Hinshelwood. Dotted section never 
investigated. 


only because it would not materially affect the 
final result.) 

The chain is started by the thermal dissocia- 
tion of He in the gas phase 


H2+M—-H+H+M. XIII 


This set of reactions successfully describes a 
very complicated behavior; the steady reaction 
at varying pressure bordered on both sides by 
explosions. 

By means of a detailed computation Lewis and 
von Elbe explain the shape of the curve of the 
second explosion limit (Fig. 1) and the effect of 
varying the partial pressures of Ho, Os, and 
inert gases.® 

In a letter to the editor,® we reported meas- 
urements of the third explosion limit, which also 
are in substantial agreement with the predictions 
given by the present theory. 


III. DIFFICULTIES 


While this theory explains the facts presented 
in the preceding section, it leads to some diffi- 
culties’ with regard to other experiments which 


5 Lewis and von Elbe identify the shape of this curve 
with the temperature dependence of reaction II (book, 
pp. 32 and 35). So they attribute an energy of activation 
of 24.2 to 26.0 kcal. per mole to this reaction. This agrees 
sufficiently well with the value of 32 kcal. per mole derived 
by Kantorowa and Nalbandjan on the basis of entirely 
“gaa experiments (Physik. Zeits. Sowjetunion 4, 758 

1933)). 

6Q. Oldenberg and H. S. Sommers, Jr., J. Chem. 
Phys. 7, 279 (1939). 

7As an outstanding difficulty Lewis and von Elbe 
(book, p. 68; J. Chem. Phys. 7, 711, footnote 10 (1939)) 
mention the surprisingly strong effect of a silver surface 
in suppressing the reaction. Recently they succeeded in 
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limit. Pyrex vessel coated with KCI. T=568°C. 


will be discussed in the present section. Our 
disagreement with the theory is largely based on 
recent experiments which were not available to 
its authors. 


A. The rate-pressure curve 


An effect difficult to interpret is manifest in 
the rate-pressure curve (constant temperature, 
varying initial pressure, the rate measured either 
at the beginning of the reaction or, instead, after 
the inhibition period). Schematically such a 
curve was given by Thompson and Hinshelwood.® 
It shows a gradual decrease of rate with diminish- 
ing pressure to a very low value, suddenly 
terminated by the explosion (second explosion 
limit). As the abruptness of this change is very 
striking we checked it in greater detail. (We shall 
describe our experimental arrangement in a paper 
on the third explosion limit, to be published. 
For low pressure the mixture 2H.+Oz cannot be 
directly admitted into an evacuated vessel be- 
cause it would explode during admission ; in order 
to prevent this a total of about 50 cm was 
admitted, beginning with He, and then quickly 
pumped ‘ut to the desired pressure.) The re- 
sulting rate-pressure curve is given in Fig. 2. 
The minimum rate was only 0.1 cm/min. The 
experiments were made in a KCl-covered Pyrex 
vessel, 14 cm long, 4.5 cm wide, at 568°C. 

The same behavior was observed, without 
accurate measurements, in very many experi- 


explaining this effect themselves by the presence of silver 
in the vapor phase, exceeding the regular vapor pressure, 
so that this difficulty has been removed. 

8H. W. Thompson and C. N. Hinshelwood, Proc. Roy. 
Soc. London, A122, 610 (1929). 
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ments with clean Pyrex as well as Pyrex covered 
with KCI: at a certain range of low pressures the 
reaction is very slow indeed, the more so the 
lower the pressure, while at a pressure slightly 
lower still, below a sharply defined limit, the 
second limit, explosion occurs. 

This second explosion limit is reproducible, 
within a mm, in a large group of successive 
experiments in the same vessel. In some types 
of vessel it may depend on the rate of manipula- 
tion (see Section III D2). There is no conspicuous 
induction period ;° as soon as gradual pumping 
brings the pressure to the limit, explosion occurs; 
if in another experiment the pumping is stopped 
a millimeter above the limit, an exceedingly slow 
reaction follows which never goes over into an 
explosion. 

The explosion at the second limit is certainly 
not caused by heating, as it takes place after a 
minute change of pressure from another con- 
dition at which the rate is exceedingly low and, 
therefore, the heat production negligible. The 
sudden change of rate is explained, in general 
terms, by Semenoff’s* theory of branching chains 
as follows. The explosion occurs as soon as the 
branching exceeds the breaking in one chain. 
If we assume the breaking to be caused by 
triple collisions, the average length of the indi- 
vidual chain increases with decreasing pressure 
until finally one of the chains gets out of control 
and leads to the explosion of the second limit. 
This picture, in which the explosion is attributed 
to the development of an individual chain, could 
be compatible with an over-all yield decreasing 
with decreasing pressure right above this limit 
(Fig. 2) only if we assume that the other factor 
affecting the over-all yield—the number of 
chains per sec.—decreases very strongly with the 
pressure. 

While this general picture of the second ex- 
plosion limit may seem satisfactory, the specific 
picture of the chain reaction given in the present 
theory must furnish such chains (number per 
sec., continuation, branching, breaking) as a 
function of the pressure as just described in the 
general picture. In the theory, the rate dH,0/di 


* A very brief induction period of 0.1 sec. was discovered 
by A. Kowalsky (Physik. Zeits. Sowjetunion 1, 595 (1932); 
4, 723 (1933)) with the help of a photographically recorded 
membrane manometer. 
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THERMAL REACTION 


is expressed by a complicated equation (66 in the 
book by Lewis and von Elbe)." It is claimed that 
for any combination of partial pressures of 
hydrogen, oxygen, and a foreign gas, this equa- 
tion represents the rate, in particular, in the 
various limiting cases in which the rate may 
approach infinity, that is, the explosion limits. 
For the second explosion limit this equation is 
simplified by neglecting the chain continuation 
through the radical HOs, that is, discarding the 
reaction X: HOs+H2—-H2,O+OH. The reason 
for discarding this reaction is that if HO» should 
be able to continue the chains, the second 
explosion limit would depend upon the vessel 
size. This would be contradictory to the observa- 
tions of Hinshelwood. Another obvious simpli- 
fication of the general formula is given by the 
fact that we are dealing with 2H2+O, without 
foreign gases (total concentration M=2H,.+O.). 

The explosion limits are derived from the 
equation for the rate by the requirement that 
the rate become infinite. As a matter of fact, the 
rate is derived from the theory as a smooth 
function of the pressure curving up, as expected, 
on both sides toward the vertical (Fig. 3). The 
minimum rate (caused by the counteraction of 
increasing chain length and decreasing number 
of chains) is derived from the same formula; it 
occurs at a pressure twice as high as the second 
explosion." Hence the theoretical curve does not 
contain the sharp kink given by the experiments 
at a low value of the rate near the minimum of 
the curve (dotted line of Fig. 3). 

To this argument the objection has been made 
that the left wing of the curve discussed (left of 
the dotted line) represents unstable states. This 


10 Although Lewis and von Elbe state it as their im- 
mediate aim to discuss experiments from 560°C upward, 
any explanation of the second explosion limit would be 
unsatisfactory which would fail to apply to the whole 
temperature range in which this type of explosion occurs, 
that is, down to 450°C. According to Hinshelwood and 
Thompson, the high order gas reaction comes into promi- 
nence already from 520-530°C; this applies to a con- 
ventional vessel size. 

"This minimum is computed with the same approxi- 
mation, neglecting reaction X. It is true that this must 
be applied only to the second explosion limit, not neces- 
sarily to higher pressures. But since the theory claims to 
represent the real curve in which minimum and second 
explosion limit nearly coincide, no such distinction must 
be made. Moreover, it is obvious that by discarding the 
approximation we may obtain a shifted position of the 
minimum, but by no means the sharp kink in the curve 
as it is actually observed. 
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means that the curve accounts only for the 
initial reaction rate of the mixture; by the 
progress of the reaction the partial pressure of 
the reactive mixture decreases, so that the rate 
increases until the explosion occurs. 

This intrepretation of the present theory seems 
to us not valid. First of all, the initial rate for a 
pressure slightly above the second explosion limit 
is so very low that even a considerable increase 
of rate would not affect the stability. A reliable 
criterion for the stability, however, is derived 
from the general equation for the reaction rate 
d(H2O)/dé (Eq. (66) in the book). This equation 
is supposed to express the rate for any concen- 
trations of He, Oz, and a foreign gas (the latter in 
this case would represent the H,O gradually 
developed by the reaction). But now, in contrast 
to the preceding application of the same formula, 
we want to discuss the progress of the reaction 
(not only the initial rate for which M =2H2+0Oz). 
Therefore, the variables Hz, Oo, and M are con- 
nected by the obvious relation expressing the 
gradual decrease of M by the reaction, 2H.O 
replacing 2H2.+O:. Furthermore, dealing with 
the second explosion limit, we neglect reaction X 
(see above). 

The result of a simple computation is that for 
the initial pressure which belongs to the mini- 
mum of the curve in Fig. 1, the rate not only 
remains finite throughout, but actually fails to 
increase during the reaction. This result of the 
theory contrasts with the experiment in which 
the explosion occurs right in the beginning of the 
reaction for conditions very close to those of an 
exceedingly slow reaction. Hence the theoretical 
formula does not lead to the abrupt change which 
is so striking in the experiments.” 











Rate 





i@) 2 4 6 
Pressure ~(units arbitrary) 


Fic. 3. Rate versus pressure. The solid line represents 
the theory of Lewis and von Elbe. The broken line indicates 
the pressure where the second explosion actually occurs. 


12So far we have considered only the case of chain 
initiation in the gas phase, which is the case principally 
treated by Lewis and von Elbe. These authors also 
consider chain initiation at the wall (their Eq. (82)). 





472 


Another possibility for reconciling the theory 
with the abrupt transition from a very low rate 
to the explosion is as follows: In the theory, the 
chains get longer under decreasing pressure 
because the breaking requires triple collisions. 
It may happen that an individual very long 
chain develops a considerable amount of heat 
within a small volume and so causes the explosion. 
Such a local concentration of one chain can 
occur only if all steps of the chain take place 
without many intervening collisions, which would 
cause diffusion of the carriers and distribution of 
the heat over a wider volume. Hence a high 
heat of activation of one or the other step would 
militate against such local origin of the second 
explosion. The heat of activation of II is supposed 
to be high (30,000 cal./mole, see Section II). 
For such a high value, local heating is incom- 
patible with the fact that the second explosion 
can be traced down to temperatures as low as 
450°C where the consecutive steps of the chain 
would not be expected to follow each other 
without very many intervening collisions. The 
result is that we are unable to reconcile the 
present theory (chain breaking by formation in 
triple collisions of the radical HO.) with the 
abruptness observed in the change from the 
lowest rates to the explosion. 

This argument is directed against the specific 
picture of the reaction discussed in Section II. 
One still may hope to preserve the picture of the 
chain proper, that is, its continuation, breaking, 
and branching, and only mend the theory by 
introducing another starting process which would 
decrease very rapidly with the pressure. Although 
it may be mathematically possible to invent such 
a pressure function, it seems safe to guess that 
it would not represent a reasonable process of 
chains starting in the gas phase. 


B. Effects of vessel size on second explosion 
limit and rate 


Another difficulty arises from the effect of the 
vessel size on the second explosion limit and the 
rate. The second explosion limit is observed to be 


But after a thorough mathematical treatment they reject 
this hypothesis (p. 68) for several reasons. Anyway, it is 
of interest that this modification of their theory (although 
not so easily treated mathematically) does not lead either 
to the striking sharp kink near the minimum pressure 
which is actually observed. 
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independent of the vessel size." Theoretically, 
this leads to the conclusion that, under the con- 
ditions of the second explosion limit, reaction X 
(continuation of the chain by HOz) is rare as 
compared to XII (destruction of HOz at the 
wall). Apart from the mathematical reasoning 
given by Lewis and von Elbe (book, p. 42), the 
exclusion of reaction X under the conditions of 
the second explosion limit may be understood as 
follows: According to the theory, the explosion 
is essentially a process in the gas phase. The only 
wall effect appearing in the reaction scheme 
(Section II) is the destruction at the wall of the 
HO: radicals. The observation that the second 
explosion limit is independent of the vessel size 
leads to the conclusion that the HOs is certain to 
be destroyed at the wall, irrespective of whether 
it travels a longer or shorter distance before it 
reaches the wall. Hence any chance for HO, to 
continue the chain must be relatively rare, that 
is, reaction X must be negligible under the 
conditions of pressure and temperature charac- 
teristic of the second explosion limit. 

This explanation, however, necessarily involves 
an effect of the vessel size on the rate as is seen by 
the following argument. We understand that 
under the conditions of the second explosion 
limit the triple collision VI is absolutely chain 
breaking, since reaction X is excluded, inde- 
pendent of the vessel size. Hence the reaction X 
must be excluded as well at a pressure only very 
slightly higher, that is, in the region of the 
steady reaction, just above the second explosion 
limit. Hence the rate of this reaction must be 
independent of the vessel size, the wall having 
only the effect of destroying each HO» formed 
by VI. This theoretical conclusion is surprising, 
as experimentally the vessel size is known to 
affect the rate of the reaction strongly. 

One may try to explain the discrepancy by 
the fact that rates are generally measured at 
pressures considerably higher than the second 
explosion limit. Such higher pressures obviously 
increase the chance for the continuation of the 
chain by X and so may explain the effect of 
vessel size on the rate. This explanation, how- 
ever, does not agree with the numerical values, 
as is evident from the following comparison. 


18G, H. Grant and C. N. Hinshelwood, Proc. Roy. Soc. 
London A141, 31, Table III (1933). 
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Fic. 4. Rate versus temperature. Quartz vessel clean and 
coated with KCl. 


Grant and Hinshelwood found the second ex- 
plosion limit to be 125 mm (at 552°C for 1:1 
mixture), independent of the size, in silica bulbs 
up to 7.5 cm diameter. At a slightly higher 
pressure, say 130 mm, a steady reaction takes 
place, the rate of which, if the theory is correct, 
must also be independent of the vessel size 
(HO: being 100 percent destroyed at the wall). 
Unfortunately no data are available obtained 
under these conditions, but nevertheless the 
following extension of the argument indicates a 
disagreement between experiment and theory. 
For higher pressures a correspondingly smaller 
vessel will give the HO, the same 100 percent 
chance of being destroyed at the wall, and so 
lead to the prediction of a rate independent of 
the vessel size. (The relation between the two 
cases—low pressure, large vessel and high pres- 
sure, small vessel—results from the idea that the 
probability of an HOz radical being created at 
random in the volume and diffusing to the wall 
depends only on the total number of intervening 
collisions; it is independent of whether at low 
pressures these collisions spread over a long 
distance or at high pressures are concentrated 
within a small distance.) The result is that, for 
example, at a pressure of 300 mm the rate 
ought to be independent of the vessel size for 
bulbs with a diameter of at least 7.5125/300 
=3.1 cm. This conclusion fails to agree with the 
rates given by Hinshelwood and Williamson. 
Silica bulbs of diameters 1.7 and 3.2 cm (identical 
lengths, p= 30 cm, 560°C) showed rates as widely 
different as 0.18 and 0.50. We conclude that the 
effect of the vessel size on the rate is pronounced 


in a range analogous to that in which the corre- 
sponding effect of the vessel size on the second 
explosion limit vanishes. This observed fact is 
contradictory to the theory. On this point, more 
extensive measurements are planned. 


C. Effects of KCl surface 


The effect of a KCI surface on the thermal 
reaction has to be explained by the theory from 
various points of view. 

(1) Effect of (KCl) on the rate.—As discovered 
by Pease," KCI covering the walls reduces the 
rate by a large factor. The experimental results, 
however, show some inconsistency, as Hinshel- 
wood and Williamson’ briefly report that KCl 
in quartz (as distinguished from Pyrex) loses its 
effect of reducing the rate. Lewis and von Elbe 
offer the plausible assumption that the experi- 
ments in quartz were performed at higher tem- 
peratures than those in Pyrex. They argue that 
under this condition the wall loses its effect on 
the rate, since the reaction takes place com- 
pletely in the gas phase. This argument leads in 
particular to the hypothesis that chains start in 
the volume, at least at high temperatures. 
(Chain breaking may be the same in quartz 
and KCl.) 

In order to find the reason for this incon- 
sistency, we performed such experiments in 
Pyrex as well as quartz vessels. In a group of 
experiments with quartz we found a very pro- 
nounced effect of KCl. We varied the tempera- 
ture over as wide a range as is permissible, just 
avoiding the explosion (60 cm pressure). This 
limited the temperature to about 538° in clean 
quartz. The comparison between quartz clean 
and quartz covered with KCI is given in Fig. 4. 
The effect of KCI at a certain temperature would 
be given by the ratio of the ordinates. In another 
quartz vessel, however, the rates were so much 
more erratic that the effect of KCl may be 
almost obscured when comparing extreme values. 
There is an additional effect which may help 
explain the failure of Hinshelwood and William- 
son to find an effect of KCl applied to quartz; 
Prettre'® found that ‘fused quartz changes by 


4 R, N. Pease, J. Am. Chem. Soc. 52, 5106 (1930). 

15 C, N. Hinshelwood and A. T. Williamson, book, p. 50. 
See the book by Lewis and von Elbe, p. 50. 

16 M. Prettre, Mémorial des Poudres 27, 255 (1937). 
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contact with KCI so that the effect of the KCl 
cover rapidly disappears.” But our main ob- 
jection to the argument just mentioned of Lewis 
and von Elbe is still another. The experiments 
of Hinshelwood and Williamson, comparing clean 
quartz and quartz with KCl cover, were cer- 
tainly not performed at higher temperatures 
than those in Pyrex, because in any case the 
third explosion limit is below the softening point 
of Pyrex. For Pyrex and Pyrex with KCl the 
same relation between the curves holds as given 
in Fig. 4, up to the third explosion. 

While these experiments remove one argument 
for the starting of the chains in the gas phase, 
they put the burden on the theory of explaining 
the high effect of KCI on the rate. Such a strong 
reduction of the rate by the wall is not easy to 
explain on the basis of the theory of Lewis and 
von Elbe. In this theory the wall affects the rate 
only by destroying chain carriers. It is assumed 
that the carrier HOe, when colliding with Pyrex, 
is destroyed with nearly 100 percent efficiency 
(see Section II B). Thus, there is no margin left 
for another material to have a still higher de- 
stroying effect on HO». Therefore, we may 
assume that KCl has the power of destroying 
another chain carrier, O, OH, or H. O atoms do 
not have a very long life,!? and therefore have 
not much chance of diffusing to the walls be- 
cause of the rather low heat of activation of III. 

If OH were the carrier with a life so long that 
it might be destroyed at the wall, a very large 
concentration of OH would be built up (in Pyrex 
still larger than in KCl), since each reaction II 
and III would produce one OH which would 
diffuse to the walls. This disagrees with our 
complete failure to detect any OH by its absorp- 
tion spectrum in a layer as deep as 90 cm.!® 
There remains H as the chain carrier that may 
be destroyed at KCI walls. This seems to be an 
exceedingly implausible picture, assuming that 
each link of the chain lives so long that the 
carrier H has an excellent chance to diffuse to 
the surface, where it may be destroyed if the sur- 
face consists of KCl, and may survive if it con- 
sists of glass. (This assumption, which cannot be 
completely excluded, would contradict the re- 

17 P, Harteck and U. Kopsch, Zeits. f. physik. Chemie 


B12, 327 (1931). 
18 QO. Oldenberg and H. S. Sommers, Jr., to be published. 
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mark of Lewis and von Elbe (book, p. 40) that 
‘“‘a chain carrier spends most of its life as HO:.” 
But this statement seems to be unessential to 
their theory.) 

(2) Effect of the surface on the second explosion 
limit.—Lewis and von Elbe base their discussion 
on the idea that only gas phase reactions need 
be considered for the theory of the second 
explosion limit, since it is practically unaffected 
by the walls. This had been observed for Pyrex, 
quartz, and porcelain. It is not confirmed by our 
results on the second explosion limit, comparing 
clean Pyrex with Pyrex coated with KCl (diam- 
eter 2.4’’); coating raised the second explosion 
limit materially. Furthermore, the speed of 
manipulation had some effect. We determined 
such an effect in our KCI covered vessel by the 
following procedure: First, the explosion limit 
was determined to be 18.2 cm. Next, an initial 
mixture of 52 cm pressure was rapidly evacuated 
to 18.4 cm and then let stand. A very slow 
reaction took place. After six minutes we con- 
tinued the pumping until the explosion occurred. 
Although during the waiting period the progress 
of the reaction was slow (only 4 percent of the 
mixture reacting) the second explosion limit was 
lowered from 18 to 12 cm. In general, the longer 
the period of waiting, the lower is the pressure 
of the second explosion limit. The amount of 
the water produced during the waiting period is 
so small that the present theory (attributing to 
H.O only the part of an inert molecule) would 
not explain such an effect. Similar experiments 
performed in a clean Pyrex vessel also indicated 
significant effects of the rate of manipulation. 

(3) Effects of KCl on the rate and the second 
explosion limit.—In the present theory, the effect 
of KCl, strongly reducing the rate, must be 
explained by the great power of KCI to break 
chains by destroying some type of chain carrier. 
This leads to the prediction that KCI should 
make the mixture far less explosive. Instead, 
the opposite is observed. As published in a 
letter to the editor,® KCl makes the explosion 
peninsula expand into the territory which, in 
clean Pyrex (diameter 2.4’’, 540°C), would be 
occupied by the nonexplosive reaction. It seems 
difficult to reconcile this observation with the 
theory. 

One may try to explain the fact that KCl 
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THERMAL REACTION 


strongly reduces the rate by a slow reaction of 
KCl with the glass walls so that the K is 
gradually built into the glass and Cl, is liberated, 
which in the gas phase may cause the effects 
described. That this explanation is wrong follows 
from the fact that the effect of Cl, on the 
second explosion limit is known to be the opposite ; 
it reduces the size of the explosion peninsula. 
Furthermore, the effect of Cle on the rate, 
although it is of the same type as the effect of 
KCI, is very much less pronounced.!® We con- 
clude that the liberation of Cl: is not responsible 
for the effects discussed with KCI covering 
the walls. 

The vapor pressure of KCI is so small that 
such effects cannot be explained by the presence 
of its vapor. At 540°C this pressure is only 
0.00028 mm.” This gives a mean free path of 
such a magnitude that collisions with a vapor 
molecule are not more probable than collisions 
with the wall. Furthermore, in a quiescent gas, 
when only an exceedingly slow reaction is going 
on, one cannot very well expect the vapor 
pressure of KCI to build up to a value far ex- 
ceeding the equilibrium value.”! 


IV. CONCLUSIONS 


The great successes of the present theory make 
it highly desirable to modify it as little as 
possible. Nevertheless, we did not succeed in our 
attempts to adapt it to the experimental results 
discussed in the preceding section. The following 
brief argument is intended only as a preliminary 
account of conclusions which, obviously, must be 
checked by more experimental work. 

The main difficulties are summarized as 
follows: (1) the theory fails to describe the sharp 
kink in the rate-pressure curve; (2) the vessel 
size has no effect on the second explosion limit, 
but a strong effect on the rate in a pressure 
range where, according to the theory, it is not 
expected ; (3) for the surface covered with KCl 
the effect on the rate (slowing down the reaction) 


® Hinshelwood and Williamson, book, p. 76. 

20 Joseph E. Mayer, J. Chem. Phys. 6, 301 (1938). 

*t The situation is different in the experiment on a silver 
vessel, There, according to Lewis and von Elbe, a vapor 
pressure of Ag exceeding the thermal value builds up. 
It cannot be caused by the thermal reaction of 2H:+O2 

ocking out silver atoms, because in a silver vessel this 
reaction is infinitely slow. Presumably it is due to the 
oxidation of silver at 550°C. 


BETWEEN H: AND O; 475 


would theoretically lead to an effect on the 
second explosion limit (reducing the tendency 
toward explosion) which contradicts the facts. 


A. The second explosion limit 


In the search for another theory one may be 
guided by the following consideration which 
would make a surface effect responsible, or 
partly responsible, for the second explosion limit. 
It is known that in the thermal hydrogen-oxygen 
combination two processes are superimposed, a 
gas phase reaction and a surface reaction. The 
gas phase reaction has a very high temperature 
coefficient, the surface reaction a low one.! Con- 
sequently, at high temperatures the gas phase 
reaction is much faster; on the other hand, “‘at 
temperatures below, or not much above 500°C, 
the reaction takes place almost entirely on the 
walls of the containing vessel’ (Hinshelwood and 
Williamson). The second explosion limit has been 
reported by Thompson and Hinshelwood for 
temperatures as low as 450° where the gas 
phase reaction has almost completely vanished. 
Hence it is plausible to attribute this limit to 
the surface reaction. 

According to Garstrang and Hinshelwood,” it 
is true that the surface reaction does not vary 
with the pressure in the required manner, 
showing no sign of increasing as the pressure 
drops from higher values toward the explosion 
limit. But this does not force us to attribute the 
second explosion limit to the volume reaction, as 
this reaction does not show such an effect 
either.” 

Furthermore, it seems difficult to mend the 
present theory—attributing the second explosion 
to a volume effect—by assuming more involved 
steps in the volume reaction. The reason is that 
the treatment of Kassel and Storch seems to 
exhaust the possible volume reactions. 

The main argument of Lewis and von Elbe 
for attributing this limit to a volume reaction 
had been the fact that it seemed practically 
unaffected by the walls. This argument is in- 
validated by the observation that KCI covering 


2 W. L. Garstrang and C. N. Hinshelwood, Proc. Roy. 
Soc. A130, 642 (1931). 

% Alyea’s (J. Am. Chem. Soc. 53, 1324 (1931)) detailed 
scheme for the explosion as originating from the surface 
has been disproved by Garstrang and Hinshelwood (Proc. 
Roy. Soc. A134, 1 (1931)). 





476 


the walls materially changes this limit.2* (With 
KCI the vessel size has a considerable effect on 
the limit, the larger size favoring the explosion. 
This suggests that the ‘‘second explosion’”’ is not 
a straight surface phenomenon, for then it would 
be entirely unaffected by the vessel size. These 
observations will be discussed in a subsequent 
paper.) 

Grant and Hinshelwood* presented some argu- 
ments for the theory that a process in the volume 
is responsible for the ‘‘second explosion.”” The 
present discussion, however, makes us believe 
more in an earlier argument of Thompson 
and Hinshelwood® which leads them to the 
“hypothesis of two simultaneously operating 
mechanisms.” 

The interpretation of the second explosion as 
caused, or partly caused, by the surface would 
straighten out the difficulties mentioned above: 
(1) the sharp kink of the rate-pressure curve 
would be understood by the thesis that, for a 
certain low pressure, a gas phase chain reaction 
is interfered with by some foreign substance 
coming from the surface ; (2) the second explosion 
limit would lose its objectionable correlation with 
the rate as a function of vessel size; (3) the fact 
that KCI slows down the volume reaction would 
not lead any longer to a conclusion regarding 
the second explosion, because this would be 
attributed to the surface. 

It is admitted that by attributing the second 
explosion limit to a surface effect we do not 
explain this explosion as satisfactorily as it 
apparently is done by the present theory; 
instead we shift difficulties to something the 
details of which remain unknown. 


B. The chain reaction 


How would this step modify the present 
theory? The change would be profound—far 


24Q, Oldenberg and H. S. Sommers, Jr., reference 6. 
See also G. H. Grant and C. N. Hinshelwood, reference 13. 
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more so than we should wish. In the present 
theory, the complicated part played by the 
hypothetical radical HO: is largely introduced 
for the very purpose of explaining the second 
explosion limit. Once we attribute this limit to 
another process, there is not much reason left 
for introducing HO: and attributing a consider- 
able chemical stability to it (ability to diffuse 
to the walls) at temperatures as high as 580°. 
Then the Bonhoeffer-Haber reaction may seem 
simpler and, therefore, preferable. In recent 
papers, Prettre*>. investigated the order of the 
thermal reaction and found that it agreed well 
with the Bonhoeffer-Haber scheme. 

One difficulty is caused by the very powerful 
quenching a KCl surface exerts on the rate. 
This observation seems to be hard to explain 
by a gas phase reaction starting in the gas phase, 
so that for the surface no other effect on the chain 
would be left but breaking it. It may be prefer- 
able—in contrast to Prettre’s hypothesis—to 
assume that the chains start at the surface. 
(This will be accessible to an experimental test.) 
Chains starting at the surface have also been 
assumed by Neumark, Kulashina and Poljakow’* 
on the basis of other evidence. The theory of 
Lewis and von Elbe, however, cannot simply be 
modified by assuming that the chains start at 
the surface (instead of reaction XIII). The 
reason is that this modification would not de- 
scribe the facts concerning vessel size and 
rare gases.’ 

Since new experiments are suggested by the 
present discussion, it seems futile to elaborate 
the theory any further before the experiments 
are carried out. The experimental work is being 
continued. 


25 M. Prettre, J. Chim. Phys. 33, 189 (1936). 

26 E. Y. Neumark, L. P. Kulashina and M. V. Poljakow, 
Acta Physicochim, U. R. S. S. TIX, 733 (1938). 

27 Lewis and von Elbe, book, p. 68. 





present 
by the 
oduced 
second 
imit to 
on left 
nsider- 
diffuse 
s 580°. 
y seem 
recent 
of the 
ed well 


owerful 
e rate. 
explain 
; phase, 
e chain 
prefer- 
‘sis—to 
surface. 
il test.) 
o been 
jakow"® 
eory of 
iply be 
tart at 
). The 
10t de- 
re and 


by the 
iborate 
‘iments 
s being 


‘oljakow, 


JUNE, 1940 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 8 


Theoretical Isotherms for Adsorption on a Square Lattice 
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A method which is not self-consistent has been used by 
J. K. Roberts to calculate the adsorption isotherm of 
adatoms, which are larger than the spacing of the adsorp- 
tion sites on the substrate so that each adatom excludes all 
others from neighboring sites. It consists in completely 
formulating the internal statistics of a sample of the film 
surface and allowing for the effect of the rest of the film 
by introducing one or more ‘‘environmental” parameters 
as undetermined functions. These are determined by the 
requirement that each site of the sample shall be an average 
site. The failure of this method is traceable to the fact that 
these functions, when determined, are found to vary with 
surface concentration in a manner inconsistent with the 
physical requirements. The same difficulty is found to be 
latent in Peierls’ treatment of adsorption with interaction 
between adatoms and in Bethe’s treatment of super- 
lattices. The adsorption isotherms for two cases with no 
interaction between neighbor adatoms are here calculated 
by another method. The one case is that of Roberts, the 
other that in which 9 sites are excluded by each adatom. 
The method requires separate treatment of dilute (small @) 
and concentrated (@ nearly unity) films. At each extreme 


it involves successive approximations. For the dilute film’ 
the number of free sites is determined as a function of the 
number of adatoms by enumerating the number of ways 
that adatoms occur (1) without overlapping of exclusion 
patterns, (2) with the exclusion patterns of only two 
adatoms overlapping, etc: This yields the adsorption 
isotherm. For the concentrated film, the interactions of 
the free sites themselves are treated in an analogous but 
necessarily somewhat different manner to give the ad- 
sorption isotherm near 6=1. Some fairing of the relations 
to bridge the intermediate values of @ is necessary. The 
results are: for the 5-site exclusion pattern 

6=2A pli — (5/2)0+(3/2)0+3/467], 0050.5 

6=Ap[1—0—3(1—6)?+5(1—6)*], 0.55051 
and for the 9-site exclusion pattern 

6=4A p[1—(9/4)0+6?], @ near zero 
6=2Ap(1—6)*, 6 near unity. 

These cases are compared with the limiting cases, first 
the 1-site exclusion pattern and second the infinity-site 
pattern which is the same case as adsorption on a siteless 
surface. 





1. INTRODUCTION 


I‘ conjunction with more general investiga- 
tions by Langmuir! on the behavior of 
monatomic films of adatoms on an underlying 
crystal lattice, I undertook an investigation of 
the statistics of simple cases in which the 
adatoms were located on sites on the substrate 
and were large enough to prevent adsorption on 
neighboring sites. This was in 1936. 

Since then several similar investigations have 
appeared?* which found their inspiration in 
Bethe’s treatment of the order-disorder problem‘ 


! J. Chem. Soc. (London), April or May, 1940. 

*R. Peierls, Proc. Camb. Phil. Soc. 32, 471 (1936) 
first applied the general method to surface films. He was 
concerned with the case where the density of the close- 
packed film is equal to the density of sites, and the inter- 
action energy between neighbors is negative. 

*J. S. Wang, Proc. Roy. Soc. A161, 127 (1937) and 
J. K. Roberts, Proc. Camb. Phil. Soc. 34, 577 (1938). 
J. S. W. treats the case where the interaction energy is 
positive, and hence can lead to exclusion—see curve a of 
his Fig. 2. J. K. R. treats exclusion plus interaction 
between permitted neighbors. The method used by these 
Investigators will be designated as the PWR method. 

*H. A. Bethe, Proc. Roy. Soc. A150, 552 (1935). 


and are based on the use of a partition function. 
Such treatments have the appearance of deriving 
the adsorption isotherm of a surface film to a 
good approximation over the whole range form 
6=0 to @=1. On the other hand, I found it 
necessary, with my method, to treat the sparsely 
populated and the densely populated films as 
two separate cases, and to fill an intervening 
gap by “‘intelligent’’ guessing. The calculation 
became rapidly more complicated as_ higher 
orders of approximation were sought, and, I fear, 
leads at best to semiconvergent series. 

But the difference between my results and 
those achieved by the newer method has led to 
a critical examination of the situation which 
exposes fundamental inconsistencies in the latter 
when it is applied to a case where the number 
of sites exceeds the number of close-packed 
adatoms. 

In the following, I shall first analyze the 
newer method so as to exhibit its paradoxes, and 
second describe and apply the older more 
rigorous method. 
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Fic. 1. Site arrangement of Roberts. Tonks’ sample of 
film is outlined by the dashed lines. 


2. THE PEIERLS-WANG-ROBERTS 
(PWR) METHOD 


This method consists in choosing a small 
section of the adsorbing surface as a sample 
and writing down its partition function, taking 
account of the effect of the unanalyzed environ- 
ment by undetermined functions as parameters. 
These are fixed by the condition that every site 
of the sample shall be an average site, that is, 
that the probability that it is occupied by an 
adatom is 0, the average coverage. 

The model in which Roberts* was interested 
involves an adatom whose exclusion pattern 
covers the adjacent sites. Thus an adatom at 
site 0 on the 110 crystal face shown in Fig. 1 
would prevent adsorption on 1, 2, 3 or 4, that is, 
it has an exclusion pattern which covers these 4 
sites. In addition, adatoms forming a pair, as, 
for instance, two on 0 and 5 or 0 and 6 or 2 and 
4, but not 0 and 7 or 0 and 8 (because they are 
further apart) have an interaction energy V. 
The adsorbed layer is in equilibrium with the 
gas, the adsorption energy per atom being x. 

As the sample area the seven sites, 0—6, were 
chosen. The partition function consists of a sum 
of terms of the type 


Zn Ci Conzgnotnitns (1) 


one for each of the possible configurations of the 
sample area. ~ is identified with the partition 
function B of the gas* per molecule, the adsorp- 
tion energy x, and the partition function® v(T) 
for the adsorbed states 


QanmkT\i kT 
t=0(T)ex!*"/B, 5-(——) ~-4,(T); (2) 
h? p 


5 R. H. Fowler, Proc. Camb. Phil. Soc. 31, 260 (1935). 
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Mo, M1, M2 are, respectively, the number of 
adatoms in the central site (0 or 1), the first 
ring consisting of sites 1-4, or the second ring 
consisting of sites 5 and 6; 


n=evlkT 


and p is the number of interacting pairs in the 
configuration. {; and ¢2 are the environmental 
parameters which can be looked upon as 
corresponding, respectively, to energies V; or V2 
of interaction between ring 1 or ring 2 and the 
surrounding film. g is the multiplicity of the 
configuration. 

In order to maintain uniformity with Lang- 
muir’s notation. @z will be used for the PWR 4, 
and @ will here denote the fraction relative to 
close-packing that the surface is covered. It is 
apparent that 

0=206r. (3) 


By formulating successively the probability 
of occupancy of the central site, a ring 1 site, 
or a ring 2 site in terms of the probabilities in 
the partition function, and then setting them 
each equal to 6z, three equations result. From 
these ¢; and {£2 can be eliminated, leaving a 
relation between @ and £é which, in effect, 
constitutes the adsorption isotherm by virtue of 
Eq. (2). Roberts has plotted such isotherms up 
to £=400 for »=0.01, 0.1, 1.0 to cover a range 
of interaction energies. 
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Fic. 2. Adsorption isotherm characteristics calculated 
by Roberts’ method with zero pair interaction energy. 
R is for Roberts’ sample of film. T is for Tonks’ sample 
of film. 
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Henceforward we confine our attention to the 
case where the pair interaction energy is zero. 
Then »=1 and by Roberts’ Eqs. (9) and (2), 


f1=f2(=$). 


At high values of ~ (high gas pressure, low 

temperature) it is easily found that asymp- 
totically : 

c=F4, Or=2; 6=1 

1 

g=1-—6=1-26rg=— =f 

, &f 


whence, in the limit as 6-1, 
0/(2§)—¢"/2. (4) 


By a general line of reasoning it can be shown 
that for isothermal changes in the adsorbed film® 


6=(2Ap)6;, (5) 


where A is a constant, p the gas pressure, and 
6; the fraction of the total number of sites which 
are free to adsorb. When no sites are covered by 
adatoms on adjacent sites, it is evident that 
6;=1-—6 and we have Langmuir’s adsorption 
isotherm. From Eq. (5) it thus appears that in 
general 0/(2Ap) is the fraction of ‘‘free’”’ sites, 
and A is fixed by the requirement that at 
6=0, 6/(2Ap) shall be unity. 

Now in Eq. (4), & is proportional to p and 
6/2 is unity at 6=0 as can be verified from 
Roberts’ equations. Accordingly 6/2¢ is the 
§/(2Ap) of the previous paragraph. Thus Eq. (4) 
becomes 

0;= ¢°/2. (6) 


This equation appears as ‘‘R asymptote”’ in 
Fig. 2. Hence the analysis yields the result that 
for an almost complete film the relative number 
of free sites is proportional to the square of the 
film’s incompleteness. 


TABLE I. Component partition functions of sample area. 














nm Pn, Pn, 
0 (E+1)(y+1)*=Po 
1 4x(y+1)? =P 
2 (adj.) 4x*(y+1) =Poe 
2 (opp.) 2x? =L 20 
3 4x3 =P, 
4 x =P, 











® See Eq. (41) of reference 1. 
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TABLE II. Calculation of isotherm according to Roberts 
using Tonks’ sample. 








y 0.1 1.0 3.0 10 30 





x 0.1089 1.500 5.355 21.30 77.4 





g 0.1189 2.033 6.357 16.92 40.86 











26r | 0.155 0.572 0.78 0.92 0.98 
v1 0.91 0.74 0.84 1.26 1.90 
$2 0.84 0.49 0.47 0.59 0.73 























Once pair interaction energy has been aban- 
doned, however, the choice of the 7-site sample 
appears somewhat arbitrary. The sites 0, 1, 2, 3, 
4, 7, 8 are in every way equivalent to 0, 1, 2, 3, 
4, 5, 6 and one should reasonably expect a more 
exact result if all the sites 0-8 were included in 
the sample.* This result should not differ signifi- 
cantly from Eq. (6). 

The following procedure to handle this case, is, 
I think, somewhat less unwieldy than Roberts’. 
With him the “environmental” energies V; and 
V2 are combined with the “‘site’’ energy to give 
the parameters 


xaeVtVis/kTB-1 yaeVtVaik? B-1, 

The x and y are his «; and es. The partition 
function for the sample can now be formulated 
from a table each row of which contains the 
partition function for one possible configuration 
of the first ring (Table I). For the probabilities 
of finding the central site occupied, or a given 
first ring site occupied, or a given second ring site 
occupied, we then have: 


£0(EPn)/dE (x/4)0(EPny)/dx 

. i. =P 
(y/4)0(2Pn1)/dy 
r EPn1 , 





R 





The last two equalities lead to two relations, 


*Note added in proof.—This analysis has recently been 
made by Roberts in an article (Proc. Camb. Phil. Soc. 36 
(1), 53 (1940)) which reached the author after the present 
paper was submitted for publication. [n that article the 
nature of the failure of the statistical method is not 
analyzed, but Roberts resorts to an empirical method to 
handle values of 0g greater than 0.25. He does not seem to 
recognize that this technique does not lead to true equilib- 
rium conditions (reference 1). 
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Fic. 3. Behavior of environmental parameters ¢, and £2: 
(1) Qualitatively according to physical considerations, 
curves marked ‘“‘odd” and “‘even.”’ (2) As given by Roberts’ 
method, curves marked T. 


one between x and y: 
x(1+x)?—y(1+y)*+xy(x —3y—xy— 2xy*) =0, 
and the other for é in terms of x and y: 
é=y[1+x(2+x+2y)(1+y)-*]. 
Then from the first equality: 


Or m &(1+~y)4 
1-On (1+y)*+(1+x)*—14+4xy(2+a+y) 





Choosing a set of values of y, these equations 
were solved in order for x, £ and 6p, and the result 
is tabulated in Table II, and represented in 
Fig. 2 as Curve T. For the larger values of x it 
was necessary to use an expansion valid for this 
region which, in the limit, gives 


O;=¢. 


Comparing this with Eq. (6) we are thus 
presented with the anomaly that an additional 
step toward closer approximation leads to a 
highly divergent result! 

The nature of the trouble can be seen from the 
course of ¢; and {2 which are also tabulated in 
Table II. The only influence of the environment 
on the sample lies in the possibility that the 
exclusion pattern of an external adatom covers a 


sample site. Accordingly this influence is always 
unfavorable to ring 1 and ring 2 occupancy. 

The behavior of Bethe’s € is a guide to the 
reasonable behavior of the ¢’s. He found that « 
(comparable with ¢) was double-valued. The 
value lying between unity and infinity corre- 
sponds to a predominantly “‘right’”’ environment, 
the value between unity and zero corresponds to 
a predominantly “wrong” environment. The 
situation here is similar. A complete film will 
occupy either the 0, 5, 6, 7, 8, 9, 10, 11, 12 etc. 
sites or it will occupy the 1, 2, 3, 4, 13, 14, 15, 
16, 17, 18, etc. sites. In Langmuir’s terminology 
it will either be an even film, which we can 
arbitrarily identify with the former, or an odd 
film. Thus in place of “‘right”” and ‘“‘wrong’”’ we 
have ‘“‘odd” and ‘‘even.” It is apparent that ¢, 
for an odd film decreases to zero with increasing 
6. On the other hand, fe finally increases to 
unity because the odd environment excludes ring 
2 sites less and less. Accordingly we should 
expect the situation shown qualitatively by the 
f10aa and {£2 oaa Curves of Fig. 3. For an even 
film the situation would be approximately re- 
versed, as shown by the £1 even and £2 even Curves 
in the same figure. Therefore, any solution which 
is valid for the whole range of @ must have two 
pairs of values for the ¢’s, which behave in 
accordance with the foregoing analysis. It follows 
that the equality of ¢; and £2 found by Roberts 
proves the incorrectness of the method, as does 
the variation which the newer calculation gives 
for the ¢’s. These values are plotted in Fig. 3. 

It seemed possible that if a large enough 
sample were taken so that the environment did 
not directly affect sites whose probability of 
occupancy was set equal to @, then a valid result 
might be obtained. I tried this with a square 
sample of 13 sites (0-12 in Fig. 1). The equations 
were 


6= Prob. for site 0= Prob. for site 1. 


This meant that only one parameter could be 
introduced. Noting that site 5 could be covered 
from the outside in two ways, whereas 9 could 
be covered in three, I chose 


Scorner _ (+ stiend. 


This, I think now, was not a good guess, but the 
attempt has convinced me that the method is 
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unworkable, for the nature of the equations is 
such that with finite ¢’s, & cannot approach 
infinity. It remains possible that by replacing ¢ 
with some function of — or 6 the method might 
be salvaged. 

The environmental parameter does not lead to 
these difficulties in the range where Bethe’s' 
results are of interest; namely, in the neighbor- 
hood of the transition to complete disorder. His 
two successive approximations do not give 
diverging results, and, as has been discussed, the 
variations of e’ and 1/e’, the two possible values 
of €, are consistent with the physical require- 
ments. But as complete order is approached ¢’ 
goes to infinity. However, the most favorable 
environmental effect will occur when all the 
environmental neighbors, approximately Z/2 in 
number, are right. This means that the largest 
physically attainable value of ¢’ is of the order 
«4/2, A similar argument applies to 1/e’. This is 
a limitation of the theory where the ordering 
tendency is high—just the range where the 
adsorbed film behavior is being analyzed. 

Peierls? has dealt with an adsorbed film 
without encountering the present difficulties. 
He considers adatoms whose exclusion pattern 
includes only the site which the adatom occupies. 
It is interesting to note what happens when an 
exclusion pattern covering nearest neighbors, 
which is latent in Peierl’s treatment, is accentu- 
ated. This can be brought about by putting his 
interaction energy, — V, equal to +. Then 7, 
his Eq. (4), becomes zero, and solving 


e= FE 
and his Eq. (10) with »=0: 
t=e(1+6)2-! 


simultaneously for ¢ for large we have 
c='-(2-D/z_,9, 


This is absurd, because, by assumption the 
central sample site is an average site so that a 
valid method would yield a second value of ¢ 
which approaches unity. 

Wang has brought out the latent difficulty in 
P’s equations by applying them to the case 
where 7 is small. So doing simply gives another 
in the set of divergent results. A replot of Curve a 
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of W’s Fig. 2 (n=5.53X10~') above £=20 onto 
the coordinates of the present Fig. 2 gives the 
curve W. 

Faced with these failures to solve the problem 
before us by Bethe’s ingenious stratagem, it is 
pertinent to examine the older method. 


3. CALCULATION OF ADSORPTION ISOTHERMS— 
ASSUMPTIONS AND NOMENCLATURE 


The cases which will be discussed are some of 
those listed by Langmuir.' A square array of 
sites on a substrate is assumed, each one of 
which can hold an adatom with a certain energy 
of attachment. The adatoms are supposed to 
have a diameter d,, the sites to be a distance d, 
apart. Thus depending on the ratio d,/d,, one 
adatom will prevent any other atom from 
occupying a total of 1, 5, 9, 13, 21 etc. sites. 
Any such array of sites constitutes an exclusion 
pattern and any site in the array is excluded. 
The central site is occupied and the rest of the 
excluded sites are covered. Any site which is not 
excluded is free. 

There is the further assumption that there is 
no interaction energy between adatoms (save 
that which maintains the exclusion pattern). 

We now define: E=number of sites in the 
exclusion pattern, ¢,=number of sites per unit 
area, og =number of excluded sites per unit area, 
o;=number of free sites per unit area, o,=num- 
ber of adatoms per unit area when close-packed, 


7 =Cs; ‘o1, (7) 
o=number of adatoms per unit area, 
O;=a;/0s, (8) 


and, as before, 
6=0/04. (9) 


The method to be used depends upon Eq. (5), 
used here in the form 


The procedure will be to calculate 6; in terms of 
6, or what amounts to the same thing, oy, in 
terms of o. 

Incidentally, from the previous discussion we 
see that 


A=t/p 
with & given by Eq. (2). 


(10.5) 





482 


4. ADSORPTION ISOTHERM FOR E=5, y=2 FOR 
Low CONCENTRATION 


This is the case the treatment of which by 
Roberts has already been discussed. Each adatom 
excludes its own and the four neighboring sites. 
The lattice need not be square but can, within a 
certain range of side ratio, be face-centered 
rectangular. If the sides of the elementary 
rectangle are d; and d2, di<d2, the range of 
adatom diameter within which the present case 
applies is 


3 (d?>+d2?)!<da<di. 


The limiting proportion for the rectangle occurs 
when the range for d, vanishes, i.e., when 


d2/d,=3}. 
For the square array, 
d,=d,=2'd,, 
so that the range of d, is 
dy<da<2'd,. 


For convenience we shall think invariably in 
terms of a unit area of surface. 

We shall be repeatedly interested in the proba- 
bility that a given site is occupied. Evidently 
this is o times the chance that it is occupied by a 
particular adatom, which is 1/¢,. This chance has 
been lessened by the fact that there is some 
probability that the particular site lies within an 
exclusion pattern. If it is known that this is not 
the case, that is, that the site is not excluded by 
any of the o—1 other adatoms, then the proba- 
bility of finding the particular adatom on it is 
1/o;. The relation between the two at low 
concentrations, i.e., when oo, is apparent from 
the fact that each adatom excludes 5 sites so 
that 
(11) 


This cannot be the general relation for o,, for it 
assumes that every adatom excludes 5 sites, none 
of which are excluded by any other adatom. 
Whenever two (or more) exclusion patterns 
overlap oa; is greater than the Eq. (11) value. 
Thus the problem of calculating o; is that of 
enumerating the ways in which exclusion pat- 
terns can occur, singly and overlapping, and the 
frequency of each configuration so that the 


os=0,—Se. 
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excluded areas contributed by all of them may 
be summed for subtraction from @s. 

The calculation can be made by successive 
approximations. The first is given by Eq. (11). 
The smaller o the better an approximation it is, 
for the smaller relatively is the number of over- 
lapping patterns which it neglects. The next 
approximation considers the doublets, that is, the 
pairs of adatoms whose exclusion patterns over- 
lap; and the third takes triplets, three adatoms 
which have a continuous exclusion pattern, into 
account. Since each doublet contains two singlets, 
and each triplet contains 3 singlets and two or 
three doublets, there is a possibility of over- 
counting. To avoid this, only those adatoms can 
be counted as singlets, whose exclusion patterns 
are isolated from other exclusion patterns, so 
that it is known that they are not parts of 
doublets, triplets, etc. These will be called 
isolated singlets. Similarly, only isolated doublets 
will be counted. 

In this connection it becomes necessary to 
calculate the probability of finding an adatom at 
a site when something is known about neighbor- 
ing sites. For instance, given an adatom at A, 
Fig. 4, what is the probability of finding one 
at B? It is greater than o/o, because B certainly 
is not covered by an adatom located at either 
b or d, and it is less than o/o; because it may be 
covered by an adatom at L or K. 

Consider the probability that a particular 
atom, X, is located at B. It is 1/o,. But alter- 
natively it is the probability that B, b, d, Land K 
are unoccupied by any other atom and that X is 
on the otherwise free site A, namely, 


1 o\?> 
—{1-——)}. 
os Os 
Since every adatom can have the role of X, 
we have then two expressions, necessarily equal, 
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for the probability of finding a site occupied: 


oc o@ a\§ 
—=— 1-=) (12) 


Os Of Os 


This equation is identical with Eq. (11) to the 
first-order approximation, which is as far as 
either of them are valid. Formulated in terms 
of X, the right member of Eq. (12) appears 
reasonable enough, but considered in terms of 
the probability of finding an adatom in general 
at B, it appears somewhat paradoxical, for it is 
the product of the six probabilities that b, d, L, K 
and B are unoccupied and that B as a free site 
is occupied. 

The original question can now be answered. 
The probability of finding an adatom at B when 
one is already at A is the probability that L, K 
and B are unoccupied and that B is occupied as 
a free site, that is, 


(1 ~*) efor 


Since the approximation we shall make is to 
take triplets into account, the relative number 
of which is of the order o*, all probabilities must 
be correct to terms of this order. 

We first consider isolated singlets, an adatom 
at A is isolated if B, C, D, E, F, G, H, J, are all 
unoccupied. On this basis the number of isolated 
singlets would be 


(-)-(-) 


This is valid only to o*, thaf is, as far as it has 
been expanded. To justify another term in the 
expansion the o/o, needs correction. The rule 
already established is applied to occupancy of B. 
The probability of an adatom at B is 


(1—0/0,)°a/o; 


rather than o/c, since we know that 6 and d are 
unoccupied. For C, since of its neighbors only d 
is certainly unoccupied, the probability is 


(1—a/0,)*a/o;. 
The probability that B and C and the other 6 
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sites are all unoccupied is then 
[1—(0/0s)(1—0/0,)* [1 — (0/04) (1—0/o,)*}*. 


Since the final approximation does not distin- 
guish between o, and o;, os can be substituted 
for o,. Putting 


¥=0/0;=0/y6; (13) 


for convenience and noting that there are o 
adatoms like A, the total number of isolated 
singlets is seen to be 


o[1—y¥(1—y)* #L1—y(1—y)*}', 


each excluding 5 sites. : 

There are two types of doublet, the d type 
represented by adatoms at A and B, and the s 
type by adatoms at A and J. 

In the former case, the occurrence of an 
adatom at B, as we have already seen, is governed 
by the probability ¥(1—y)*. For this doublet to 
be isolated requires further the unoccupancy of 
sites V, C, D, E, F,G, H, J, Q, M, ten in number, 
introducing the probability factor (1—y)!®. 
Other d doublets involving A are AD, AF, 
and AH. The last two will be counted with sites 
F and H, respectively, the first requires a 
multiplicity factor 2 in the probability. Thus the 
number of d doublets is 


2oy(1—y)" each excluding 8 sites. 
Similar reasoning for s doublets gives 
2o¥(1—y)!* each excluding 9 sites. 


Of triplets there are 5 types, all shown in 
Fig. 5. Their common probability factor is y*. 
This is to the required degree of approximation 
as it stands. The number of type a with mul- 
tiplicity 2 is 2cy” excluding 13 sites each. 

The information for all five types is given in 
Table III. From this it appears that the whole 
number of cells excluded by triplets is 


240cy”. 
r-2d5+ 

b 
TERR 


Fic. 5. Triplets—arrangements of three atoms having 
overlapping exclusion patterns. 
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TABLE III. Types of triplets and their characteristics. 

















NUMBER OF TOTAL 

EXCLUDED EXCLUSION 

TYPE SITEs MULTIPLICITY FOR TYPE 
a 13 2 26 
b 12 8 96 
c 13 4 §2 
d 11 4 44 
e 11 2 22 

Total Multiplicity 20 

Total Exclusion, all Types 240 











By way of a check on the method used here 
to calculate the occurrences of the various con- 
figurations, we note that the sum of the isolated 
singlets plus twice the isolated doublets plus 
thrice the triplets should equal o identically 
through terms in y*. It does. 

The sum of the excluded sites from (5), (6), 
(7) and (8) plus the free sites is equal to the 
total number o;: 


o,=0;+5o[1—y(1—y)* #L1—-v(1—¥)* }* 
+1609(1—y)¥+ 180y(1—y)"©4+240cy%. (14) 


Dividing Eq. (14) through by o;, expanding 
the right member through y’, and noting Eqs. 
(7), (8) and (9), we find 


05-*=1+58/y8;) — 6[8/78s) P+24[0/76;) P. (15) 


Formally, the problem is now solved, for this 
gives the required functional relation between 0 
and 6; which, with Eq. (10) establishes the rela- 
tion between # and @. 

To express 6; in terms of 8 we assume the 
result with undetermined coefficients: 


6;=1+00+86?+56° 


and fix them by substitution in Eq. (15). As a 


result 
6;=1—30+$30?+ 36°. (16) 


Thus for the adsorption isotherm we have 
6=2Ap[1—$0+$6+30°]. (17) 


Curve A; in Fig. 6 is a plot of 6; against 1—0@ 
down to 1—6@=0.5. 

The whole treatment is dependent upon treat- 
ing ¥ as a small quantity. Eq. (14) is seen to 
contain powers of (1—y) up to the 16th. In 
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expanding to 
1—16y+120y?—560y’ 


all terms but the first two were necessarily 
neglected. Rigorously this limits the present 
result to values of ¥ of the order of #5 or less, so 
that @ is limited to less than 7g by Eq. (13). 
This is a severe restriction inherent in the 
method which would make further approxima- 
tions along the present line futile. The fact that 
the coefficients 6 and 24 in Eq. (15) are differ- 
ences between much larger numbers suggests, 
however, that this limitation may be more 
severe than necessary. Besides, the succession of 
coefficients in Eq. (17) points even more strongly 
in the same direction. It is for this reason that 
the curve has been carried down to 6;=0.22 
at 6=0.5. 


5. ADSORPTION ISOTHERM FOR E=5, y=2 
HiGH CONCENTRATION 


In order that full coverage of the field should 
occur, that is, that c=o,=0,/2, it is necessary 
that all adatoms lie in one of the two possible 
but identical lattices. Langmuir! discusses how 
the coexistence of both leads to long time effects 
in reaching a true equilibrium. In the present 
analysis I confine my attention to a nearly 
complete film, which requires that it be pre- 
dominantly of one type, let us say the even one. 

At first every adatom which evaporates from 
a complete lattice leaves only one free site 
behind it. But when four adjacent adatoms in a 
square have evaporated, case 4a of Fig. 8, it is 
seen that there are five free sites; for the central 
site is uncovered when the last of the four corner 
sites loses its adatom. It is thus possible that o; 
exceeds oi:—¢. On the other hand, an atom 
condensing into the central site, case 3a of Fig. 8 
(forming a single adatom in the odd lattice) 
excludes four other previously free sites. Thus 
a; can be less than o:—c. In general this will be 
the case because there is a strong tendency to 
form the odd lattice as o decreases. 

We shall denote the deficiency in o; by oa thus: 


Ca=01—o—o; (18) 


and shall then calculate o4 in terms of the other 
quantities in Eq. (18). Finally, by eliminating a, 
a; can be expressed as a function of (1— 8). 
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This problem could probably be approached in 
the same way as the dilute film except that we 
consider free: sites rather than adatoms as dis- 
tributed over the surface. There is, however, the 
essential difference that the adatoms are con- 
served, whereas the number of free sites depends 
upon the configuration of the film. 6 fixes o, but 
the same @ can yield many values of o;. 

We shall again in this case deal in isolated 
configurations in order to avoid overlapping and 
multiple counting. To assure isolation it is 
necessary and sufficient that the configuration 
shall be bounded by a single close-packed line of 
even adatoms. See Fig. 8. The omission of a 
corner adatom would permit another configura- 
tion to have a covered site in common, the 
omission of a side boundary atom would allow a 
common free site. As before we can use the terms 
singlets, doublets, etc. but now they apply to 
free sites. 

Since the equilibrium condition of the surface 
is independent of the mechanism of evaporation 
and condensation, provided they are not incon- 
sistent, the analysis can be based on the assump- 
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Fic. 6. Various adsorption isotherm characteristics. 
A,, y=2, Eq. (17), 0<0.5. As, y=2, Eq. (29), @>0.5. 
B,, Bz, larger adatoms, y =4, Eqs. (32) and (31). C. Smaller 
adatoms, y=1. Eq. (33). U. Adatoms on a uniform field. 





tion that adatoms evaporate at random from 
the occupied sites per element of time, and that 
the same number condense at random on free 
sites. The probability of evaporation P, of a 
given adatom is, therefore, 


P.=n/o=no;/0 (19) 


and the probability P,. that condensation occurs 
on a given free site is 


P.=n/oy;. (20) 


Consider a sample area of the surface con- 
taining a certain number of sites which are 
grouped around site S, and denote the configura- 
tion of the adatoms as type a. Let the total 
number of @ configurations which exist on the 
surface in equilibrium be o.. We suppose that 
the evaporation of any one of p, adatoms from 
a will give a new configuration 8, whose total 
occurrence is og, and likewise that the condensa- 
tion of an adatom into any one of pg free sites in 
6 will give an a configuration. 

Accordingly, the rate a—£ will be 


P pala 
and the rate B—a will be 


P .pgog. 


The principle of detailed balancing requires that 
these rates be equal, whence, using Eqs. (19) 
and (20), 


0 8=a(Pa/ps)0-'(o;/0;). (21) 


This relation enables us to calculate successively 
more complicated configurations. It is convenient 
to put 


o;/01=20;= oy. (22) 


Let the a and £6 configurations be those of 
Fig. 7A. The multiplicity of @ is unity. The 
probability that the 9 lattices points shall be 
simultaneously filled is @°, and the number of 
sites which could be S is o;. Accordingly, 


Co=1XPXo,=068. 
Also, 
Pa=1, ps=1. 


Therefore by Eqs. (21) and (22) 


op= 8 a;=0108 gy 
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and this is the number of isolated free sites, to be 
denoted by on, 


os1= 0108 yy. (23) 


For isolated doublets, d type, we have the a 
and £6 configurations of Fig. 7B. Evidently o. is 
directly related to o,; with a multiplicity of 4: 


a= 40 710° =40:0"9,;. 
Also, 
Pa=1, pg=2. 
Therefore 


CO s2a= 2010 y,?. (24) 


Likewise, for isolated doublets, s type, the 
a and 6 configurations are shown by Fig. 7c and 


O 52s = 2010'79,;?. (25) 


The general rule for the occurrence of these 
configurations is that 


o=Ao 0" 9;7'-=, (26) 


where A is the multiplicity, 2 the number of 
adatoms in the configuration, and 2, the number 
of even sites. 

The validity of this method of calculating has 
been verified for first, second and third orders 
by adding the number of free sites due to iso- 
lated singlets, isolated doublets and triplets using 
Eqs. (23), (24) and (25) and 5 similar expressions 
for the triplets, and checking that the sum is 


@ 


involving isolated free sites. 
@: -@ B Dots show sites, small circles 

show free sites, large circles show 
adatoms. 
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equal to co; through terms in o,°. Beyond this the 
configurations for quadruplets etc. are so many 
in number as to make wider verification un- 
practical. That already carried out is sufficient 
to demonstrate the soundness of the counting 
method. 

Accordingly, from here on we confine ourselves 
to configurations which contribute to og. The 
first that we consider is that of 4a of Fig. 8 
where 4 adjacent free sites in a square yield a 
fifth free site in the odd lattice. That the rule 
already stated in Eq. (26) applies in this and 
future cases will be assumed; in any case it can 
be confirmed by the method of Eq. (21). 

Using a subscript for. o from the representa- 
tions of Fig. 8, 


4g = 010!" 054. 
For each occurrence 4a contributes an amount 
Dato op: 
2a=21-—2-—2,, 


~, being the number of free sites. This expres- 
sion for the individual configuration is strictly 
analogous to Eq. (18) for the film as a whole. 
Since 
2;=5 and 2,-—2=4, 


it follows that 
Za= —1 for 4a. 


TABLE IV. Configurations in incomplete film. 
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MULTIPLICITY SITES 
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No. oF 
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Fic. 8. Configurations contributing to deficiency of free sites to order of approximation ¢/. 
The cross-hatching distinguishes the adatoms in the odd from those in the even film. 


In Fig. 8 all of the configurations for which 
~i—2 is 3 or 4 and which contribute to 4 are 
shown and the information regarding them is 
tabulated in Table IV. 

The successive approximations consist in 
counting configurations whose o’s contain suc- 
cessively higher powers of g;. At this point all 
those for which its exponent, 2,—2, is less than 
5 have been counted. Breaking off the process 
here, and adding the contributions to og we have 


oa=oi¢s*[ 3013+ (— 02+ 12617 + 249154 89!) v, J. 
From Eq. (18) 
g4=01(1—0— gy) 


and oq can be eliminated. We now replace 6 by 
1—y, g small. The bracketed expression is valid 
to the order ¢ or gy of small quantities. Accord- 
ingly 6 is expanded to 1—13¢ whereas 6'”, 61, 
#° and 6'’, being multiplied by ¢g;, are each re- 
placed by unity. We then have 


gs=9—¢7*[3(1—13¢)+439;]. (27) 


It is worth noting that this equation will give 
the relation between ¢ and ¢; more accurately to 
higher values of g than would the series expan- 
sion for gs in terms of yg. This is because the 
exactness of the ¢,* coefficient of the bracket 
would be lost in the transformation. Eq. (27) 
solved for ¢ is 





1+39;2+439,3 
14399 


ee | 


Table V gives a few pairs of values of ¢ and ¢y. 
Values of 0; are also included. This checks well 
with Eq. (16) at ¢=0.5. 

A series for 0; in terms of g which gives even 
better agreement with Eq. (16) at g=0.5 and 


reproduces Table V in the range below g=0.5 is 
0;=3(9—3¢°+5¢%). (28) 


Some values of 6; calculated by this formula are 
listed in the table under the symbol 6723. 
Thus the isotherm 


6=A,[1—0—3(1—0)®+5(1—8)*] (29) 


can be adopted for the range 
1=6=0.5. 


This is plotted as curve Az in Fig. 6. 

To carry the approximation another step 
would be difficult, although not impossible, for 
30 different configurations would have to be 
taken into consideration. 


6. ISOTHERMS OF ADATOMS OF OTHER SIZES 


For larger adatoms the same methods that 
were used here for low values of @ can again be 
applied, but near 6=1 new complications arise. 
For instance, in the next range of atom size, 
2!d,<d,<2d,, four different close-packed lattices 
(instead of the previous two) become possible; 
and ten pseudo-close-packed lattices for which 
6=0.8 can be found. All the adatoms in the 
latter are related to each other by the knight’s 
move in chess. As the atom size increases, such 
alternatives multiply; so that the range of ¢ 
becomes ever smaller within which a fairly 
simple analysis can be applied which is based on 


TABLE V. 
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Fic. 9. Close-packed film of adatoms (crosses) excluding 
9 cells each. Solid squares show exclusion pattern, vertical 
crosses show alternative filling of sites. 


the single ultimate lattice that must exist for 0 
to be unity. 

There is probably still some significance to a 
low ¢ analysis for this case, particularly as it 
has certain features not found with the smaller 
adatoms. 

A regular close-packed arrangement of the 
adatoms is shown in Fig. 9. The exclusion 
pattern of a single adatom, F, is shown by the 
black squares. Each adatom excludes 9 sites. 
One-quarter of the sites are occupied so that 
y=4. Evidently all the adatoms in any number 
of columns can be displaced one site without 
violating any exclusion conditions. For example, 
in column 6 of the figure the atoms could occupy 
the positions indicated by the vertical crosses. 
Alternatively, any number of rows could shift in 
the same way. But the displacement of either, 
rows or columns, locks the adatoms against dis- 
placements at right angles. 

The possibility of such shifts is important, 
because it means that for @ almost unity the 
number of free sites varies with (1—6)? instead 
of 1—@. This can be seen as follows : Suppose that 
from the completely covered surface the adatom 
at A evaporates. The free site A can then be 
excluded by the inward motion by one site of 
B and D, singly or both, or by C and E similarly, 
followed by any number of other adatoms in the 
particular row or column. This number is limited 
only by the locking already mentioned, which 
can occur when other adatoms besides A have 
evaporated. As a result, the single free site can 
be split in many ways into two gaps each 3d, long 
in a line of adatoms whose normal spacing is 2d,. 

Each such gap in an adatom row represents a 
deficiency of $ atom below saturation. Hence 


TONKS 


Ap 


Fic. 10. Adsorption isotherms. Designations as in Fig. 6. 


their total, s, numbers 
(30) 


With a number of free sites the interlocking 
becomes much more complicated, because one- 
site shifts in the direction of both rows and 
columns become possible. As a result, single gaps 
are to be found in a chain terminated at both 
ends by locking. (The pseudo-close-packed @=0.8 
lattices which have already been mentioned 
represent the ultimate in this direction.) This 
means that the problem is one of finding the 
number of adjacencies of two such gaps to form 
a free site on a surface rather than in a row. 

We note that each gap is a gap with respect 
to a row or column but not both. In this sense 
it is polarized. Only two similarly polarized gaps 
can combine to give a free cell. Since half the 
gaps are polarized horizontally, half vertically, 
we have to consider the combinations of only 
s/2 gaps at a time. Each of the s/2 horizontal 
gaps finds either another one of the s/2 hori- 
zontal gaps or one of the o adatoms adjacent to 
it on the right. Thus the chance that it finds 
another gap is 3s/(0+s/2), and for all s/2 gaps 
there are {s?/(0+5s/2) gap pairs (without duplica- 
tion) and hence that number of free sites. 
Including the vertical gaps, 


os =35°(0+5/2). 
Elimination of s with Eq. (30) gives 
os =201(1—8)?, 
whence the isotherm Eq. (10) takes the form’ 
9=4Ap(1—6)2/2, (31) 


s=2(0,—0). 


6 near unity. 


The analysis of the dilute film taking isolated 


7 The difference in coefficient of Ap between Eqs. (31 
and (32) here and Eqs. (50) and (51) of reference 1 arises 
from a difference in convention regarding A. Here 4A 15 
defined by Eq. (10.5). 





(30) 
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singlets and the three possible types of doublet 
into account by the method already outlined 
gives’ 

6=4Ap(1—2.250+ 6"), @near zero. (32) 


In Fig. 6 curves B; and Bz are plots of 6; in 
these equations. 

Two other 6; plots have been added to Fig. 6 
in order to fill out the picture of the change in 
the isotherm as y increases. Curve C is 6; for 
7y=1 in which case each adatom excludes only 
the site it occupies, and for which the isotherm is 


0=Ap(1—8) (33) 


over the whole range. Curve U is for adsorption 
on a uniform field,* corresponding to y= ~. 

The isotherms have been plotted in an 
unusual way in order to exhibit their behavior 
the better with respect to the extremes of con- 
centration. In Fig. 10 these same isotherms have 
been plotted in the usual fashion @ against 
pressure. 

In conclusion, I want to thank Dr. Langmuir 
for having proposed this investigation and for 
consulting with me during its progress. 


8L. Tonks, Phys. Rev. 50, 955 (1936), Eq. (28); or 
Eq. (35) of reference 1. 
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Potential Function of Acetylene Molecule. II 
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From a study of the interaction between the vibration and the rotation, the convergence 
factor of the rotational lines of the vibration-rotational bands of C2H: is obtained as a function 
of the various vibrational quantum numbers and the constants in the first-order anharmonic 
potential. A tentative calculation of the first- and second-order anharmonic potentials is 
carried out for the C,H» molecule on the basis of the classification of the absorption bands 


given in Part I. 


INTRODUCTION 


N Part I of this paper,! the coefficients 

V1, *** V5, X11, ***X45 in the vibrational energy 
expression of C.H»2 have been obtained in terms 
of the constants in the first- and the second-order 
anharmonic potentials. It was found that there 
are 30 unknown constants 1, --+w5, a, °-°k, 
8;, --+Bi; in the most general quadratic, cubic 
and quartic potentials that are consistent with 
the geometrical symmetry of the molecule while 
there are only 22 relations connecting these con- 
stants and the quantities 1, -+-v5, X11, °° *%45. 
In the following, an attempt is made to obtain 
five more relations between these constants and 
the convergence factors of the vibration-rota- 
tional bands from an analysis of the interactions 
between the vibrational and the rotational 
motions of the molecule, and to obtain an es- 
timate of the potential constants. 


_'T. Y. Wu and A. T. Kiang. J. Chem. Phys. 7, 178 
(1939), referred to hereafter as I. 


INTERACTION BETWEEN VIBRATION AND 
ROTATION 


The general theory of the energy of a vibrating- 
rotating molecule has been given by Wilson and 
Howard,” and for a linear Y-X—X-Y molecule 
in particular by Sayvetz.* The kinetic energy 
part of the Hamiltonian is 


1 Mas 
eV bbe — Pu FH3V/ DL Na—Ag 
er aa eer man 


+3 2 Hap MaMs— 2 HapMoAs 
1 Mas 

—3VuL Mal As—}, (23) 
a8 Ve 


in which M,, M,, M. are the components of the 
total angular momentum, A,, A,, Az are the com- 
ponents of the angular momentum arising from 

2E. B. Wilson, Jr., and J. B. Howard, J. Chem. Phys. 


4, 260 (1936). 
3 A. Sayvetz, J. Chem. Phys. 7, 383 (1939). 
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the vibrational motions, along the x, y, z axes of 
a coordinate system rotating with the molecule, 
and the ~’s are the momenta conjugate to the 
normal coordinates. The first term represents the 
kinetic energy of vibration and has been treated 
in I. The second term represents the energy of 
interaction between the parallel and the per- 
pendicular vibrations and has been given by 
Sayvetz.* These two terms give rise to energies 
which are functions of the vibrational quantum 
numbers alone and hence do not affect the rota- 
tional structure of the individual bands. The 
third term represents the energy of rotation of 
the molecule with moments of inertia 1/pz2, 1/yy 
etc. which are functions of the normal coordi- 
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nates. The last two terms involving the products 
of the M’s and the A’s represents the coupling 
between the vibrational and the rotational 
motions of the molecule. Here we shall be chiefly 
concerned with these last three terms of the 
Hamiltonian (23). 

With the notations of J with the exception of 
R; and R: which are now defined as (u;)! and 
(u2)! times, respectively, the Ri and Rzin J, where 


2mM 2mM (a+b)? 
~(M+m)  Mb®+m(2a-+b)? 
a=r(C—H), b=r(C—C) 





M1 


it can be shown that the A’s are given by 








Az=C1a(u1)*(yabs— Osh.) +5(¥sh1— QiPys) HH Vsp2—Qahus), 
Ay= —C13(u1)*(xaps— OsPz.) — $(x5P1 — QiPzs) — t(xsp2—Qeps,), 
Az= (XaPy,—YaPrs) — (X5Pys— Yoh), 

- (bc11— C21) (ue)? c (bc12— C22) (ue)? 


§ ’ t ’ 


a+b a+b 


P+P=1, C13241=1 (25) 





and x, is the x component of R,, etc. 

To obtain the moments and products of inertia as functions of the normal coordinates, it is 
necessary to picture the perpendicular vibrations as involving changes in the bond angles alone, the 
bond distances remaining unchanged during these vibrations. Denoting by ¢: and ¢:2 the angles 
which R; and R,2 make, respectively, with the x axis, we obtain 


m(2a+b) 


A 1 R,? (sin? 
2 }- 141 +aQi+6Q2+70?+ eee 
0 


ual Io \cos? gi 





Mb+m(2a+5b) at , R,? (sin? “| 


WP rcomene 
2Mb aly Ig cos? go 


C=R,?+ R:?, D=R; sin ¢1 COS gitR? sin ¢g2 COS ¢2, 
$1 Sin ¢2 
+(sQs-+10) Re} 


’ 


E «an - 
aa 2° 


(COS ¢1 


(COS ge 


where 


aly=2m(2a+b)ci+[Mb+m(2a+d) Jeo, 
BIo=2m(2a+b)c1.+[Mb+m(2a+5d) |coo, 
ylo= 2me11(611 +621) +3(M+m)co1?, 
619 = 2mc12(C12+C22) +3(M+m)c22?, 


and Jo>=3[Mb?+m(2a+5)*] is the equilibrium value of the moment of inertia. Calculation gives for 
the coefficients y’s in (23) 


Mez=1/A’, My =1/B’, Mry~9, Myz=Mez=0, M,—A.,=0, 








oducts 
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(24) 
(25) 


3, it is 
ne, the 
angles 


in? * 
OS? ve 


(26) 


ives for 
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where 

s? t? 2st 
A'=B'=1 1+00:+80:+(1-—)ovt+(s-—)or-— 0,0 

Io Io Io 


m(2a+b) Mb+m(2a+b) (2a+6) 


Ri? Re (27) 
malo 2Mb aly 








The contributions to the energy from the terms Ag(utas/+/“)Ag are, up to the second approximation 
of J, 


h2 





ey, [ (ws, wa) (¥s+ 3) (va +1) +57(01 +4) (05 +1) (wr, ws) +27(02+43) (V5 +1) (we, ws) —2], (28) 


where 


(wi, w;) =(w;/wj;+;/w;). (29) 

The effect of these contributions can be taken account of in the energy expressions (17) and (18) in I 
simply making the following replacements, with By=h/8x?J0oc, 

Bo Bo+s*Boc(w1, ws); Bi2z—>BietlBoc(we, ws); Bis —Bist+Boc(ws, wa). (30) 


The term 2yegM.M¢ gives, in a first approximation in which the yzz, wyy are given the equilibrium 
values, the energy of rotation of a rigid rotator 


(h?/82?Io) (J (J+1) — (it.)?]. 


The interactions between the vibration and the rotation of the molecule appear in the next approxi- 
mation. First, taking into account the variations of yz:, uy, with the normal coordinates gives the 
energy change 








h? s? {2 
AE oe= ———[J(J +1) (+ 1)*Y 001+802+( 1-at——)orr+( i—p——)os 
827Jo Io Io 
st 2a+bsm Mb+m(2a+5b) 
~2(a8+~)o.0.- (— 2 Rt) (31) 
I alo \m 2Mb 


The quadratic terms introduce energy changes that can be calculated by means of the perturbation 
theory. The linear terms in Q; and Q2 change the equilibrium separations of the atoms by amounts 
determined by the equations 0V/0Q:=0V/0Q.2=0 where 


V=} > i AQ? + AE rot, 
1 


namely, - oi 
a ieee 
4Q1=— [J(J+1)—(itds)?], AQe=— 


d; 82?Io 2 81*lo 








LJ(J+1) —(ith)?]. 


These stretchings necessitate a transformation to the new normal coordinates Q*, and Q*: defined by 
g 


Q*1=Q:—AQi, Q*2=Q2—AQ:. (32) 


As a consequence of this transformation, the first-order potential (14) in I contributes the following 
terms which are quadratic in the normal coordinates 





h? a B B a a B 
av= [II+1)— (ht )*] (3k.-+k— )ort+ (hha )On+2(k— +h) ONO 
827] At Ae de vi 1 Ae 


a B a B a B 
+ (tthe )ort+ (a +h—)ot+ (tothe) ot (33) 
A Ae M1 Ae Al he 
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in which 
ka=8r'w1(w1/h)'a, ky = 87? wo(we/h)'d, ke = 825w1(we/h)'c, ka= 82 wo(w1/h)'d, 
ke =8rw3(wi/h)te, kp=S8m wayr(wi/h)*f, kg = 82? wsu2(w1/h)'g, (34) 
kh, =81rw3(we/h)i, ky =SB8mwayr(we/h)'j, ke = 82? w5u2(we/h)R. 


Next we have to consider the contributions from the last two terms in (23). The last term vanishes 
since prz/(u)!=p,,/(u)t=1 and the operator A operates only on quantities within the parentheses. 
The diagonal matrix elements of M.Aq also vanish; but in the next approximation, the nondiagonal 
elements contribute the following energy change 


AE=— >’ (VR| M.A.+M,A,| V'R’)(V’R'| M.A. +M,A,| VR) 
Io? v’R’ Ey®— Ey 
2(vst+3)—(vat1) = 2(ar+3)—(Us +1) — 2(v2+3)—(s+1)] hh 
+s? +72 
8271 





h? 


— [I +1)—(h+i*] 


8171 





W3—-W4 W1i—W5 Wea—-W5 
+terms independent of J. (35) 


In obtaining these second-order contributions, the rotational energy has been neglected in com- 
parison with the vibrational energy in the denominator and use is made of the relation Zz/(R| M,| R’) 
< (R’| M,|R)=(R| M22|R), etc. The total energy change is, on adding (31), (33) and (35), 


h2 3 5 
Reet 1)—(itde)* ILE (itd) ait dD (vit ai, (36) 


h s*h 


4r*w, 42?Io(wi—ws) 





s* a B 
a= («4-438 +h) 


0 v1 Ae 





t? B a h th 
c= (s+ — 543 +h~) + 
0 


4n*w, 42?Io(we—ws) 


he Ai 


a B h h 
a=(t—+h—) + ’ 
M1 Ne 47r*w3 4m? I o(w3—w4) 





m 2a+b a B h h 
a= — +kp—+k;— ~ , 
Mi ao Aur =— Ager 4rr2wq 82D (w3 — wr) 


Mb+m(2a+b) 2a+6 a B h h s? t? 
as=( ’ +k ( ). 
w1 








o—— +k - + 
2Mb alo Are opel 422w; 827Io —W5 We—-Ws 





The frequencies of the vibration-rotational lines arising from transitions from the normal vibra- 
tional state to a state v4, v2, v3, V4, 11, v5, 12 are given by 


3 5 5 
is Vyib+Bol 2+ >, (viti)a:t+> (v;+2)a;|N+Bo(> v;a;) N?, (38) 
1 ° 4 1 


where vi, is determined by (17) in I together with (30). The a’s can be determined empirically from 
five properly chosen and accurately analyzed bands. Eqs. (37) then furnish five more relations for 
the determination of the potential constants in (14) in I. 





vibra- 


(38) 


y from 
yns for 
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TABLE I. 








ACCORDING TO I 


1 v2 v3 Bi—Bo oss.| Bi CALc. 





—0.0179 . 
— 0.0238 ° 
—0.0250 | 1.1618 
—0.0254 3 

—0.0294 | 1.1467 
—0.0329 | 1.1435 
—0.0406 | 1.1360 
—0.0378 | 1.1376 


eee OO 




















-OOonNn CK KOS 








CALCULATION OF POTENTIAL CONSTANTS 


The transformation coefficients c’s and the 
quantities a, B---, ka, -:-k, can be calculated 
from the zeroth-order frequencies and potential 
constants in I. 


¢4=1.125, €12=5.57, C1i3> — 5.70, Co = 2.81, 
Coo= —1.53X10"; s?=0.037, t?=0.963; 
a=4.06X10'%, B=0.689X10!9; 
vy =4.25X10*8, 6=4.2510*; 
kg=4.3 X10*%a, ky = 10.3 X 104%, 
k.=5.8X10**c, ka=7.7 X10*8d, 
k.=7.35 X10**e, ky =1.56X 104% uf, 
kg=1.30X10*% og, k:=9.8X 104%, 
k;=2.07 X10*%u17, k;,=1.74X10** yok, 
where the a, b, ---k are in cm and the other 
quantities in c.g.s. units. 

To obtain the a’s, consider first those bands 
involving only the parallel fundamentals »;(v), 
vo(¥5), ¥3(va). The values of B, for the upper vi- 
brational states are known from the analyses of 
Funke and Lindholm.‘ Table I gives eight such 
bands. Of these, the B,’s of the 9640, 11,600, 


12,676 cm~ bands are probably most accurately 
known. From them we obtain® 


a2By= —0.0075, 
a3By)= —0.0029. (40) 


a,By= — 0.0059, 


The B,’s of the other bands calculated with these 


*G. W. Funke, Zeits. f. Physik 104, 169 (1937); Funke 
and Lindholm, ibid. 106, 518 (1937). 

*These values of the a’s can be compared with the 
empirically chosen values a;Bo=a2Bo=a3Bo= —0.0051, 
«4By=asBy=0.0020 of Funke and Lindholm, reference 4. 


values are given in Table I. Except for the 
11,783-cm™ band, the agreement between the 
calculated and the observed values is satisfactory. 

From the bands 73 +2v2+72+2v;=12,711 cm", 
B,=1.156 and 27;+2ve+73+2v;=14,597 cm“, 
B,=1.151 and the values in (40), we obtain the 
average value® 


a4By = 0.00193. (41) 


With the values in (40) and (41) and the follow- 
ing bands 


Vyyt2vetvt v,= 9906 cm, Byo=1.162 
Vit3votvgt v5 = 13,033 1.1501 
vit 2vgtvet+ v5s= 9745 1.166 
2vit2v3stvst+ v5=11,663 1.156 
3yyt+ votvyt3v;= 11,586 1.1584 
3vy;+ vot vgt+5v;5 = 12,632 1.1643 
we obtain the average value® 


a;By= 0.00245. (42) 


The agreement between the calculated and the 
observed B,’s for bands involving the perpendicu- 
lar fundamentals is not as satisfactory as in the 
case of those bands involving the parallel funda- 
mentals alone. This may be due to the uncer- 
tainty in the assignments of some of the bands 
or to the lower accuracy with which the B’s are 
determined. 
Substituting the values in (39)—(42) into (37) 

gives 

0.363 a+0.0090 c= —30.6 

0.048 b+0.216 d=-—46.1 

0.206 e+0.0152 «= —20.6 (43) 

0.437 f+0.0321 j=— 4.06 

0.366 g+0.0272 k= — 8.34. 


Eqs. (18) in I and (43) are still not sufficient for 


TABLE II. Anharmonic potential constants of C2H2 in cm. 








= —110 for 


b=— 200 d= —167 
— 168 


e-— 85 i~ 
S- 11 


g—— 25 
B= 54 
By — 23 
Bie— 17.5 
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the evaluation of the potential constants in (14) 
and (15) in I. However, by making certain 
reasonable assumptions about some of them, a 
rough estimate’ can be made of the others as 
shown in Table II. These values are admittedly 
very rough, but any accurate calculation will not 
be possible without additional empirical data. 
The order of magnitude of these constants shows 
that the effect of anharmonicity is not quite 
negligible in C2.H». This is in fact what we would 
expect from the small mass of the hydrogen atoms 
in the molecule. 

The writer is grateful to Professor D. M. 
Dennison, of the University of Michigan, for his 
interest and advice concerning these calculations. 


TABLE III. Classification of C2He bands. 











ba 
2va+vs+5a 
Ya t2K. +é5s 


? 
2va+vot+batds 
2va+2v0+5a +5 
va tvs t+vo 3va +v0+25¢ 
va +2vs+v0 3vqa+2v0+25s 
3va+vo votda tds 
3vatvstvo va t+vs+v0+25s 
Svatvo Va tve +20 +255 
3va+2vo ? 
va +5 va +e +v0+5s 
vetba 3va+vs+v0tds 
vot+dba Sva+vot+is 
ba +5s - - 2va+vo0+5a 


vs +ba+5s Va vo 

3vat+5s —5s vg t+3v0+6a+35s 
3vq +s +5s —5s ¥g +3v0+5a +55 
va +23 +53 —5s 80 3va+vs tis 

5vat+5s —5s 15521* 























6 The way in which we arrive at these values is as 
follows. From the second equation in (43) and the equation 
for x22, we obtain b and d for a few small values of Be. 
It is found that b and d are not very sensitive to B2. From 
the first equation in (43) and the equation for xu, it is 
found that there are no real solutions for a and ¢ unless 
81 <1.575, and then c is very sensitive to the value of 6: 
while a is not. In a similar manner, e, i, f, j, g, & are esti- 
mated, it turning out that 8;>20/3 for g and k to be real. 
While the values of e, f, g are not very sensitive and their 
values in Table II may be taken as correct within 10 cm“, 
the values of 7, j7, & are very sensitive to the values of the 
6’s assumed and are hence quite uncertain. 


TABLE IV. 








v1 = 1987 xu = —13 xz=— 5 
v2 =3358 xe2= 13.6 xs=— 5 
v3 =3318 x33 = —28 xu= 
vatvua= 729 xs= 9 xus=— 17 
vs= 603 x333= — 2.3 xe3 = —112 
x555 = — 0.55 








TABLE V. 
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15,600 
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APPENDIX 


In I, it was stated that with the values of the coefficients 
v1, ***¥5, X11, ***X45 given by Herzberg and Spink,’ the 
agreement between the calculated and the observed fre- 
quencies is not satisfactory. However, by modifying the 
assignments of those bands involving the perpendicular 
fundamentals, it is possible to obtain good agreement for 
all except 5 or 6 bands, as shown in Table III. The bands 
marked with * are those employed in the calculation of 
the coefficients v1, -++vs5, X11, ***Xas which are given in 
Table IV. It is seen that the agreement between the ob- 
served and the calculated frequencies is good except in 
the case of the 6500, 8450, 15,081 and 9151-cm™ bands. 
The discrepancy in the case of the 6500-cm™ band, 
namely, ca. 50 cm™, must, however, be regarded as a 
serious weakness of the classification, since the agreement 
is expected to be best for such a simple combination band. 
Also the two parallel type bands at 9906 and 13,033 cm™ 
and the perpendicular band at 9366 cm™ do not find 
reasonable assignments in accordance with the selection 
rules. Another reason for our preferring our classification 
given in I to the above modified Herzberg-Spink system is 
the following: According to the latter classification, the 
convergence factors of the bands in Table I now become 
as shown in Table V. It is seen that the agreement between 
the calculated and the observed B,’s is not as satisfactory 
as that obtained in Table I. It seems that to decide among 
the various classifications for the C2H: bands that have 
been proposed, further experimental work, preferably on 
C2Dz, is necessary. 

7G. Herzberg and J. W. T. Spink, Zeits. f. Physik 91, 386 (1934). 
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It has been found that positive ions are emitted by filaments of tungsten, tantalum, and 


molybdenum when heated to temperatures near the melting points. Sodium and potassium ions 
were found in addition to those of the filament metal. When alkali halides were heated to 
800-1000°C on platinum they emitted the singly charged ions of the corresponding alkali 
metal. No negative ions were found and no ions were obtained from the chlorides of calcium, 


magnesium, strontium, barium, or zinc. 





HE emission of positive ions has been 

studied by several experimenters. The early 
work is summarized by Richardson! and a later 
summary has been given by Jones.? The most 
detailed studies using a mass spectrograph to 
identify the ions have been by Wahlin* and by 
Barnes.‘ They have established that in addition 
to the ions of sodium and potassium which had 
been observed earlier, ions of the filament ma- 
terial are emitted. We have made some observa- 
tions on the variation in intensity of the alkali 
ions and ions of the filament metal with time and 
also on the nature of the ions produced by 
heating a number of salts. 

The apparatus used was essentially the same as 
described by Hogness and Lunn.® The magnetic 
field was calibrated by measuring the change in 
resistance of a bismuth spiral. It was found that 
starting from zero current in the magnet coils, 
the current-field curve was quite reproducible. 

In order to prevent contamination by fluxes, 
the filaments were not brazed or welded to the 
electrodes but were clamped between heavy 
tungsten leads in the following way. Each lead 
was made of two heavy tungsten rods ground 
flat on one side. These were brazed two inches 
above the ends with flat sides together. The 
filament could then be inserted by springing the 
ends apart and these, when released, held the 
filament tightly clamped. 





‘O. W. Richardson, Emission of Electricity by Hot 
Bodies (Longmans, Green, 1921). 

?T. J. Jones, Thermionic Emission (Methuen and Com- 
pany, London, 1936). 

*Wahlin, Phys. Rev. 34, 164 (1929); 35, 653 (1930); 
45, 886 (1934). 
‘Barnes, Phys. Rev. 42, 487, 492 (1932). 
— and Lunn, Proc. Nat. Acad. Sci. 10, 398 
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The ratio for the charge to mass of the ions was 
determined by keeping the magnetic field con- 
stant, varying the accelerating potential applied 
to the ions, and noting the electrometer deflec- 
tions. The peaks in the curve obtained by 
plotting the electrometer deflections against the 
accelerating voltages give the values of e/m by 
application of the formula e/m=2V/H?R? in 
which H is the magnetic field, R the radius of 
curvature of the path of the ions in the field, and 
V the accelerating potential applied before 
extrance of the ion into the magnetic field. The 
resolution of the apparatus for ion currents of the 
strength obtained was not sufficient for the 
separation of the isotopes but was adequate for 
the determination of the nature of the ions. 

Positive ions were obtained from tungsten, 
tantalum, and molybdenum filaments when these 
were heated to temperatures slightly below their 
melting points. All of the ions were found to be 
singly charged. Sodium and potassium ions were 
obtained from all of the filaments in addition to 
the ion of the filament metal. The tungsten 
filaments also gave some thorium ions. The 
intensity of emission of the heavy metal ions was 
constant at a constant temperature for any given 
filament or like filaments. The sodium and 
potassium ion currents, however, showed a falling 
off as time progressed so long as the ions were 
removed by application of a potential. When a 
filament, from which the emission had decayed so 
far as to be unmeasurable, was heated over night 
at a lower temperature without the application 
of a potential the current regained nearly its 
initial value. The rate of decay varied with the 
different metals. By selecting filament tempera- 
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ture such that the intensity of the ionic current 
was the same, it was found that the times 
required for the emission to fall to its half-value 
were ten, fifteen, and one hundred and forty 
minutes for molybdenum, tantalum, and tung- 
sten, respectively. These half-time values cannot 
be interpreted in any simple manner since it was 
found that the decay of the emission did not 
follow any simple rate law. It is probable that the 
decay is due to an exhaustion of the sodium and 
potassium on the surface of the filament. The 
diffusion of more alkali atoms from the interior of 
the filament is too slow to supply a measurable 
current. 

A large number of experiments were tried with 
salts heated on a platinum spiral to temperatures 


of 800 to 1000°C. It was found that the alkali 
halides gave the corresponding positive ions. For 
example, the chlorides of sodium, potassium and 
rubidium gave ions the ratio of whose e/m was 
1/23 : 1/39 : 1/85 within two percent. None of 
the halides of magnesium, calcium, strontium, 
barium, or zinc gave positive ions. No negative 
ions were obtained from any of the salts. The 
sharpness of the peaks indicated that all of the 
ions fell through the same potential drop and 
hence it may be concluded that they were formed 
either at or very close to the filament. Under such 
conditions the work function of the platinum is 
high enough to remove the extra electrons from 
any of the negative ions which might have 
formed. 
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Entropy of Vaporization and Restricted Molecular Rotation in Liquids 


R. S. HALFoRD 
Department of Chemistry, University of California, Berkeley, California 
(Received December 18, 1939) 


Deviations among entropies of vaporization to constant vapor volume are attributed to 
restricted rotation of molecules in the liquid. An arbitrary measure of the amount of restriction 
is defined and the deviations are shown to exhibit a remarkable regularity when compared on 
this basis. Deviations due to perturbation of vibrations in the liquid must either parallel those 


due to restricted rotation or be quite small. 


HEN a liquid is caused to vaporize iso- 
thermally against its own vapor pressure 
the entropy of vaporization is given by the well- 
known relation AS=AH/T. AH, the heat ab- 
sorbed, is a characteristic property of each 
liquid which has approximately the same value 
over a wide range of temperature for any given 
liquid. A characteristic temperature can there- 
fore be defined for each liquid such that, at 
their respective characteristic temperatures, all 
liquids have the same entropy of vaporization. 
Three important attempts have been made to 
relate the characteristic temperature with other 
properties of the liquid. The widely used rule of 
Trouton! identifies it with the boiling point of 
the liquid. The inadequacy of this definition is 
evidenced by the fact that the corresponding 


1F, Trouton, Phil. Mag. (5) 18, 54 (1884). 


entropies of vaporization show a trend when a 
wide range of boiling points is considered. 
A second definition of the characteristic tem- 
perature is deduced from the theory of corre- 
sponding states. According to this definition, 
liquids at their characteristic temperatures have 
the same ratio of vapor volume to liquid volume. 
A recent discussion, given by Pitzer,? makes 
evident the fact that this rule is not generally 
obeyed although a few simple substances show a 
very striking agreement. Large differences in 
entropy of vaporization are obtained with sub- 
stances selected at random and Pitzer has dis- 
cussed the significance of these deviations. 

The third definition of a characteristic tem- 
perature was proposed by Hildebrand* as 4 


2K. S. Pitzer, J. Chem. Phys. 7, 583 (1939). 
3 J. H. Hildebrand, J. Am. Chem. Soc. 37, 970 (1915). 
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modification of Trouton’s rule. Liquids are in 
corresponding states when all are in equilibrium 
with the same concentration of vapor. This rule 
meets with considerable success in practical 
application; there is no trend with temperature 
and deviations among so-called ‘‘normal”’ liquids 
are relatively small. Several attempts‘ have been 
made to furnish the Hildebrand rule with a 
satisfactory basis in theory. It appears, however, 
that theoretical considerations always lead to 
the second definition given above. Excepting the 
few simple substances which were shown by 
Pitzer to be in perfect accord with theory, this 
definition of characteristic temperature cannot 
compete in practical value with the empirical 
rule of Hildebrand. A very attractive suggestion 
concerning the solution of this paradox for the 
case of complicated polyatomic molecules was 
made to the author by Professor Pitzer. The 
theoretical treatments leading to a definition of 
corresponding states have assumed that the 
potential energy can be represented as a function 
of the distance between molecular centers. This 
assumption will lead rigorously to the definition 
in terms of equal vapor to liquid volume ratios. 
However, the most polarizable electrons in the 
molecules will be the nonbonding valence elec- 
trons on outside atoms and the major contribu- 
tion to the potential energy will come from these. 
This suggests that the attraction between neigh- 
boring molecules may arise largely in the inter- 
action of adjacent outside atoms; in which 
case, the potential energy should be represented 
as a function of the distance between these 
atoms. It can be seen qualitatively that this 
modification of the theory should lead to some- 
thing similar to the Hildebrand definition of 
characteristic temperature. I trust that Pro- 
fessor Pitzer will find the opportunity to elaborate 
on this suggestion soon. In the meantime it will 
be tentatively assumed that the comparison of 
entropies of vaporization at equal vapor volumes 
is justified. 

Recently, Hildebrand® has called attention to 
certain regularities among deviations from his 





*(a) R. F. Newton and H. Eyring, Trans. Faraday Soc. 
33, 73, 80 (1936); (b) O. K. Rice, J. Chem. Phys. 5, 353 
(1937); (c) H. Eyring and J. Hirschfelder, J. Phys. Chem. 
41, 249 (1937); (d) J. E. Lennard-Jones, Proc. Roy. Soc. 
London, A165, 1 (1938). 

°J. H. Hildebrand, J. Chem. Phys. 7, 233 (1939). 
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rule and has suggested an interpretation for 
them. Assuming that the rule would be obeyed 
for the vaporization of a completely random 
liquid to a completely random vapor, he points 
out that any tendency toward order in the 
liquid (with consequent lower liquid entropy) 
would increase the entropy of vaporization to a 
random vapor. Among “normal” liquids the 
factors tending to produce ordering would be 
polarity and geometrical shape. Hildebrand 
showed that the higher entropies of vaporization 
were observed for those liquids which were 
most likely to be ordered for just these reasons. 

In the subsequent discussion these deviations 
will be examined from the point of view of 
hindered rotation of the molecules in the liquid. 
The suggestion that this phenomenon may be 
the cause of deviations is not a new one.? 4(%») 
The relation to be demonstrated should be re- 
garded as empirical. The concept of restricted 
rotation may be synonymous with the ordering 
proposed by Hildebrand but it is not necessarily 
so. It is possible to imagine molecules, so packed 
that dipoles or other identifying axes are oriented 
at random, but prevented by their shape from 
executing a complete rotation without inter- 
ference from their neighbors. 

If rotation and vibration, respectively, are not 
equally excited in liquid and vapor a contribu- 
tion to the entropy of vaporization can be 
associated with these degrees of freedom. The 
contribution from rotational motion can be 
described qualitatively in terms of a simple 
model consisting of an assemblage of molecules 
occupying a certain volume which is subdivided 
into cubical compartments. The molecules are 
distributed with one in each compartment and 
the center of gravity of each molecule is fixed to 
coincide with the center of its compartment. 
Initially it will be supposed that the molecules 
are rotating freely with a Maxwellian distribu- 
tion of rotational energy in the assemblage. 
Each compartment is thus large enough to per- 
mit a molecule to sweep out a complete rotation 
within it. Under such conditions the rotational 
factor for the assembly is the same as that for a 
gas. Now, if this system is slowly compressed so 
that the volumes of all compartments are 
diminished similarly, there will be no change in 
the rotational factor until a certain critical 
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volume is reached. At this point the side of each 
cubical compartment will be twice the length of 
the longest rotating arm of the molecule. Any 
further compression will result in interference 
among the rotators and a loss of rotational 
entropy. Torsional oscillation about the diagonals 
of the cube will persist in diminishing degree 
over a range of volume and a continuous loss of 
rotational entropy will accompany the com- 
pression. This phenomenon has been discussed 
quantitatively by Pitzer.* The model used here 
is of course only qualitative and highly simplified 
but it furnishes a basis for the comparison of 
actual liquids. 

One method of comparison would be to com- 
pare the entropy of vaporization with the 
volume of a liquid at various temperatures, 
always correcting to the same vapor volume. 
However, the amount of expansion occurring 
over the range of temperature where the vapor 
remains dilute would correspond to only a small 
portion of the process discussed in connection 
with the model. The region of concentrated 
vapor introduces vapor imperfections as an 
undesirable complication. 

On the other hand, if the notion of hindered 
rotation is to be applied quite generally the 
phenomenon should occur in very much the same 
way for all liquids, allowing for individual 
differences. An estimate of the importance of 
these differences can be taken from the work of 
Pitzer who showed that over a wide range of 
moments of inertia the entropy loss due to a 
given amount of restriction of rotation is about 
the same. The individual differences are com- 
parable in magnitude to the errors in the best 
available entropy values for the liquids in which 
we are interested. Thus we may expect that in 
corresponding states of restriction, in the sense 
of the model, the various liquids will have lost, 
to a good approximation, the same amount of 
rotational entropy. If then, the various liquids 
are found to cover a wide range of restriction 
under the conditions that their entropies of 
vaporization have been determined it will be 
possible to carry out a comparison on the basis 
of the model. 

In order to do so it is necessary to be able to 
define suitable corresponding states of restriction 


*K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). 
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for various liquids. In terms of the model these 
states refer to the volume required by a molecule 
in order that it may sweep out a complete 
rotation. It is difficult to define this volume 
exactly. For a rigid molecule, the distance from 
the center of gravity to the center of the extreme 
atom along the longest rotating arm can be 
given quite accurately. In many cases, accurate 
electron diffraction data are available; otherwise, 
distances can be estimated satisfactorily from 
covalent radii. However, the additional extension 
of the extreme atom away from the center of 
gravity is not so simply determined. It should be 
greater than the covalent radius and possibly 
somewhat less than the ionic radius for the same 
atom as a singly charged negative ion. Estimates 
of such distances are available but they are not 
sufficiently accurate for the purpose at hand. 

In view of this difficulty, corresponding states 
of restriction will be arbitrarily defined in terms 
of a volume calculated by assuming that the 
extension of the extreme atom away from the 
center of gravity is equal to its covalent radius. 
The length, so defined, in turn defines a sphere 
corresponding to complete rotation of this length 
about the center of gravity. The molal volume 
calculated for cubic closest packing of these 
spheres will be called the molal volume of the 
free covalent rotator. The ratio of the calculated 
molal volume of the free covalent rotator to the 
actual molal volume measures the degree of 
restriction of the liquid. Liquids having the same 
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degree of restriction are said to be in correspond- 
ing states of rotational freedom. 

The entropy of vaporization.to a vapor volume 
of 49.5 |. has been plotted against the degree of 
restriction for a number of liquids in Fig. 1. 
The close grouping of the points about a single 
simple curve seems highly significant. The devia- 
tions from the curve as drawn are less than 
might be expected from the errors in the values 
of the entropies, making no allowance for the 
individual differences due to different moments 
of inertia. The horizontal line at a value of 
20.1 E.U. is the average for several monatomic 
substances which agree among themselves within 
0.2 E.U. The fact that this entropy value is 
achieved at a degree of restriction of 0.7 rather 
than unity is due merely to the arbitrary 
definition chosen for the length of the rotating 
arm. In calculating these lengths all data have 
been taken from Pauling, The Nature of the 
Chemical Bond.? When in a terminal position 
methyl and methylene groups have been assigned 
an extension of 1.37A, equal to the covalent 
radius of C plus the covalent diameter of H. 
The entropy values were taken from the vapor 
pressure curves generously provided by Pro- 
fessor Hildebrand and used by him in the 
preparation of his paper.’ The entropy values 
given by the author differ slightly from those 
given by Professor Hildebrand and the differ- 
ences, about +0.1 E.U. represent the limit of 
accuracy in reading the plots. The circles en- 
closing the points in Fig. 1 correspond to this 


7 Cornell University Press, Ithaca, New York. 
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uncertainty in the entropies and to an uncer- 
tainty of about 0.01A in the length of the rotating 
arm. The actual experimental uncertainties may 
be several times larger in both quantities. These 
figures are cited to call attention to the truly 
striking consistency of the points in Fig. 1. 

A few remarks may now be made concerning 
the curve which has been drawn.to fit the points. 
Pitzer’s calculations® lead to a curve of this 
shape when the entropy lowering due to re- 
stricted rotation is plotted against the value of 
the restricting potential. Since the relation be- 
tween the restricting potential and the degree 
of restriction, as defined here, is unknown, a 
quantitative application of Pitzer’s theory cannot 
be undertaken. However, it is to be expected 
that the two quantities would be related in a 
rather simple way and as a result the shape of 
the lines in the two types of plot should have 
some resemblance. While these remarks can 
hardly be taken as evidence for the interpretation 
given in the previous discussion they nevertheless 
add to its plausibility. 

A description of the effect of the vibrational 
factor is not so readily given. One observes 
however that for the substances included in the 
figure it must either closely parallel the rotational 
effect or be very small. 

It is interesting to note that the several polar 
substances included in the plot appear to be 
displaced upward to a small extent. This may be 
indicative of the contribution of polarity to the 
restriction of rotation but the accuracy of the 
data hardly warrants such a conclusion. 
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Absolute Reaction Rates and the Photovoltaic Effect* 


EDWARD ADLER 
Department of Chemistry, Columbia University, New York, New York 


(Received April 2, 1940) 


N several recent papers!” the theory of abso- 
lute reaction rates has been applied to the 
study of electrode processes. By an extension to 
the photovoltaic effect, the experimental de- 
pendence of the photo e.m.f. upon the intensity 
of illumination can be deduced. In a later paper 
we shall discuss some of the factors complicating 
the photovoltaic effect, such as ‘‘dark”’ or ther- 
mal e.m.f. and variation of resistance with light. 
We are concerned here only with the primary 
photo e.m.f., and shall consider only the case in 
which the photo e.m.f. is positive. (The illumi- 
nated electrode becomes increasingly positive.) 
Following the treatment in,! the current from 
the electrode for the process A—A*+¢e is 


I)! =eki [A Je“(Ve-V 8/7, (1) 


where [A] is the activity of atoms A in the 
surface of the electrode or adsorbed on it, ¢ is 
the charge of the ion, Vz is the electrical potential 
of the electrode (or in the adsorbed atomic layer), 





Vc the potential at the ‘‘saddle point,” and 
kT 


k= —e-AFi*/kT_ 


h 


Similarly, for the discharge of ions in the 
process A+-+e—A, the current due to thermal 
transfer is 


Tp!" = ehg[ At Jaase@*V 0-0/4 ?, (2) 


where [At ],a, is the activity of ions A+ in the 
adsorbed ionic layer, and V, is the potential in 
the adsorbed ionic layer. 

In addition, if the electrode is illuminated by 
radiation of the proper frequency, electrons may 
be transferred from the electrode to the ionic 
layer across the potential hill, resulting in the 
same process, namely, At+e—A. If J is the 
intensity of illumination, then the rate of the 
photo-discharge is k3J and I,’’=ek;J. 

The total current, accordingly, is 


I =ehye(V c-V BIETL A] —ehge(V CV PET At] .4, — eg J. 


At equilibrium, J = Jp the diffusion current, where 


eD 
Ip=—{LA* Juas—LA*}}. 


Setting these two equal, and solving for [A* eas 


Rye~*Ve-VB/kTT A | — kg J+ (D/l)[A*] 





[A* Jaas= 


koe *(Ve—-VF)/kT 4 (D/l) 


and substituting in (4) 


eD {kA Je~*(Ve-V B)/kT — ko J — kof At eV e-V rd /kT 
I=— 


(5) 





l Ree~*(Vc—VF)/kT 4. (D)/]) 


Now, if we assume that Vg= V, the potential of the electrode, that Vr=0 and let Ve=aV, and 


note further that at equilibrium, 
V=V., 


then, 


I=0, [A ]=[A*t], 
k3J = [A Jem eaV elk? | RyeeV el kT — fo} _ 


* Publication assisted by the Ernest Kempton Fund for Physical Research of Columbia University. 


1 Kimball, J. Chem. Phys. 8, 199 (1940). 


2 Eyring, Glasstone and Laidler, J. Chem. Phys. 7, 1053 (1939). 
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Now when J=0, 


ky 
—=e VokT (Vy is “dark” e.m.f.) 
2 


k3J = ko A Jem eaV el kT f ee(Ve-Vo)/kT 1 a (6) 


and 


This equation may be simplified in two limiting cases. 

Case 1.—Where the active electrode consists of an oxide, sulfide, or halide of copper or silver, 
selenium, etc., then (V.— Vo) is of the order of one-tenth of a volt even at moderate light intensities. 
In this case, at room temperatures e*Ye—V0)/kT>1 


ond k3J = kel A Jet eraVe-Voy/kT 
or, taking logarithms and combining the constant terms, at constant T 
V.=C+B In J. 


This is in agreement with the experimental dependence of V, upon In J for active electrodes at 
moderate illumination intensities. 
Case 2.—At very low intensities, or for pure metal electrodes, giving potential differences of the 
order of microvolts or millivolts, 
e(V.— Vo) ea(V,— Vo) 
——__—_—_- <1 sania 
kT 


. 
’ 


kJ = kof A je~*a(Ve—Vo) pew 9 
: kT 





- kof A JemeaVolkr. 


e(V.— Vo) | , ea(V.— | 


kT kT 
and, combining constant terms, 


Ve=C'+B'J. 


Again, this is in agreement with the experimental linear dependence of V, upon J at small values of 
(V.— Vo), i.e., for very small light intensities or for slight photo-activity. 
Thanks are due Professor G. E. Kimball for a very helpful discussion. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section should not exceed 600 words in length and must 


Lengths and Strengths of Atomic Bonds 


It has been claimed that it is possible to construct a 
table of bond energies applicable to wide classes of mole- 
cules, which gives agreement withone one or two kcal./mole 
with the heats of combustion.' Unfortunately ‘‘resonating” 
molecules are found to require special treatment and they 
are assigned supplementary energies, arbitrarily chosen so 
that agreement with heats of combustion is obtained. This 
scheme is essentially concerned with relating molecular 
structures, as determined by x-ray and electron diffraction 
studies and spectroscopic observations, to heats of com- 
bustion. Consequently, to the individual atoms single- 
bond, double-bond and triple bond radii are given, which 
are so chosen as to give, when added together, the best 
fit with the interatomic distances found in the restricted 
class of molecules of nonresonating types for which such 
measurements are available. 

The accumulating data seem now to call for a revision 
of this scheme. The assumption that a C—O bond or a 
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C=O bond e.g. have invariably the same energy and the 
same length, even if it were permissible as a first assump- 
tion, can no longer be justified. The C—O bond in tri- 
methyl carbinol apparently is a number of kcal./mole 
stronger than in dimethyl ether or methyl alcohol.? And 
the C=O bond is apparently shorter in formaldehyde,’ 
acetaldehyde, and glyoxal,5 than in diketopiperazine. 
There are similar difficulties with the other bonds. Even 
within the class of saturated hydrocarbons, there is some 
evidence that the C—C and C—H bonds vary in strength.’ 
That the C—C and C—N bonds vary in length, even in 
cases in which none of the usual types of resonance is 
applicable, is shown by the cases of diketopiperazine® and 
glycine.’ On the whole, the data indicate that we have to 
anticipate gradations in bond energies and gradations in 
bond lengths associated with one and the same bond. 

It is therefore suggested that the essential objective is 
to construct strength/length curves for each pair of atoms, 
alike and different, based upon the single assumption that 
there is a 1-1 relation between the strengths and the 
lengths of each such bond. In a more general scheme of this 
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kind the difficulties mentioned above disappear. Further, 
we have an approach to the problem which is at once more 
direct and more comprehensive in that, regarding a mole- 
cule, in the first instance, as a set of atoms at definite 
distances apart, we view these definite lengths as connoting 
definite energies and by summing these energies obtain a 
heat of formation directly related to the structure, without 
any discussion as to the nature of the bonds. It is more 
comprehensive since, in such a scheme, benzene and other 
“resonating’”’ molecules find a place, no arbitrarily chosen 
resonance energies having to be added so as to tally with 
the heats of combustion. 

As an illustration, Fig. 1 shows a tentative carbon- 
nitrogen bond curve, based upon the bond energies and 
atomic radii assumed by Pauling.! A curve of this type for 
the C—C bond has already been attempted.® The idea of 
strength/length bond curves may also be applied to atoms 
bonded by hydrogen, and Fig. 2 shows a tentative OH: --O 
curve based upon the lengths and strengths given for these 
bonds in ice, alcohol and formic acid.! 

Full details of the work, which includes C—N, C—C, 
C—O and other strength-length curves, based upon a new 
analysis of the data, will be published shortly. 

DorotHy WRINCH 


Davip HARKER 


Department of Chemistry, 
The Johns Hopkins University, 
Baltimore, Maryland, 
May 8, 1940. 


1 Pauling, The Nature of the Chemical Bond (1939). 

2 Pauling, J. Am. Chem. Soc. 54, 3570 (1932). 
as —_— LaValle and Schomaker, J. Am. Chem. Soc. 61, 2508 
oo Burnham and Schomaker, J. Am. Chem. Soc. 61, 2922 

5 LaValle and Schomaker, J. Am. Chem. Soc. 61, 3520 (1939). 

6 Corey, J. Am. Chem. Soc. 60, 1598 (1938). 

7 Fajans, Ber. 53, 643 (1920); 55, 2826 (1922); Rossini, J. Research 
Nat. Bur. Stand. 13, 21 (1934); Rossini, J. Chem. Phys. 3, 438 (1935). 

8 Albrecht and Corey, J. Am. Chem. Soc. 61, 1087 (1939). 

*Fox and Martin, J. Chem. Soc. 2106 (1938); 884 (1939). 





The Dissociation of Iodine 


G. M. Murphy! calculated the general equation for 
AF°/T for the reaction I2(4)=2J(,) and compared the 
results with the experimental results of Starck and Boden- 
stein? De Vries and Rodebush* and part of those by 
Braune and Ramstetter.* For J(,) - 

P-—E, 3 5 
———— == R In M+=R In T+R In Q—7.267. (1) 
T 2 2 
Or using M=126.932, R=1.9869 cal. mole deg.~!, and 
Q=4 this reduces to 


_— 0 
_ PEt _ 11.4375 log T+9.92315. (2) 


By use of the calculated values at five temperatures, 
Murphy represented this function in the empirical form 


_P-ES 
T 


with a= —28,052, b=11.1292, c=1.907X10-4, d= —4.09 


=7+0 log T+cT+dT?+i (3) 
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X10-8, t=10.7277, whereas by the theoretical equation, 
a=0, c=0, and d=0. The values tabulated in Table II of 
his paper are too high by 0.001, 0.002, 0.001, 0.001, and 
0.0 cal. mole“ deg.—, and this is typical of the relation of 











TABLE I. 
No. No. MEAN 
INVESTI- POS. NEG. MEAN NET 
GATORS A(LOGK) Zpos.A A(LOGK) ZNEG.A DEV. DEV. 
(a) From Eq. (3) and AEo® =35,796 cal. 

B&R 9 0.489 11 1.071 0.078 —0.029 
DeV & R 3 .399 3 0.084 .071 + .043 
A 4 -221 1 .048 054 + .035 

eV&R 
S&B 7 .560 4 -132 .063 + .039 
B&R 
DeV & R 16 1.048 15 1.203 .073 — .005 
S&B 

(b) From Eq. (2) and AFo® =35,396 cal. 

B&R 1 0.134 19 2.295 0.1215 —0.108 
DeV & R 1 143 5 0.420 094 — .066 
ae: 1 012 4 .065 O15 — O11 

eV&R 
S&B 2 155 9 485 058 — .030 
B&R 
DeV & R 3 .289 28 2.491 .099 — .080 
S&B 








Eq. (3) to Eq. (2) at intermediate points. Of course, the 
solution of the five simultaneous equations, if carried out 
exactly with the correct values of —(F°—Epo°)/T, would 
lead to the constants of Eq. (2). That is, values of a, c, and 
d would result, which are of such small magnitudes, that 
when multiplied by the 1/7, T, and T? terms, the contribu- 
tions to — (F°— E,°)/T are negligible, and they are properly 
set equal to zero. 

Since one of the purposes of the paper was the com- 
parison of the experimental equilibrium constants with 
those determined from the spectroscopic data, a more 
fundamental objection may be made to the use of AEp° 
(=35,795+100 cal.) derived from the experimental data. 
The spectroscopic value by Brown,’ D=1.535+0.001 ev. 
or AEo®=35,395+23 cal. appears to be a satisfactory 
determination. The consistency of the experimental results 
decreases in the order (1) Starck and Bodenstein (S & B), 
(2) De Vries and Rodebush (De V & R), (3) Braune and 
Ramstetter (B & R). Table I shows that the deviations of 
the experimental values of log K by Starck and Bodenstein 
from log K derived from Eq. (2) and AEo®=35,395 cal. 
are less than a third of the corresponding deviations from 
log K derived from Eq. (3) and AE»®=35,795 cal. The 
deviations of the results of the first two groups of investi- 
gators from the data of purely spectroscopic origin are 
also less than those from the hybrid equation involving 
AE)? =35,795 cal. 

HucH M. SPENCER 


Cobb Chemical Laboratory, 
University of Virginia, 
Charlottesville, Virginia, 
May 4, 1940. 
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Transition in Nickel Nitrate Hexammoniate 


The low temperature crystal structure of nickel nitrate 
hexammoniate has been investigated by Yii,! who found 
evidences of very high thermal agitation of the nitrate 
groups in the crystal. We have measured the heat capacity 
curve for this substance from 54°K to 300°K, using a 
modified Nernst vacuum calorimeter. The data, shown in 
Fig. 1, indicate the existence of two regions of high energy 
absorption, one a typical lambda-type transition extending 
from 173°K to 247°K, with a maximum value at 243°K, 
and the other a region of anomalous energy absorption 
below 85°K. 

The lambda transition can be supercooled to liquid-air 
temperatures and held there indefinitely without trans- 
forming; the rate of transformation becomes appreciable 
only upon warming to about 200°K. Cooling-curve meas- 
urements indicate a hysteresis effect of several degrees. In 
one cooling experiment, the minimum cooling rate occurred 
at 239°K. 

The height of the lambda-peak at 243°K is not deter- 
minable; our highest C, measurement in this region was 
1270 calories/mole/degree, certainly a lower limit. The 
total heat absorption in excess of that represented by the 
“‘normal’’ heat capacity curve was obtained by “heat of 
transition’? measurements extending from 227—247°K, thus 
including the area which could not be integrated graph- 
ically. The excess heat effect (above the “normal” curve 
from 173—247°K) is 1818 calories/mole. The corresponding 
entropy is 7.65+0.5 E.U. 

This high entropy of transition seems to invalidate any 
explanation involving rotation of only the nitrate groups 
in the crystal lattice. The assumption that both the 
ammonia and nitrate groups are executing rotational oscil- 
lations gives fair agreement with the experimental entropy, 
and with the value of C, just above the transition. This 
assumption is in accordance with the theory of Frenkel,? 
and that of Kirkwood,? which assume that the active 
groups are librating both below and above the transition, 
and that the transition is simply a process of disorientation, 
or the disordering of these librations with respect to the 
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Fic. 1. Heat capacity in calories per degree per mole of nickel nitrate 
hexammoniate. 
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crystal axes. A more precise check of the assumption 
could probably be made from an investigation of C, values 
in the transition region.‘ If the ammonia groups are librat- 
ing, as proposed, the prsence of nitrate or similar groups 
would seem to be required, since cooling-curve measure- 
ments on nickel chloride hexammoniate® show no indication 
of a transition. We are investigating ammonium chloro- 
stannate and other substances of similar crystal structure 
to the nickel hexammoniates, in order to establish more 
definitely the conditions for this phenomenon. 

The region of anomalous energy absorption below 85°K 
is of much smaller magnitude, the C, values being about 
five percent high. The dotted line in Fig. 1 shows the results 
of measurements made on a sample which had been par- 
tially supercooled through the 243° transition. The pres- 
ence of the form stable above 243° almost completely 
eliminates the high energy absorption below 85°. Measure- 
ments at temperatures considerably below 54°K will be 
necessary for a full description of the phenomenon. 

A detailed report of the work will be published later. 
This investigation has been supported by a grant from 
the University of Missouri Research Council. 

Ear A. LonG 


FreD C. TOETTCHER* 
Department of Chemistry, 
University of Missouri, 
Columbia, Missouri, 
May 3, 1940. 
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The Mercury Photosensitized Hydrogenation of Ethylene, 
Propylene, and n-Butylene 


A convenient method of preparing a free aliphatic 
radical is by the addition of a hydrogen atom to the corre- 
sponding olefine; if the hydrogen atoms are produced by 
mercury photosensitization their concentration will be 
low compared with the radical concentration, and sub- 
sequent reactions of the radicals with each other may be 
studied. Thus Taylor and Jungers! reported that butane 
was the predominant product from the mercury photo- 
sensitized hydrogenation of ethylene. 

In an apparatus similar to that of Taylor and Jungers 
we have compared the reactions of ethyl, propyl, and buty! 
radicals formed in a similar fashion. In each case olefine 
at about forty millimeters pressure, mixed with a sixfold 
excess of hydrogen, was illuminated at room temperature 
in the presence of mercury vapor with 2536.7A radiation. 
The analytical method employed was a suitable adaptation 


TABLE I. 








PERCENTAGE OF OLEFINE TO SATURATED PRODUCTS 





* OLEFINE 
USED CHa 


CoHa 2 
C3He 2 1 
CsHs 3 1 


CeHe | C3Hs | CaHio 


14 _ 84 
26 5 
1 29 


CsHi2 | CoHus | CsHis 
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of that previously described.? The results are given in 
Table I as the means of two concordant experiments with 
each olefine. In some of the experiments a trap cooled to 
— 78° was used to decrease if possible the amount of lower 
hydrocarbon products formed, by preventing more exten- 
sive secondary reactions of higher hydrocarbons with 
hydrogen atoms. 

There is a significantly greater formation of propane 
and butane from propylene and butylene than of ethane 
from ethylene. A possible explanation may be that addition 
of a hydrogen atom to propylene and n-butylene may form 
the isopropyl and isobutyl radicals as well as the normal 
radicals. Our earlier work with mercury diethyl indicates 
that the disproportionation reaction 2C,Hoen41=CnzHon+2 
+C,He, is not important with the normal radicals, C2Hs. 
The reaction may become of greater importance when iso- 
radicals are involved. It is of interest in this connection to 
compare the observation of Norrish* that considerable 
amounts of propane and propylene are formed from the 
photolysis of di-isopropyl ketone, while neither of these 
products is obtained from di-n-propyl ketone. 

WALTER J. Moore, Jr. 


HuGu S. TAYLor 
Frick Chemical Laboratory, 
Princeton University, 
Princeton, New Jersey, 
May 10, 1940. 


1H. S. Taylor and J. C. Jungers, J. Chem. Phys. 6, 325 (1938). 
2W. J. Moore, Jr., and H. S. Taylor, J. Chem. Phys. 8, 396 (1940). 


3R. G. W. Norrish, quoted in Pearson and Glazebrook, J. Chem. 
Soc. 1777 (1936). 





The Thermodynamic Properties of the PN Molecule* 


In 1938 Moureu and Wetroff! showed that equi- 
molecular quantities of phosphorus and nitrogen, when 
heated together, form a compound (PN). In a more 
recent paper Moureu and his co-workers? show the exist- 
ence of an equilibrium between (PN), and PN, and also 
one between PN and P2+Nbe. In view of their results, we 
thought that it might be of interest to calculate from 
spectroscopic data the thermodynamic properties of a 
perfect gas composed of PN molecules, and the equilibrium 
constant for the gaseous reaction 


By combining the values for the equilibrium constant of 
reaction (1) with the values of the equilibrium constant 
for the reaction 


P,=2P2 (2) 
given by Dushman,? the equilibrium constant for 
Ps +2N2=4PN (3) 


can be calculated. 


The results of these calculations are listed in Tables I 
and IT. 


In the calculation of the thermodynamic functions, the 
total partition function Q is given by 
Q - Qtrans ss Qinternal ‘ Qcorrcetion: 


Qtrans has been calculated by the usual method assuming 
a perfect gas. Qinternal assumes an harmonic oscillator 
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TABLE I. The thermodynamic functions of the PN molecule. (One mole of 
the perfect gas at one atmosphere.) 











—(F —Fo) (H —Ho) Ss Cp 
Temp. °K CAL. CAL. CAL./DEGREE| CAL./DEGREE 
100.0 3588.0 694.6 42.824 6.9544 
298.1 12959.8 2079.1 50.448 7.0973 
500.0 23553.6 3565.1 54.235 7.6474 
750.0 37538.7 5546.6 57.450 8.1759 
1000.0 52209.0 7633.8 59.842 8.4774 




















TABLE II. The equilibrium constant for P2+N2=2PN. 











Temp. °C Kye t4E0/RT e— 4Eo/RT Kp 
800 6.09 7.87 47.9 
900 5.94 6.58 39.1 

1000 5.82 5.68 33.1 











and a rigid rotator. From the value of Qgorrection due to 
anharmonicity in the oscillation and to the lack of rigidity 
in the rotation given by Mayer and Mayer,’ we have 
calculated the following formulae for the corrections in 
the thermodynamic functions. The symbols used are de- 
fined in reference (4). 


RT 
Foorrestion - -—|8y+ 





bu 2xu? 
ev—1 at 

167 6 4xu*e“ 
re Pat aie 


8y u(de"—2x) 4xuter 
A ecorrection =rr|— | 


due” 





Scorrection _ R{ 








(ev—1)? (ew—1)8 
l6y ude" 
Ceorrection =R|— ~ (e«—1)2 
2ute"(de"—4x—2ux) 1 =) 
(ex—1)3 (ex—1) 


The molecular constants which we have used were taken 
from the literature.6 From the given values of the heats of 
dissociation, AE for reaction (1) is calculated to be —0.2 
electron volt. 

According to the work of Moureu and his co-workers,? 
PN at 800°C is almost completely dissociated into P2 and 
No. The difference between this work and the values we 
have obtained may be attributed to the fact that the Eo 
for PN has been determined by an extrapolation method,” 
and is probably an upper limit. A change in Eo for PN of 
only a few tenths of an electron volt would have the effect 
of lowering the calculated value of K,. 


KENNETH J. McCALLUM 


EpGAR LEIFER** 
Havemeyer Laboratories, 
Columbia University, 
New York, New York, 
May 13, 1940. 
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506 LETTERS TO 
The Entropy of Dimethyl Sulfide from Low Temperature 
Calorimetric Measurements. Restricted Ro- 
tation of the Methyl Groups 


In order to arrive at a satisfactory theory for the poten- 
tial restricting the rotation of methyl groups in many 
molecules, it is desirable to determine how the magnitude 
of the barrier depends on the kind of atom to which the 
methyl groups are bonded. Recently we have obtained an 
estimate of the barriers in dimethyl sulfide by comparing 
the entropy obtained from calorimetric measurements ex- 
tending to low temperatures with that computed from 
molecular data. 

The heat capacities of dimethyl] sulfide in the tempera- 
ture range 13.4 to 291°K, heat of fusion, and heat of 
vaporization at 291°K were accurately determined in an 
adiabatic calorimeter equipped with a calibrated platinum 
resistance thermometer. In addition, precise measurements 
of the vapor pressure were made from 250 to 293°K. The 
size of the sample in the calorimeter was 1.4432 moles, 
with an estimated impurity of 0.007 mole percent. Com- 
plete details will be published later. 

In Table I are given the entropy changes calculated 
from our calorimetric data. The final temperature and 
pressure are those for the directly measured heat of 
vaporization, 6696 cal./mole, which may be compared 
with 6690 cal./mole calculated from the Clapeyron relation, 
assuming that the gas follows the Berthelot equation of 
state, with critical temperature 503°K and critical pressure 
54.6 atmos.' This same assumption was made in the correc- 
tion for gas imperfection. 

The calculation of the entropy from molecular data is 
summarized in Table II. The translational entropy comes 
from the familiar Sackur-Tetrode equation. The following 
frequencies and multiplicities,? as well as the usual approxi- 
mation of harmonic oscillation, were used to obtain the 
vibrational entropy: 284, 695, 742, 1030 (2), 1262 (2), 
1323 (2), 1445 (4), and 2900 (6) cm™. The calculation of 
the rotational entropy was made by means of the formulas 
of Crawford,* in the derivation of which it is'assumed that 
the potential energy of each methyl group has tye form 


TABLE I. Molal entropy of dimethyl sulfide from calorimetric data. 








0-15°K, Debye extrapolation, hv/k =110.4 
15-174.855°K, graphical 

Fusion, 1908.4/174.855 

174.855-—291.06°K, graphical 

Vaporization, 6696/291.06 


0.369 
21.156 
10.914 
13.821 
23.006 


Entropy of actual gas, 291.06°K, 365.5 mm 69.27 +0.1 
Correction for gas imperfection 0.08 


Entropy of ideal gas, 291.06°K, 365.5 mm 69.35 cal./deg. mole 
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TABLE II. Molal entropy of gaseous dimethyl sulfide from molecular data 
at 291.06°K and 365.5 mm. 








Strans 
Svib 


Srot 


(Harmonic) 
(V =2015 cal./mole) 


Total entropy (V =2015) 
(V =1800) 


69.35 cal./deg. mole 
69.60 


71.03 
64.08 


(V =0) 
(v= 








(V/2)(1—cos 3a), where a@ is the angular displacement of 
the methyl group from its equilibrium position. The bond 
angles and interatomic distances taken were <C—S—C 
=102°, C—S=1.82A,A C—H=1.09A, and tetrahedral 
angles for the methyl groups. These values lead to principal 
moments of inertia of 45.49, 116.1, and 151.0 10~* g cm? 
and to a moment of inertia 5.30310-* and reduced 
moment of inertia 4.818X10-* for each methyl group. 
Since the ratio of these last two is nearly unity, use of 
Pitzer’s tables does not cause great error at the tempera- 
ture of interest. Thus the entropy for V=2015 cal./mole 
calculated with the aid of Pitzer’s tables is 69.28 cal./deg. 
mole rather than 69.35 cal./deg. mole by the more accurate 
method of Crawford. 

It is seen that if the values assumed for the structural 
parameters and for the vibrational frequencies are suffi- 
ciently accurate, a potential barrier of about 2000 cal./mole 
for each methyl group is required for dimethyl sulfide. 
The fact that the barriers in ethane® and in the methyl 
amines’ are of about this same magnitude, while that in 
dimethyl acetylene? is negligible, suggests that the restric- 
tion is mainly due to the proximity of the methyl groups 
rather than to the kind of atom to which the methyl groups 
are attached. 

We wish to express our thanks to Professor Alexander 
Goetz for supplying us with liquid hydrogen. 

DarRRELL W. OsBORNE* 
RussELt N. DoESCHER 


Don M. Yost 
Gates Chemical Laboratory, 
California Institute of Technology, 
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May 11, 1940. 
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